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SUMMARY 

A  system  ot  equations  has  been  developed  tor  t  In'  simulation  ot  an  airerut t's 
motion  in  real  time  using  a  digital  eomputer.  Those  parts  ot  the  maihematieal 
model  common  to  all  aircraft  have  been  created  as  a  set  o!  Fortran  subroutines, 
leaving  the  user  to  create  only  a  small  group  ot  routines  specifically  to  describe 
his  aircraft.  The  equations  employed  are  defined  and  the  computet  implement  at  ion 
described  in  detail.  The  Report  can  be  used  as  a  handbook  and  'user  guide'  hut 
as  the  routines  described  are  not  specific  to  real-time  simulation  thev  could  he 
used  as  a  basis  for  a  general  mathematical  model  of  an  aircraft  for  use  on  anv 
computer  which  supports  Fortran. 
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INTRODUCTION 


This  Report  describes  a  system  of  equations  for  the  simulation  of  an  air¬ 
craft  ' s  Vigid-body'  mot  ion  in  real  time  using  a  digital  computer.  Known  as  SESAME 
a  System  of  Equations  for  the  Simulation  of  Aircraft  in  a  Modular  Environment  - 
it  may  be  described  as  opening  the  door  to  flexible  use  of  the  simulator.  Those 
parts  of  the  mathematical  model  which  are  common  to  all  aircraft,  such  as  the 
equations  of  motion  and  axis  transformations,  have  been  created  as  a  set  of 
standard  modules  (written  in  Fortran),  leaving  the  user  to  create  only  a  small 
group  of  routines  specifically  to  describe  his  aircraft.  The  two  sets  of  modules 
are  then  linked  together  to  produce  a  complete  model  program. 

The  simulator  for  which  SESAME  has  been  created  changed  from  pure  analogue 
to  hybrid  computing*  in  1974.  The  computer  system  consists  of  a  Xerox  Sigma  S 
digital  computer  linked  to  an  Applied  Dynamics  AD4  analogue  computer.  This  Report 
is  concerned  mainly  with  the  digital  computer  aspects.  In  the  early  days  of 
using  this  new  computing  facility,  the  aircraft's  mathematical  model,  computed 
entirely  digitally,  was  coded  using  SI. I,  a  simulation  language  for  solving  sets 
of  ordinary  differential  equations.  With  increasing  experience  and  increasingly 
complex  simulations,  the  single  large  program  imposed  by  the  SLI  language  was 
found  to  be  more  and  more  inconvenient  and  the  decision  was  taken  t  restructure 
the  program  into  small,  self-contained  modules. 

Such  restructuring  became  feasible  largely  because  a  newly-created  software 
package,  for  on-line  parameter  variation  and  inspection,  provided  the  means  to 
access  and  change,  on-line,  variables  contained  in  subrout inos . 

Although  the  basic  modules  are  themselves  coded  in  Fortran,  SLI  itself  has 
been  retained  to  provide  the  overall  framework  in  which  the  modules  sit.  This  is 
because  SLI  provides  facilities  for  numerical  integration  and  for  synchroni sat  ion 
with  real-time. 

The  structure  of  the  Report  is  as  follows.  After  a  discussion,  in  section 
of  the  objectives  to  be  achieved  by  the  modular  nature  of  SESAME,  section  1  gives 
a  brief  theoretical  background  to  the  system  equations,  including  such  topics  as 
choice  of  axes  and  definition  of  the  equations  of  motion.  This  is  followed  in 
section  4  by  an  outline  of  how  the  equations  to  be  solved  are  distributed  among 
the  specific  routines.  This  is  supplemented  (in  Appendix  A)  by  a  detailed  defini¬ 
tion  of  each  routine  and,  as  the  ultimate  specification  for  reference  purposes, 
by  a  listing  of  the  routines  (Appendix  B) .  Section  5  then  describes  the  overall 
program  structure  in  which  the  routines  are  embedded.  Section  h  discusses  how 
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communication  among  routines  is  achieved  using  the  Fortran  named  COMMON  facility, 
and  how  the  user  should  adopt  the  same  technique  for  his  routines.  Hybrid  input/ 
output  (analogue  to  digital  and  digital  to  analogue  conversion,  and  single  bit 
discretes)  is  also  described  in  section  6  in  terms  of  what  the  user  receives  and 
must  provide,  while  section  7  is  a  general  "User's  Guide",  summarising  what  the 
system  expects  from  the  user,  defining  the  set  of  routines  and  the  data  the  user 
must  provide,  and  outlining  the  choices  the  user  may  make  regarding  the  atmos¬ 
phere,  turbulence  and  other  features. 

The  Report  has  been  written  with  three  principal  aims:  to  provide  primary 
simulator  users  with  a  description  of  the  standard  computer  model  and  an  outline 
of  how  to  build  it  into  a  full  simulation;  to  provide  simulator  customers  with  an 
idea  of  how  the  model  is  formulated  and  what  must  be  done  to  represent  an  indivi¬ 
dual  aircraft,  so  that  they  can,  if  necessary,  contribute  to  the  creation  of  a 
complete  simulation  model;  and  to  provide  a  general  description  of  a  computer 
model  for  others  interested  in  doing  something  similar. 

2  OBJECTIVES 

Specific  objectives  to  be  achieved  by  the  modular  system  of  equations  con¬ 
cern  size,  testing,  transferability,  communication  and  execution  time,  and  are 
intended  to  overcome  detailed  weaknesses  of  the  simulation  language  SLI.  However, 
as  explained  later  in  section  5.1,  SLI  is  retained  as  the  appropriate  language 
Tor  the  overall  program,  containing  the  individual  modules. 

2 . 1  Size 

A  complete  model  program  is  large,  often  exceeding  2000C  words  of  computer 
storage.  If  coded  entirely  in  SLI  source  language  this  means  many  lines  of  code. 
The  process  of  translating  the  SLI  source  code  into  Fortran  and  then  compiling  is 
slow  and  tedious  and  must  be  repeated  in  its  entirety  whenever  a  change  is  made 
at  source  level.  The  present  scheme  drastically  cuts  the  size  of  the  SLI  source 
program,  which  in  itself  will  speed  up  the  translation/compilation  process.  In 
addition,  changes  at  SLI  level  will  be  rare,  since  the  detail  of  a  simulation  is 
relegated  to  a  lower  level  of  routine.  The  existence  of  subroutines  also  provides 
the  potential  ability  to  create  overlays,  in  order  to  reduce  main  store  occupancy. 

2.2  Creation  and  testing 

Creation  and  checking  of  modules  will  be  easier  than  handling  a  large  SLI 
program  since  routines,  such  as  for  the  generation  of  the  aircraft's  aerodynamic 
forces,  can  be  defined,  created  and  tested  in  isolation  before  being  included  in 
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the  main  program.  Work  could  proceed  in  parallel  on  several  routines,  and  be 
done  by  external  users  of  the  simulator,  ie  people  who  are  not  specialists  in 
simulation  as  such. 

Modifications  will  also  be  easier,  since  they  should  involve  only  revision 
of  a  single  subroutine,  its  compilation  and  then  creation  of  a  new  load  module 
without  having  to  repeat  the  SL1  translation. 

Using  Fortran  rather  than  SLI  as  the  primary  medium  does  not  mean  that  all 
the  useful  SLI  facilities,  such  as  function  generation  by  table  look-up,  must  be 
discarded.  These  can  still  be  used  at  the  Fortran  level,  since  each  one  involves 
a  Fortran-callable  subroutine  or  function.  However,  the  calling  sequence  will, 
in  general,  be  more  complicated  in  Fortran  than  in  SLI.  SLI  operators  involving 
integration  are  not  usable  in  subroutines  without  considerable  effort. 

2.3  Transferability 

Having  a  modular  structure  means  that  simulation  of  a  new  aircraft  requires 
only  the  creation  of  the  routines  specific  to  that  aircraft.  The  standard  parts 
can  be  picked  up  from  libraries  and  special  facilities,  such  as  simulation  of 
guidance  beams,  can  readily  be  moved  from  simulation  to  simulation  in  subroutine 
form.  With  SLI  the  only  way  of  transferring  across  an  existing  piece  of  code 
is  by  explicit  inclusion  of  the  source  lines  and  punching  of  the  necessary  cards. 

None  of  the  modules  includes  any  compromising  feature  which  limits  its 
employment  to  real-time  simulation.  Hence  the  routines  described  could  also  form 
the  basis  for  a  conventional  digital  model  of  aircraft  dynamics  for  use  in  non- 
real-time  studies. 

2.4  Communication 

A  defect  of  the  SLI  language  is  the  way  most  of  the  significant  variables 
are  lumped  into  one  large  block  of  labelled  COMMON,  the  name  of  which  is  not  fixed 
and  the  contents  of  which  can  change  if  the  SLI  source  code  is  varied. 

In  this  new  system  communication  of  variables  among  all  the  modules  is 
achieved  by  using  Fortran  labelled  COMMON,  but  under  user  control.  This  retains 
the  advantage  of  using  COMMON  storage,  that  variables  are  conveniently  accessible 
by  any  subroutine  without  using  arguments,  thus  saving  execution  time. 

2 . 5  Execution  time 

The  aircraft's  motion  is  computed  by  repeatedly  integrating  the  differen¬ 
tial  equations.  If  a  multi-pass  integration  technique  is  used,  such  as  Runge- 
Kutta  fourth  order,  there  are  some  sections  of  the  equations  which  are  not  part 
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of  the  dynamic  loop  and  so  need  only  be  executed  once  per  time  step.  Examples 
are  calculation  of  ILS  guidance  and  TV  position  signals.  By  having  these  elements 
as  subroutines,  it  is  easy  to  arrange  that  they  are  executed  only  as  often  as 
necessary  and  so  save  time  for  more  important  functions. 

To  date,  aircraft  model  programs  have  been  solved  with  only  one  basic  loop 
or  frame.  SLI  does  provide  a  means  to  execute  one  part  at  a  different  iteration 
rate  from  another.  Should  circumstances  arise  where  more  than  one  frame  time  is 
necessary,  recasting  of  the  model  program  (at  SLI  skeleton  level)  will  be  aided 
by  the  simplicity  of  merely  moving  subroutine  calls  from  one  'derivative'  section 
to  another,  rather  than  blocks  of  code  which  would  be  necessary  in  a  program 
coded  completely  in  SLI.  Communication  between  'derivative'  sections  is  also  no 
problem  when  variables  are  in  COMMON  areas  controlled  by  the  user  rather  than  by 
the  SLI  translator. 

2.6  Data  logging 

Access  to  variables  for  data  logging  purposes  is  easy  and  convenient.  A 
self-contained  data  logging  routine  can  be  written  knowing  that  all  variables 
may  be  accessed  without  special  action  and  without  run-time  inefficiencies. 

3  SYSTEM  EQUATIONS 

3. I  Choice  of  axes 

The  first  requirement  is  to  choose  sets  of  axes  in  which  to  solve  the 
.  .  .  2-4 

fundamental  equations  of  motion.  In  classic  texts  a  set  of  axes  fixed  in  the 
aircraft  is  generally  chosen  but,  being  a  rotating  frame  of  reference,  such  body 
axes  result  in  the  translational  accelerations  including  angular  velocity  terms, 
as  in 


X  -  mg  sin  6  =  m(u  +  qw  -  rv)  . 

For  efficient  computer  solution  of  the  equations  of  motion,  it  is  desirable  to 
uncouple  translational  motion  from  rotational  motion,  a  point  which  was  made  many 
years  ago  by  Howe"*  and  reiterated  recently^.  This  may  be  achieved  by  suitable 
choice  of  axes.  Then,  since  rotational  motion  of  aircraft  is  intrinsically  more 
rapid  than  translational  motion, solut ion  of  the  rotational  equations  may  be 
performed  more  frequently  than  the  translational  equations  without  an  excessive 
computing  load,  and  while  maintaining  overall  accuracy. 

Body-fixed  axes  still  provide  a  natural  frame  for  the  solution  of  rotational 
equations  of  motion,  with  the  advantage  of  constant  moments  of  inertia.  Choice 
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of  axes  therefore  reduces  to  selecting  an  appropriate  frame  for  the  translational 
equations  of  motion. 

So-called  'flight  path  axes'  could  be  a  suitable  choice,  with  the  origin 
at  the  aircraft's  centre  of  gravity  and  the  x  axis  aligned  with  the  aircraft's 
velocity  vector  with  respect  to  the  ground.  However,  in  the  presence  of  winds 
and  turbulence  this  introduces  complications,  so  that,  following  Ref  7,  the 
choice  falls  on  earth-based  axes. 

3 .  2  L)e  f  i  n  i  t  i  on  o  f  axe  s 

All  axes  systems  used  are  orthogonal,  right-handed  triads. 

3.2.1  Earth  axes 

Earth  axes  are  an  inertial  frame  assuming  a  flat,  non-rotating  earth. 

(See  McFarland  for  the  case  of  a  spherical,  rotating  earth.) 

The  origin  is  at  a  datum  point  on  the  visual  model*  in  use,  typically  at 
the  runway  threshold  and  on  the  centreline.  The  x-axis  points  northward 
(suffix  N) ,  the  y-axis  points  eastward  (suffix  E),  the  xy  plane  being  parallel 
to  the  earth's  surface,  and  the  z-axis  (suffix  H)  points  down  to  the  centre  of 
the  earth. 

3.2.2  Geometric  body  axes 

Geometric  body  axes  have  their  origin  at  the  aircraft's  centre  of  gravity 
and  are  located  with  respect  to  the  aircraft  by  some  geometric  feature  such  as 
the  longitudinal  fuselage  datum  line.  Once  defined  they  are  fixed  in  the  aircraft 
The  x-axis  points  forward,  the  y-axis  to  starboard  and  the  z-axis  'down'. 

3.3  Aircraft  attitude 

£ 

Aircraft  angular  orientation  with  respect  to  the  earth  is  defined  by  a 
conventional  trio  of  Euler  angles  ip,  9,  .  The  heading  angle  v  is  measured 

from  north  and  lies  in  the  range  0,  360°;  the  elevation  angle  or  pitch  attitude 

,  .  ,  o 

t)  is  measured  from  the  horizontal  plane  and  lies  in  the  range  -90,  +90  ;  the 

roll  (or  bank)  angle  is  measured  from  the  horizontal  plane  and  lies  in  the  range 

-180,  +180°. 

Rates  of  change  of  these  angles  may  be  related  to  the  components  of  angular 
velocity  of  the  aircraft  by 

*  Outside  world  cues  are  provided  by  a  closed-circuit  TV  system  viewing  a 
physical  model  of  an  appropriate  terrain. 
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0  «  q  cos  $  -  r  sin  $ 

♦  *  p  *  (q  sin  $  ♦  r  cos  $)  tan  0 
4>  ■  (q  sin  $  ♦  r  cos  $)  sec  0 


p  ♦  v  sin  0 


and  the  angles  themselves  obtained  by  integration.  However,  a  singularity  occurs 
o  .  ... 

at  0-90  .  If  all-attitude  manoeuvring  is  desired,  an  alternative  formulation 

of  equations  is  necessary.  No  provision  has  yet  been  made  for  multiple  rolls 
or  turns. 

3.4  Transformation  from  earth  to  body  axes 

Transformation  of  a  set  of  variables  from  earth  axes  to  body  axes  (or  vice 

.  .  4  9 

versa)  is  most  conveniently  achieved  through  the  direction  cosine  matrix  ’  . 

For  example,  the  components  of  airspeed  in  body  axes  (u„,  v_,  wD)  are  related 

Boo 

to  the  components  in  earth  axes  (V^,  V^,  V^)  by 


UB  "  Vs  +  J2VE  +  >'3VD 


vB  -  m)VN  ♦  m2V£  ♦  m3VD 


w  -  n  V  +  n  V  +  n  V 
B  IN  2  E  3D 


where  tj  etc  are  the  direction  cosines,  given  by 


1 1  -  S  |  -  cos  0  cos  \p 


l,  -  Sp  -  cos  0  sin  \1> 


S  "  s!3  “  ‘  sin  9 


m^  -  S,  j  -  sin  $  sin  0  cos  ip  -  cos 


m.,  -  S,,  -  sin  $  sin  0  sin  +  cos 


♦  sin  '+> 
$  cos  ip 


ra3  "  ^">3  "  s^n  $  cos  0 


-  S3(  -  cos  P  sin  9  cos  <p  *  sin 


i  sm  ; 


n,  -  S-»  -  cos  sin  6  sin  £  -  sin 


$  COS  Ip 


n^  -  S^3  -  cos  P  cos  0  , 


In  matrix  form,  equations  (2)  may  be  written 


{M  -  s{v4  “> 

where  {ub|  *  {Vn|  are  co^umn  matrices  (after  Hopkin*)  and  S  ,  the  transform¬ 
ation  matrix,  is  orthogonal,  its  transpose  ST  being  the  same  as  its  inverse. 

The  inverse  relationship,  for  earth  axes  variables  in  terms  of  body  axis 
variables,  is  then  given  by 


Thus,  for  example,  the  total  forces  (suffix  T)  applied  in  body  axes  (suffix  X,  Y, 
Z)  may  be  transformed  to  earth  axes  (suffix  N,  E,  D)  by 


{ftn|  “  s  {Ftx| 


which,  expanded,  is 


FTN  “  S1 1FTX  4  S21FTY  +  S3lFTZ 

FTE  "  S12FTX  4  S22FTY  4  S32FTZ  k 

FTD  "  S|3FTX  4  S23FTY  4  S33FTZ  * 

Other  vector  components  are  related  in  the  same  way. 

3.5  Equations  of  motion  -  translation 

The  components  of  acceleration  with  respect  to  the  earth  are  obtained  from 


v 

KN 

"  FTN/m 

C’ke 

"  FTE/m 

-  FTD/m  ♦  g 

where  the  acceleration  due  to  gravity,  g  ,  is  assumed  constant,  and  the  forces 
etc  in  earth  axes  are  obtained  from  the  total  forces  in  body  axes 

ftx*  Fty»  Ftz  by  the  e<luations  <7)* 

Integration  of  equations  (8)  yields  the  velocity  components  (V^,  VK£»  vkd^ 
of  the  vehicle  relative  to  the  earth.  Relative  to  the  air  mass,  the  velocity 
components  of  the  aircraft  are 
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V  -  V  -  v 
N  KN  WN 


V  -  V 
KE  WE 


V  -  V 
KD  WD 


where  V^,  V^,,  are  the  components  of  the  wind  velocity  relative  to  the 

ground. 

The  body-axes  components  of  the  velocity  vector  relative  to  the  air  may 
then  be  derived  using  equations  (2)  and  (9),  from  which  the  velocity  (or  True 
Airspeed,  TAS)  is 

/  2  2  2  V 

VT  =  (UB  *  VB  +  WB /  0 

and  equivalent  airspeed  is  derived  using  the  atmospheric  density  ratio 


Angles  of  attack  and  sideslip  are  computed  from 

a  =  tan  '  wfi/uB  (12] 

(a  being  in  the  range  -180,  +180°)  and 

3  =  tan"'  VB/(UB  +  ^  (,3: 

(3  being  in  the  range  -90,  +90  and  taking  the  sign  of  vfi  ) . 

Derivatives  of  these  angles  with  respect  to  time,  needed  to  calculate  the 
aerodynamic  forces  and  moments,  are  created  by  simple  difference  equations. 
Flight  path  angles  defining  climb  (y)  and  track  (x)  are  derived  from  the 
velocity  vector  components  relative  to  the  ground,  having  first  obtained  the 
ground  speed  (VK>  as  the  horizontal  component 

VK  ■  (VL*Vk)!  <141 

'  '  VJ\)  <,5) 

*  '  ta"”'(VKE/V W.)  •  (l6) 


Positional  coordinates  x,  y,  h  of  the  aircraft's  centre  of  gravity  are 
found  by  integrating  the  velocities 
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*  =  VKN 


h  =  -V 


"  ■  'KD  •  J 

A  block  diagram  of  the  translation  equations  is  illustrated  in  Fig  1. 

3 . 6  Equations  of  motion  -  rotation 

It  is  assumed  that  the  x  and  z  body  axes  lie  in  a  plane  of  mass  symmetry 

so  that  products  of  inertia  I  and  1  are  zero.  Then  the  classical  equa- 

4  yz  xy 

tions  (Hopkin  ,  section  10.1)  may  be  manipulated  to  give 


(  *  r1 

5  '  L*T7“'[v 


(I  +  I  ~  I  )  pq 

Z  X  y 


+  (1y"V 


I2  -| 

zx 

*z  J 


*?  Vx  ' 


q  =  |m  ♦  !zx(r2  -  p2)  +  (Iz  -  Vrp \  ly 


<  N  +  L  +  £ 


dx  -  V  +  pq 


z  -TT 


p  =  CIjL  +  CI9N  +  (CI3P  +  CI4r)q 

q  =  CI5M  +  CI6(r2  -  p2)  +  CI7rp 

f  =  CIgN  +  CI2L  +  (Cl9p  +  CI|()r)q 


where  Cl j  etc  are  constants  evaluated  during  initialisation. 
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Cl 

Cl 

Cl 

Cl 

Cl 

Cl 

Cl 


t 

2 

3 

4 

5 

6 

7 

8 

9 
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1  /(l  1  -  l2  )  \ 

zl\  x  z  zx  / 

1  l( II  -  I2  ) 

ZX/y  X  z  zx  J 

‘..(‘x  '  '»  *  ‘<)/(V'  -  ‘.x) 

fl  (I  -  I  )  -  I2  )/(l  I  -  I2  \ 
l  z  y  z  zx/J  y  x  z  zx  J 

l/!y 

•xx/S 

(*.  -  ‘»)/v 

■„/(v.  -  ‘L) 

(i  (I  -  I  >  «■  I2  )j(l  1  -  I2  ) 

V  x  x  y  zx  1 1 V  x  z  zx/ 

-I  (l  -  I  ♦  I  V(l  I  -  I2  )  . 

zxyx  y  z  y  X  z  zx/  y 


(20) 


No  terms  are  included  in  these  equations  (18)  and  (19)  to  allow  for  engine 
gyroscopic  effects,  but  they  could  easily  be  added  by  the  user  when  he  provides 
the  total  moments  L,  M,  N  .  It  would  also  be  straightforward  to  extend  the 
equations  to  allow  fora  non-symmetric  mass  distribution. 

Solution  of  the  equations  of  rotational  motion  therefore  consists  in  cal¬ 
culating  the  angular  acceleration  components  from  equations  (19)  given  the  total 
moments  L,  M,  N  {  integrating  to  obtain  angular  velocity  components;  transforming, 
using  equations  (I),  to  attitude  rates;  integrating  again  to  attitude  angles  and 
calculating  the  direction  cosines  according  to  equations  (3).  This  total  process 
is  illustrated  in  Fig  2. 

3.7  Centre  of  gravity  location 

An  aircraft’s  centre  of  gravity  location  is  usually  defined,  in  a  fore-and- 
aft  sense,  in  terms  of  a  reference  length  such  as  the  mean  chord  c  .  Thus  if 
the  centre  of  gravity  is  quoted  as  being  at  0.1c,  then  it  is  10X  of  the  reference 
chord  af't  of  the  origin.  This  convention  is  in  a  direction  contrary  to  the  normal 
positive  sense  of  x,  but  will  be  retained  here.  The  vertical  location  of  the 
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centre  of  gravity,  however,  is  defined  as  being  measured  in  the  positive  z 
direction  ( ie  downwards)  from  some  origin. 

Aerodynamic  moment  data  may  often  be  quoted  relative  to  some  reference 

centre  of  gravity  position,  so  that  for  an  actual  operational  centre  of  gravity 

position,  some  correction  is  necessary.  This  is  provided  for  as  follows.  If 

the  moment  reference  position  of  the  centre  of  gravity  is  x  , ,  z  ..  and 

r  cgref  cgref 

the  actual  centre  of  gravity  position  is  x  ,  z  ,  all  defined  as  fractions 

&  v  cg»  cg 

of  some  reference  length  c  ,  then  the  differences  are  calculated  as 

°  ref 


Ax  «  x  -  x  , 
cg  cg  cgref 


Az  -  z  -  z  , 
eg  cg  cgref  . 


Then  differences  are  then  available  to  correct  aerodynamic  moments  to  the  current 
centre  of  gravity  position. 

3.8  Pilot  position  and  positional  rates 

Given  the  coordinates  x,  y,  h  of  the  aircraft's  centre  of  gravity,  the 
pilot's  position,  needed  to  drive  outside  world  displays,  is  obtained  from 


H* 

< 

N 

X  + 

(x  cos  0 

peg 

♦  z 

peg 

sin 

0) 

cos  t|> 

yTV 

y  + 

(x  cos  0 

peg 

♦  z 

peg 

sin 

e) 

sin  [p 

hTV  * 

h  + 

(x  sin  0 

peg 

-  z 

peg 

cos 

0) 

where 

X  “ 

peg 

X 

p 

+  Ax  c 
cg  ref 

i 

Z  ” 

peg 

z 

p 

-  Az  c  , 
cg  ref 

* 

and  Xp,  z  are  the  ooordinates  of  the  pilot's  eye  point,  assumed  to  be  on  the 
aircraft  centre  line,  and  defined  relative  to  the  reference  centre  of  gravity 
position. 

If  we  put 


R  »  x  cos  0  +  z  sin  0 

xp  peg  peg 

R.  -  x  sin  0  -  z  cos  0 

np  peg  peg 
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then 


16 


TV 


yTV 


TV 


x  -  R  iL  sin  ip  ~  R.  0  cos  ip 
xp  hp 


y  ♦  R  ip  cos  ip  -  R,  f-  sin  vl’  r 
xp  hp 


h  ♦  R  9 
xp 


(25) 


i .  l»  Accelerations  in  body  axes 

Linear  accelerations  at  the  aircraft's  centre  of  gravity,  at  the  pilot's 
head  and  elsewhere  U'J  for  accelerometers)  are  needed. 

Referred  to  body  axes,  the  components  (in  units  of  g)  of  the  acceleration 
of  the  aircraft's  centte  of  gravity  are 


xcg 


yep 


zcg 


F  /W 
tX 


'ty/u  > 


FT2/"  • 


(26) 


These  are  referred  to  as  the  'specific'  accelerations  and  are  the  accelera¬ 
tion  components  that  would  he  measured  by  a  set  of  orthogonal,  body-fixed, 
accelerometers  aligned  with  the  body  axes  and  located  at  the  centre  of  gravity  of 
the  aircraft.  'Absolute'  accelerations,  including  the  gravity  components,  are 
then 


a 

xacg 

5S 

a 

xcg 

sin  0 

a 

xcg 

S,3 

a 

yacg 

• 

a  + 

yeg 

sin  p  cos  0  = 

a 

yeg 

+ 

S23 

a 

zacg 

- 

a  ♦ 

zcg 

cos  4>  cos  0  = 

a 

zcg 

♦ 

S33  ’ 

I t  shou 1 d 

be  remembered  that 

a  i s  no  t 

xacg 

equal 

to 

axis  system  is 

rotat ing. 

as 

explained  in  any  text 

,  such 

as 

Ktkin 

At  an  arbitrary  location  L(x^,yj  ,z  )  ,  the  components  of  specific 
acceleration  are 


xcg 


a  + 
yeg 


a  + 
zcg 


2  2  .  • 

xL(q  ♦  r  )  -  yj  (pq  -  r)  -  z[ (pr  ♦  q) 

2  2  • 

x^ (pq  ♦  r)  -  yL(p~  +  r“)  +  Zj  (qr  -  p) 

Xt (pr  -  q)  +  yL(qr  +  p)  -  zL(p~  ♦  q‘) 


(27) 
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The  indication  of  an  accelerometer  or  a  'g'  meter  may  then  be  derived  from 
equations  (27)  given  the  location  of  the  device.  Similarly  the  linear  accelera¬ 
tion  components  at  the  pilot's  station  (x^,  z^)  may  be  derived  for  use  in 
driving  motion  systems.  Further  details  of  the  accelerations  computed  are  given 
in  the  description  of  the  SACCBOD  routine  in  Appendix  A. 


3.10  Wind,  wind  shear  and  turbulence 
Wind 

A  datum  mean  wind  is  defined  in  speed  and  direction  by  VWK1()  and  *l'w  » 
from  which  the  components  in  earth  axes  are  obtained  as 


V  =  — V  cos  lb 

WNLO  WKTO  ^W 

V  *  — V  sin  ib 

WELO  WKTO 


(28) 


bearing  in  mind  the  convention  that  when  i |/w  is  zero,  the  wind  is  from  the  north. 
The  vertical  component  of  mean  wind,  »  i-s  normally  assumed  to  be  zero. 


Wind  shear 


Wind  shear  is  obtained  as  a  multiplying  factor  f  according  to  altitude, 
and  the  wind  components  at  height  are  then 


WNL 

fV 

WNLO 

WEL 

fV 

WELO 

(29) 


This  shear  is  effective  only  in  magnitude,  not  in  direction.  Three  choices  are 
available:  no  shear  (f  »  1.0),  logarithmic  profile  and  linear  profile.  The 
initialisation  process  allows  for  actual  wind  at  the  initial  height  and  sets  up 
the  aircraft's  track  to  give  a  desired  initial  heading. 

Turbulence 

Turbulence  in  three  orthogonal  directions  can  be  added  to  the  components 

of  the  mean  wind  to  give  total  wind  components.  Gust  velocity  components 

.12  ... 

u,.,  v„,  w„  are  calculated  by  a  new  technique  ,  based  on  the  Statistical  Discrete 

°  °  °  13.  . 

Gust  theory  of  J.G.  Jones  ,  which  allows  the  intermittent  character  of  the  gener¬ 
ated  turbulence  to  be  controlled.  Scaled  components  of  turbulence  u^.,  v^,,  w^ 
are  obtained  by 
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u  .  u . 

Slg  (. 


V  .  V 
sig  0 


w  .  w„ 
SIR  G 


whore  u  .  ,  v  .  ,  w  .  are  desired  root-mean-square  intensities.  These  turbu- 
sig’  sir’  sir 

lence  components  are  considered  to  be  along  the  wind,  across  the  wind  and  verti¬ 
cal,  so  that  the  total  fluctuating  wind  components,  in  earth  axes,  are 

VWN  ”  VWNL  ~  (UT  COS  *W  '  VT  Sin  V 

VWE  “  vwkl  '  (VT  cos  *w  +  UT  sin  V  \  <3I) 

VWD  "  VWDL  "  WT  . 

3.11  Properties  of  the  atmosphere 

The  properties  of  the  ICAO  International  Standard  Atmosphere**  are  calcu¬ 
lated  up  to  a  maximum  height  of  65blb  ft  by  the  routine  SATMOS  (see  description 
of  SVEL0C2  in  Appendix  A).  Given  altitude,  it  returns  the  density  ratio,  speed 
of  sound  ratio,  temperature  ratio  and  pressure  ratio  for  a  standard  day.  A 
routine  ATMOS  is  also  available  for  hot  days  but  this  is  not  integrated  into 
SESAME.  Strictly  the  altitude  input  should  be  ’ geopotential 1  but  below  b5b)b  ft 
it  is  adequate  to  treat  geometric  and  geopotential  altitude  as  interchangeable. 

The  user  may  optionally  choose  to  use  constant  atmospheric  properties  (as 
in  many  classical  studies)  or  else  allow  standard  variation  with  altitude.  A 
software  flag  KISA  controls  this  option. 

4  TECHNICAL  IMPLEMENTATION 

Section  3  has  defined  a  set  of  equations  which  compose  the  standard  processes 
of  any  aircraft  simulation.  These  equations  are  distributed  among  a  series  of 
modules  for  solution  by  digital  computer.  This  section  briefly  describes  the 
routines  and  their  functions.  Full  details  of  each  routine  are  given  in  Appendix 
A  and  listings  are  given  in  Appendix  B.  In  general,  system  routines  have  names 
beginning  with  S  .  The  routines  fall  into  four  categories: 


Initialisation 
Rotational  motion 
Translational  motion  - 
Utility  functions 


SINIT,  SYSCOM 
SDCOS,  SEULER,  SACCROT 

SVELOCI,  SALFBET,  SACCLIN,  SPATH,  SVEL0C2 
STV,  SILS,  SWINn,  SACCBOD,  SATMOS,  SCOUNT  . 
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All  the  subroutines  described  in  this  Report  are  coded  in  Fortran.  The 
particular  dialect  used  is  Xerox  Extended  Fortran  IV,  designed  for  the  Xerox 
Sigma  range  of  computers.  It  was  not  intended  to  use  only  those  features  of 
the  language  embodied  in  Standard  Fortran  IV  and  a  few  non-standard  features  have 
in  fact  been  used.  These  include 

variable  names  of  up  to  8  characters, 

initialisation  of  COMMON  variables  within  subroutines, 
use  of  the  NAMELIST  feature  (in  SYSCOM  only). 

Conversion  of  the  subroutines  to  run  on  another  computer  should,  however,  pose 
few  problems, 

A. 1  Initialisation 

Two  routines  (SINIT,  SYSCOM)  are  specifically  concerned  with  performing 
certain  start-up,  or  initialisation,  functions.  Other  initialisation  functions 
may  be  performed  internally  in  the  other  routines. 

SINIT  calculates  various  aircraft  related  constants,  sets  up  initial  height, 
derives  initial  values  of  atmospheric  properties  and  initialises  speeds  and 
attitudes  allowing  for  wind. 

SYSCOM  sets  up  communication  with  system  variables  by  creating  a  NAMELIST 
table  and  reads  'semi-permanent'  data  changes  from  a  file. 

Execution  of  most  routines  during  an  initialisation  pass  is  organised 
intrinsically  by  the  SL1  model  program  (see  section  5.2).  Many  routines  may 
need  to  perform  preliminary  calculations.  If  their  nature  is  such  that  they 
must  be  performed  onae  only  (otherwise  errors  will  result)  then  a  local  flag  must 
be  created  and  the  routine  structured  as  in  Fig  3a,  with  the  flag  set  'off'  on 
entry  (eg  by  a  DATA  statement).  However,  if  there  are  no  reasons  why  initial 
calculations  should  not  be  performed  more  than  once  there  is  a  standard  flag 
available  in  the  system  COMMON,  JJCOMP,  which  takes  the  value  0  until  the 
'compute'  button  is  pressed,  and  thereafter  is  I  .  In  this  case  the  structure 
in  Fig  3b  can  be  used. 

A. 2  Rotational  motion 

A  diagram  of  the  information  flow  is  shown  in  Fig  A.  Rotational  motion  is 
calculated  first,  principally  to  provide  the  direction  cosines  required  by  the 
translational  motion. 
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SACCROT  calculates  the  three  components  of  angular  acceleration,  in  body 
axes,  using  equations  (19),  given  total  moments  (supplied  by  the  user)  and  inertia 
constants,  equations  (20),  (calculated  by  S1N1T). 

SKULER  calculates  the  rate  of  change  of  body  attitude  angles  from  the 
angular  velocity  components,  using  equations  (1). 

SIX'OS  calculates  the  nine  direction  cosines  S( (  etc  from  the  three 
attitude  angles  G,  $  using  equations  (3). 

4.3  Translational  motion 

A  diagram  of  the  information  flow  is  shown  in  Fig  3. 

SVELOCI  derives  components  of  velocity,  relative  to  the  air,  in  earth  axes 
using  equations  (9)  and  transforms  to  body  axes,  using  equations  (2). 

SALFBET  takes  the  body  axes  velocity  components  and  calculates  a,  6,  a,  S 
using  equations  (12)  and  (13). 

SACCLIN  transforms  total  force  components  in  body  frame  to  inertial  (earth) 
frame,  using  equations  (7),  then  calculates  the  translational  acceleration  using 
equations  (8).  It  may  sometimes  be  desirable  to  separate  the  vertical  from  the 
two  horizontal  components  so  that  vertical  motion  can  be  solved  more  frequently 
for  better  simulation  of  undercarriage  dynamics. 

SPATH  calculates  flight  path  angles  from  velocities  in  earth  axes,  using 
equations  (15)  and  (16). 

SVEL0C2  calculates  resultant  airspeed,  dynamic  pressure,  Mach  number  etc, 
using  SATMOS  for  atmospheric  properties.  Some  features  are  at  present  only  valid 
for  Mach  number  less  than  1.0,  as  indicated  in  the  routine. 

4.4  Utility  functions 

Fig  t>  shows  the  flow  of  information  among  the  utility  routines,  all  of  which 
have  something  to  do  with  the  simulation  environment,  visual  display,  motion 
cues,  wind  etc. 

STV  calculates,  from  equations  (22)  to  (25),  positions  and  velocities  to 
drive  the  TV  visual  system,  and  also  handles  the  logic  of  belt  positioning. 

SILS  calculates  one  or  two  segment  ILS  guidance  beams. 

SWIND  controts  the  generation  of  turbulence  and  adds  it  to  wind,  modified 
by  shear  effects  (if  any). 
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SACCBOD  calculates  accelerations  in  body  axes  of  the  aircraft  centre  of 
gravity,  from  equations  (26),  and  pilot  station,  accelerometers  etc,  from 
equations  (27). 

SATMOS  calculates  the  properties  (in  ratio  form)  of  the  ICAO  international 
standard  atmosphere.  It  is  an  existing  library  subroutine  and  is  not  otherwise 
described  in  this  Report. 

A  further  utility  routine,  SCOUNT,  not  included  in  Fig  b,  has  nothing  to  do 
with  the  simulation  of  aircraft,  but  assists  in  the  management  of  the  calculation 
process,  by  setting  an  'initialisation  complete'  flag,  by  keeping  track  of  the 
sub-steps  (if  any)  of  the  integration  routine,  and  by  picking  up  the  current  run 
number  from  the  data  logging  system. 

4,5  Units 

With  Aeronautics  still  retaining  the  use  of  such  units  as  knots,  it  is 
difficult  at  present  to  make  a  wholehearted  conversion  to  SI  units.  However, 

the  equations  set  up  are  self-consistent  so  that  a  force  in  newtons  acting  on 

■> 

a  mass  in  kilograms  will  produce  an  acceleration  in  m/s“  and  likewise,  a 
force  in  pounds  and  mass  in  slugs  will  yield  ft/s'.  The  routines  as  set  up  are 
intrinsically  in  the  feet,  pound,  second  system  but  alteration  of  a  few  conver¬ 
sion  constants  (gravity,  etc)  and  datum  values  (.air  density,  etc),  would  enable 
the  whole  system  to  work  in  SI  units. 

5  OVERALL  PROP  RAM  STRUCTURE 

5.  I  Introduction 

The  program  which  provides  the  general  framework  for  all  the  system  routine 
is  written  in  the  simulation  language  SL1  . 

The  basic  structure  (Fig  7)  consists  of  an  INITIAL  region,  for  start-up 
calculations,  a  DYNAMIC  region  and  a  TERMINAL  region.  (The  TERMINAL  region, 
however,  has  no  relevance  to  real-time  simulation.)  This  is  a  common  form  for 
continuous  system  simulation  languages,  like  SLI  or  CSMP.  The  DYNAMIC  region 
contains  one  or  more  DERIVATIVE  sections.  These  DERIVATIVE  sections,  each  of 
which  can  have  its  own  integration  algorithm  and  step  size,  contain  the  routines 
which  are  the  core  of  a  simulation.  During  real-time  operation,  the  code  gener¬ 
ated  from  each  DERIVATIVE  section  is  executed  repetitively  to  produce  the  desired 
solution  as  a  function  of  time.  A  listing  of  an  aircraft  simulation  program  with 
one  DERIVATIVE  section,  or  loop,  is  given  in  Appendix  C,  and  corresponding  flow 
charts  of  the  execution  sequence  are  shown  in  Fig  8.  The  reason  for  using  SLI 
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is  that  two  major  facilities  are  provided  by  the  language.  These  are  synchronisa¬ 
tion  with  real-time  and  a  centralised  integration  scheme  with  a  choice  of  five 
algorithms.  Further  details  are  contained  in  the  Reference  Manuals'"  and  in 
Ref  1. 

The  SLI  program  is  translated  into  a  series  of  Fortran  modules.  Communi¬ 
cation  is  organised  using  Fortran  labelled  COMMON.  To  permit  system  COMMON  to 
handle  the  main  variables,  however,  some  duplication  of  variables  is  necessary 
in  order  for  SESAME  routines  to  communicate  with  the  integration  operators. 

This  duplication  is  arranged  through  PROCEDURAL  blocks,  with  the  subterfuge 
that  only  those  variables  to  be  exchanged  with  the  integration  processes  are  con¬ 
tained  in  the  argument  list.  Other  variables  thus  remain  invisible  to  the  SLI 
translator  and  so  may  be  placed  in  system  COMMON.  If  they  were  visible  to  the 
SLI  translator  they  would  be  placed  in  a  separate  COMMON  area,  no  longer  under 
control  of  the  user. 

For  example,  a  PROCEDURAL  block 
PROCEDURAL  (=VV) 

V  =  VV 

END 

picks  up  W  from  the  SLI  labelled  COMMON  and  places  its  value  in  V  .  In  the 
opposite  direction 

PROCEDURAL  (WDOT-) 

WDOT  =  VDOT 

END 

transfers  VDOT  to  WDOT  ,  for  use  in  the  integration  statement 
VV  -  INTEG  (WDOT,  WIC) . 

5.2  Init ial isation 

In  the  normal  course  of  events,  the  code  generated  by  the  SLI  translator 
causes  the  DERIVATIVE  section(s)  to  be  executed  just  once  at  the  end  of  INITIAL 
to  calculate  the  initial  values  of  all  the  derivatives.  This  has  been  augmented 
by  explicit  code  (see  near  the  'END'  of  the  INITIAL  region  in  the  SLI  listing  in 
Appendix  C)  to  force  these  calculations  to  be  done  twice.  This  is  necessary 
because  routines  at  the  end  of  the  sequence  (such  as  SVEL0C2)  produce  variables, 
such  as  dynamic  pressure,  needed  early  on.  Once  the  whole  system  is  running, 
there  is  no  problem.  This  slight  idiosyncracy  needs  to  be  remembered  by  the 
user.  The  other  important  point  about  initialisation  is  that  at  the  start  of  a 
new  run,  a  fresh  copy  of  the  program  is  loaded  into  main  store  from  a  disc  file, 
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so  all  parameter  values  are  restored  to  their  datum  values.  However,  the 
'retained  changes'  facility  (see  section  7.7)  permits  revised  parameter  values 
to  be  retained  from  run  to  run. 

5.3  PV100  -  operator  dialogue 

A  set  of  programs  and  routines,  collectively  known  as  the  parameter  varia¬ 
tion  package,  exists  to  enable  the  user  to  interact  on-line  with  the  parameters 
of  the  simulation.  The  principal  item  is  PV100,  the  operator  dialogue  routine, 
which  provides  a  wider  range  of  facilities  than  the  'interpreter'  inherent  in 
SLI.  PVI00  enables  the  user,  on-line  and  during  any  phase  of  a  simulation  flight 
to  access,  display  and  amend  any  variable  within  the  scope  of  the  model  program 
and  its  subroutines.  Variables  may  be  referred  to  by  name  (the  principal 
function  of  routines  SYSCOM  and  USERCOM  is  to  provide  tables  of  names  for  use  by 
PV100),  and  in  addition  the  contents  of  any  memory  location  may  also  be  displayed 
A  display  may  be  in  alternative  styles,  eg  binary,  hexadecimal,  integer  or  'real' 
format.  Selected  changes  may  be  'retained'  to  operate  for  subsequent  runs, 
otherwise  the  change  vanishes  at  the  end  of  the  current  run.  A  readable  memory 
dump  may  be  printed,  such  as  all  the  parameters  in  system  COMMON.  There  is  also 
a  facility  for  forcing  an  'automatic  hold'  or  freeze  if  a  variable  falls  within 
defined  limits.  A  description  of  the  full  facilities  is  given  elsewhere1  . 

During  initialisation  there  are  four  explicit  calls  to  PVI00  inserted  at 
strategic  points  in  the  INITIAL  region,  as  can  be  seen  by  inspection  of  the  SLI 
listing  in  Appendix  C.  These  calls  to  PV100  are  only  executed  if  a  switch  (Desk 
Switch  1)  on  the  simulator  control  desk  is  set  on.  These  individual  calls  to 
PVIOO  may  be  identified  by  reference  to  the  indicator  ITP,  as  follows. 

ITP  *  1  After  initial  conditions  or  retained  changes  have  been 

input,  but  before  SINIT.  Enables  IC  values  to  be  altered, 
and  other  changes  to  be  inserted  before  the  main  initialisa¬ 
tion  calculations  of  SINIT.  In  particular,  if  the  frame  time 
is  to  be  altered,  it  should  be  done  here,  using  FRAMET1,  DELT1 
for  Loop  1  and  similar  variables  for  Loop  2  (if  present). 

ITP  =  2  After  SINIT  but  before  the  integrator  IC  values  are  actually 

set  up. 

ITP  =3  At  end  of  INITIAL  but  before  analogue  to  digital  converter 
(ADC)  read  and  DERIVATIVE  section(s)  initialised. 

ITP  -  4  After  DERIVATIVE  section(s)  initialised  and  ADC  read.  DAC 

values  are  calculated  but  not  yet  set. 


COMMUNICATION 


b 

6.  1  Introduction 

Communication  among  all  the  routines  which  comprise  the  aircraft  model 
program  is  achieved  by  Fortran  labelled  COMMON.  Labelled  COMMON  designates  a 
block  of  contiguous  memory  locations  which  may  be  accessed  by  any  subroutine  in 
which  the  COMMON  block  is  declared.  This  technique  provides  flexibility  and  also 
avoids  the  time  penalty  associated  with  using  arguments  in  subroutine  calls. 

Two  large  blocks  of  labelled  COMMON  are  employed.  One,  named  SYSTEM,  is  an 
inherent  part  of  the  standard  subroutines  and  contains  all  the  state  variables 
and  other  parameters  common  to  any  aircraft  simulation.  This  is  of  fixed  size 
and  has  its  variables  in  predefined  locations.  The  other,  named  USER,  is  set  up, 
in  terms  of  size  and  content,  by  the  user  to  handle  variables  specific  to  his 
particular  simulation,  eg  the  aerodynamic  force  coefficients,  engine  forces  etc. 

b . 2  System  COMMON 

This  is  defined  as 

COMMON  /  SYSTEM  /  A(1000),  L(400) 

where  the  array  A  contains  real  variables  and  L  contains  integer  variables. 

An  index  to  all  the  variables,  in  numerical  and  alphabetical  order,  is  provided 
in  Appendices  D  and  E. 

If  the  user  wishes  to  pick  up  the  current  value  of  angle  of  attack  (say) 
tie  needs  to  include  in  his  routine 

COMMON  /  SYSTEM  /  A(IOOO),  L(AOO) 

EQUIVALENCE  (ALFAD  ,  A ( I  1 2) ) 

and  then  ALFAD  can  be  used  freely  in  his  routine.  System  variables  must  not  be 
altered  by  the  user. 

b . 3  User  COMMON 

This  is  defined  for  example  as 

COMMON  /  USER  /  B(1000),  M(100) 

where  the  array  B  is  for  real  variables  and  M  for  integer  variables.  Names 
and  sizes  of  these  arrays  are  defined  by  the  user,  but  the  COMMON  name  USER  must 
be  employed.  Sizes  should  be  kept  as  small  as  possible,  in  order  to  avoid 
sterilising  areas  of  core  unnecessarily.  However,  if  the  size  is  too  small,  the 
user  will  have  to  increase  it  frequently,  with  consequent  work  updating  all 
routines  in  which  the  COMMON  statement  appears. 
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The  user  will  choose  his  own  variable  names  and  allocate  space  in  the 
COMMON  arrays.  This  will  normally  be  done  by  EQUIVALENCE  techniques,  for  flexi¬ 
bility.  (Routines  could,  in  fact,  be  created  and  checked  in  the  first  place 
without  considering  inter-routine  communication.) 

6.4  Documentation  aids 

A  utility  program  has  been  created,  called  COMMLIST,  to  assist  the  user  in 

keeping  track  of  his  variable  names  and  locations.  As  he  creates  a  variable,  the 

user  punches  a  card  defining  the  variable  name,  its  meaning  and  the  routine  in 

1 8 

which  it  is  calculated.  This  is  similar  to  a  scheme  described  by  Bean  .  The 
program  COMMLIST  can  then  produce  an  index  of  names  in  alphabetical  order,  of 
locations  in  numerical  order  or  of  names  according  to  the  routine  in  which  they 
are  created  or  used.  A  brief  guide  to  the  program  COMMLIST  is  given  in  Appendix  F. 
The  lists  in  Appendices  D  and  E,  G  and  H  have  been  produced  by  this  program. 

6 . 5  Access  to  variables  at  run-time 

In  order  to  be  able  to  access  his  variables  at  run-time,  to  inspect  and 
change  values  via  PVIOC  ^  ,  the  user  must  create  a  NAMELIST  table  of  his  variable 
names.  This  is  achieved  in  the  routine  USF.RCOM,  which  must  be  written  by  the 
user.  The  bulk  of  this  routine  will  be  COMMON  (or  EQUIVALENCE)  statements,  refer¬ 
encing  the  variables  by  name.  Its  structure  will  be  identical  to  the  routine 
SYSCOM  which  does  the  same  job  for  system  variables.  A  listing  of  SYSCOM  is 
included  in  Appendix  B. 

6.6  Hybrid  input/output  (ADC,  DAC) 

6.6.1  Introduction 

Hybrid  input /output ,  that  is  analogue  to  digital  conversion  (ADC)  and 
digital  to  analogue  conversion  (DAC),  is  integrated  into  SESAME  so  that  the  user 
merely  has  to  define  entries  in  a  set  of  tables  incorporated  in  system  COMMON. 
Conversion  and  scaling  then  occurs  automatically.  The  general  form  of  conver¬ 
sion  required  is  of  the  form 

y  =  ax  +  b 

where  a  is  a  scale  factor  and  b  is  a  bias.  For  ADC,  x  would  be  a  scaled 
value  and  y  a  value  in  'engineering'  units  whereas,  for  DAC,  x  would  be  in 
engineering  units  and  y  a  value  scaled  appropriately  for  output.  The  present 
implementation  does  not  include  any  bias. 
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Data  transfer  in  both  directions  employs  a  special  piece  of  hardware  known 
as  the  DMS12  (Direct  Memory  Sub  System)  which  is  a  processor  capable  of  operating 
in  parallel  with  the  computer's  central  processor.  Use  of  this  device  provides 
for  an  element  of  parallel  operation  which  is  exploited  in  the  software  to  save 
execution  time. 


Overall,  the  software  scheme  offers  the  user  convenience  -  because  he  only 
has  to  set  up  data  tables;  flexibility  -  because  the  tables  can  be  changed  on 
line;  and  speed  -  because  the  hardware  is  fully  exploited.  The  operating  routines, 
coded  in  assembler,  are  not  described  in  this  Report. 


6.6.2  Analogue  to  digital  conversion 


Analogue  to  digital  conversion,  and  scaling,  are  executed  prior  to  entry 
into  the  Derivative  Section  of  the  SL1  model  program,  where  the  dynamic  equations 
are  solved.  The  main  scaling  calculation  performed  is  of  the  form 
YADC(I)  =  ADC (NADC (I) )  *  AADC(I) 

where  ADC  is  the  array  of  raw  converted  values,  nominally  in  the  range  ±1.0, 
held  in  the  SL1  labelled  COMMON  Z99999. 


NADC  is  the  array  of  ADC  channel  numbers 


AADC  is  the  array  of  scaling  factors,  set  up  by  the  user 

YADC  is  the  array  of  input  values  in  engineering  units  (degrees  etc) 
ready  for  use. 


All  32  channels  are  read  in  raw  form  into  the  array  ADC,  but  scaling  and 
transfer  to  YADC  stops  when  the  first  zero  element  in  NADC  is  reached.  Because 
of  the  parallel  nature  of  the  hardware  (DMS12,  mentioned  above)  scaling  is  begun 
as  soon  as  the  first  raw  converted  value  is  available.  Conversion  of  further 
channels  then  proceeds  concurrently  with  scaling.  To  avoid  the  danger  of  scaling 
a  channel  before  it  has  been  converted,  it  is  important,  therefore,  that  ADC 
channels  are  kept  in  numeric  order  and  are  used,  as  far  as  possible,  without  large 
gaps.  For  example,  to  use  channels  1,  3,  7,  9  is  better  than  1,  17,  24,  30.  To 
use  1,  2,  3,  4  is  best  of  all.  Thus,  for  simplicity,  NADC(I)  =  I. 

The  data  structure  being  in  table  form  allows  considerable  operational 
flexibility.  For  example,  in  the  event  of  failure  of  ADC  channel  2,  the  variable 
being  input  on  this  channel  can  be  redirected  to  channel  9  (say)  simply  by 
repatching  on  the  analogue  computer  and  by  changing  the  contents  of  NADC (2)  from 
2  to  9  via  the  keyboard.  The  scaling  AADC(2)  would  not  need  to  be  altered.  The 
necessary  changes  can  all  be  achieved  on-line  and  no  program  changes  are  necessary. 
However,  the  warning  of  the  previous  paragraph  should  not  be  overlooked. 


At  present  32  ADC  channels  are  available  and  supported  by  software. 
Provision  has  been  made  for  enlargement,  however,  by  reserving  64  contiguous 
locations  for  each  of  the  arrays  AADC,  YADC,  NADC  in  the  COMMON  storage  area.1 

The  SL I  model  program  includes  in  its  source  code  (Appendix  C)  a  sequence 
of  the  form 

PROCEDURAL  (=ADC(1)) 

END 

in  order  to  ensure  that  the  key  driving  routine  Z99^1  is  called.  No  special  code 
has  to  be  created  by  the  user. 

The  user  does,  of  course,  wish  to  use  the  result  of  the  conversion  and 
scaling  processes,  as  held  in  YADC.  To  do  so  he  may  either  use  an  element  of 
YADC  directly,  as  in  the  following  example 
CL  =  CLETA*YADC (7)  +  . 

or  may  give  the  element  of  YADC  his  own  names  via  the  usual  EQUIVALENCE  declara¬ 
tions  viz 

DIMENSION  YADC (32) 

EQUIVALENCE  (YADC(l),  A (40 1 ) ) 

EQUIVALENCE  (ETAD  ,  YADC (7)) 

and  then 

CL  =  CLETA  *  ETAD  +  . 

To  make  his  variable  names  accessible  via  PVI00,  the  user  could  include  them  in 
subroutine  USERCOM. 

6.6.3  Digital  to  analogue  conversion 

Digital  to  analogue  conversion  occurs  once  per  frame,  after  the  main  inte¬ 
gration  loop.  The  DAC  software  scales  output  variables  to  be  in  the  settable 
range  for  DACs  and  DCUs  (otherwise  known  as  vari-DACs). 

In  principle,  the  same  calculation  technique  is  adopted  as  for  ADC.  The 
special  difficulty  for  variables  to  be  output  to  DACs  is  that  they  must  be 
selected  from  any  location  in  the  system  or  user  COMMON  arrays  (A  or  B)  and  must 
be  routed  either  to  proper,  four-quadrant  DACs  (the  first  32)  or  to  two-quadrant 
vari-DACs  (the  next  16).  An  additional  requirement  is  that  it  should  be  possible 
to  set  DACs  from  either  of  two  loops.  (This  description  refers  for  completeness 
to  two  loops.  For  the  moment,  however,  only  Loop  I  is  operational.) 

To  achieve  these  objectives,  data  for  four  arrays  must  be  defined  by  the 
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ADAC  array  of  scaling  factors.  A  variable  in  engineering  units  is 
multiplied  by  the  scaling  factor  in  the  appropriate  element  to 
give  a  result  in  the  desired  range  (nominally  ±1.0). 

NADC  is  the  array  of  pointers  defining  which  variables  are  to  be  con¬ 
verted.  The  pointers  may  refer  to  the  'A'  or  'B'  arrays  in  any 
sequence.  Distinction  between  the  source  arrays  is  controlled  by 
arrays  LIDAC,  L2DAC  below. 

L1DAC  defines  the  source  of  data  for  each  channel  to  be  converted  in 
Loop  1 . 

L2DAC  defines  the  source  of  data  for  each  channel  to  be  converted  in 

Loop  2.  Each  element  of  LIDAC,  L2DAC  can  take  the  value  0,  1  or  2: 

1  means  that  the  data  for  scaling  is  to  be  taken  from  array  A  and 

2  that  it  comes  from  array  B.  If  0  is  set,  no  scaling  is  performed. 

As  an  example,  if  NDAC(6)  =  112,  L1DAC(6)  =  1,  ADAC(6)  =  0.1  then  DAC 
channel  6  is  to  be  used  to  output,  at  the  end  of  Loop  1,  variable  A ( 1 12),  ie  ALFAD, 

angle  of  attack  in  degrees,  scaled  to  ±10.0. 

This  structure  provides  flexibility  by  enabling  changes  to  be  introduced 
on-line.  For  example,  by  changing  NDAC(6)  from  112  to  241  (say),  and  changing 
the  scaling  ADAC(6)  to  suit,  a  different  variable  from  the  same  array  can  be  out¬ 
put  on  a  given  DAC.  The  source  array  (A  or  B)  can  also  be  easily  changed,  via 
LIDAC. 


Operational  considerations 


(I)  Of  the  48  DAC  channels  available,  the  first  32  are  true  DACs,  while 
the  rest  are  only  two-quadrant  devices.  The  user  needs  to  remember  this 
when  allocating  channels,  so  that  33-48  are  only  used  for  variables  which 
do  not  change  sign,  eg  airspeed.  Provision  has  been  made  for  expansion  to 
64  channels  by  reserving  extra  space  for  ADAC,  NDAC  in  the  COMMON  area. 


(2)  For  those  DAC  channels  not  used,  the  corresponding  elements  of  ADAC 
should  be  set  to  zero. 


(3)  A  set  of  default  hardware  addresses  associated  with  channel  numbers 
is  defined  in  the  table  IDACAD.  DAC  addresses  can  be  reassigned  during 
initialisation. 


(4)  Timing  tests  have  given  the  following  results: 


Convert  and  scale  32  ADC  channels  1.2  ms 

Scale  and  convert  48  DAC  channels  3.8  ms  . 


All  channels  are  converted,  but  only  those  specified  are  scaled  and  stored 
so  that  time  can  be  saved  if  fewer  channels  are  needed. 
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6.6.4  User  summary  of  system  COMMON  relevant  to  ADC/DAC 


■ 

301  364 

401  464 

501  564 

-  1 

1 

1  AADC  64 

!  YADC  64 

1  ADAC  64 

1 

A  System  COMMON  array  (real) 

AADC*  Array  of  scaling  factors  for  ADC 
YADC  Array  of  converted  variables 

ADAC*  Array  of  scaling  factors  for  DAC 


| 

11  74 

I 

82  144 

■ 

305  352 

1 

1  NDAC  64 

1 

1  NADC  64 

■ 

1  LI DAC  48 

_ 

1  L2DAC  48 

L  System  COMMON  array  (integer) 

NDAC*  Array  of  pointers  defining  variables  in  'A'  array  and 
user's  'B'  array  to  be  converted 

NADC*  Array  of  ADC  channel  numbers 

L1DAC*  Array  of  DAC  source  pointers  for  Loop  1 

L2DAC*  Array  of  DAC  source  pointers  for  Loop  2 

6.7  Hybrid  input/output  (discretes) 

6.7.1  Introduction 

Discretes  are  single-bit  logic  lines  enabling  on/off  information  to  be 
communicated  into  and  out  of  the  computer.  These  lines  are  read  and  set  auto¬ 
matically  and  so  are  conveniently  available  to  the  user. 

6.7.2  Summary  of  lines  available 

A  number  of  lines  are  available  in  each  direction 

(a)  between  the  digital  computer  (Sigma  8)  and  the  analogue  computer  (AD4) 
and 

(b)  between  the  digital  computer  and  other  external  hardware,  such  as  the 
simulator  cockpits  and  the  control  desk.  The  numbers  and  types  are  sum¬ 
marised  in  the  table  below,  together  with  the  relevant  driving  routine . 

*  These  items  must  have  values  supplied  by  the  user.  It  is  recommended  that,  for 
convenience,  all  such  data  be  located  in  one  routine,  such  as  OUTSR,  or  a  self- 
contained  'information'  routine. 
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All  these  names  are  defined  in  SYSCOM,  to  be  accessible  by  PV100.  For  the 
user  to  obtain  access  to  ISLR(3),  for  example,  he  will  include  in  his  routine 

COMMON/SYSTEM/A (1000) ,  L(400) 
and  either  EQUIVALENCE  (ISLR(l),  L ( 1 5 1 ) ) 

DIMENSION  ISLR(16) 

or  EQUIVALENCE  (ISLR(3),  L(153)) 

The  former  is  more  general.  Discrete  information  is  then  available  for  use  in  such 
constructs  as 

IF  (ISLR(3).EQ.l)  CALL  XYZ 
6.7.4  Dedicated  functions 

Some  of  these  discrete  lines  are  dedicated  to  certain  functions  within 
SESAME,  as  summarised  in  the  following  list. 


Line 

Purpose 

Routine 

I SLR (4) 

Slew  forward 

STV 

ISLR(5) 

Slew  back 

STV 

ICLR( 1 6) 

TV  Desync 

STV 

IDS  ( 1 ) 

PV100  select 

SL1 

IDS(5) 

TV  cycle 

STV 

IDS(6) 

X  TV  reset 

STV 

IDS (7) 

Y  TV  reset 

STV 

IDS (8) 

Turbulence 

SWIND 

7  USER*  S  GUIDE 

7 . 1  What  the  system  expects  from  the  user 

The  user  is  required  to  provide  the  three  components  of  total  force  (FTX, 
FTY,  FTZ)  and  the  three  components  of  total  moment  (XLLTOT,  XMMTOT,  XNNTOT) ,  all 
in  body  axes.  Thereafter  the  appropriate  integrations,  resolutions  etc  occur, 
with  the  system  providing  back  to  the  user  all  the  state  variables  from  which  the 
forces  etc  are  generated.  Various  constants  must  also  be  provided.  These  are 
described  in  section  7.3. 

7.2  Typical  user  routines 

The  user  is  required  to  create  the  following  set  of  routines,  which  may 
perform  some  or  all  of  the  tasks  indicated. 
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'Input' 

CONTROLS  takes  the  raw  control  inputs  from  the  pilot,  applies  any  non¬ 
linear  gearings,  computes  autostabiliser  contributions  and 
finishes  up  with  the  total  control  surface  deflections,  for 
use  by  the  aerodynamics  routines.  It  may  also  organise  the 
computation  of  engine  performance,  including  dynamics,  and 
total  thrust,  momentum  drag  etc. 

'Calculation  loop' 

TOTF  computes  aerodynamic  force  coefficients  etc,  finishing  up 

with  three  components  of  total  force  from  all  sources. 

TOTM  computes  aerodynamic  moment  coefficients  etc,  finishing  up 

with  three  components  of  total  moment  from  all  sources. 

'Output' 

OUTSR  computes  miscellaneous  functions,  eg  scalings  of  non-linear 

instruments,  navigation  and  guidance  and  handles  data  logging. 

'  Initial ' 

USERCOM  contains  the  names  of  all  the  variables  the  user  may  wish  to 
access  via  PVI00,  and  also  reads  semi-permanent  data  changes 
from  a  file. 

*  These  routines  are  actually  called  at  SL1  level.  All  are  executed  in 
a  DERIVATIVE  section,  except  USERCOM,  which  is  executed  in  INITIAL. 

*  All  these  routines  are  likely  to  be  umbrella  routines,  in  the  sense 
that  they  are  likely  merely  to  organise  the  calling  of  other,  more 
detailed,  routines.  For  example,  TOTF  may  call  a  routine  to  compute 
all  ground  reactions,  tyre  forces  etc. 

*  Ideally,  TOTF  and  TOTM  should  be  kept  separate,  so  that  they  can,  if 
necessary,  be  executed  in  different  loops  at  different  frame  fimes. 

*  Routines  in  the  calculation  loop  may  be  executed  up  to  four  times  pet- 
frame,  depending  on  the  integration  technique  employed,  eg  Trapezoidal, 
Runge-kutta  fourth  order.  All  other  routines  are  only  executed  once 
per  frame  (see  section  5.4). 

*  Each  routine  should,  as  far  as  possible,  contain  its  own  constants. 

*  The  method  for  communication  of  variables  other  than  by  lengthy 
argument  lists  has  been  described  in  section  6. 


Data  the  user  must  provide 


In  addition  to  the  primary  forces  and  moments  (FTX,  XLLTOT,  etc)  and  initial 
values  (section  7.5),  the  user  must  provide  values  for  a  number  of  constants  and 
parameters.  These  fall  into  such  categories  as  basic  data,  hybrid  input/output 
and  external  environment.  They  are  described  briefly  below  and  listed  in  full 
in  Appendices  G  and  H.  Typically  the  values  will  be  provided  by  DATA  statements 
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in  appropriate  subroutines  ( eg  TOTM,  TOTF),  unless  default  values  are  already 
provided,  in  which  case  new  values  should  be  provided  in  the  SYSCHNO  file 
(section  7.7).  Zero  values  should  be  set  up  if  the  parameter  is  not  relevant 


Aircraft  basic  data 

Moments  of  inertia 
Centre  of  gravity 
Span,  wing  area  etc 
Pilot  location 
Coordinates  of  slip  ball 
Coordinates  of  g  meter 
Coordinates  of  accelerometer  (AX3) 
Coordinates  of  accelerometer  (AY4) 
Coordinates  of  accelerometer  (AZ5) 

Timing 

Frame  time 
Step  size 

Hybrid  input/output 

Scaling  factors  for  analogue  to 
digital  conversion  etc,  (for  full 
details  refer  to  section  6.6) 

Output  discretes  (see  section  6.7) 

External  environment 
Turbulence  rms 
Turbulence  character 
Ship  speed 


XIX,  XIY,  XIZ,  XIZX 

XCGREF,  ZCCREF,  ZCG 

SPAN,  SWREF ,  STAIL,  XLTAIL,  CREF 

XP,  ZP 

XI,  Yl,  Zl 

X2 ,  Y2 ,  Z2 

X3,  Y3,  Z3 

X4 ,  Y4 ,  Z4 

X5 ,  Y5 ,  Z5 

FRAMET1 ,  FRAMET2 
DELT1,  DELT2 

AADC,  ADAC 

LIDAC,  L2DAC,  NDAC,  NADC 
ICLR,  ICO 

USIG,  VSIG,  WSIG 
NG,  SFRACG,  SRDECAY 
VSHIPKT 


Additional  integer  parameters  which  can  take  one  of  several  specific  values  and 
thereby  select  alternative  functions  are  described  in  the  next  section. 


7.4  Options  available 

A  number  of  choices  may  be  made  by  software  flags  or  hardware  switches. 
These  are  summarised  below,  with  default  values  where  relevant. 

7.4.1  Software  flags 

(a)  Variation  of  atmospheric  properties 
KISA  default  0  set  in  SINIT 

0  constant  sea  level  conditions  at  all  times 

I  standard  atmosphere  as  function  of  current  height 

-I  constant  properties  appropriate  to  initial  height  (HIC). 
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FOR  KISA  *  1,  SATMOS  is  called  within  SVELOC2  to  calculate  the  new  proper¬ 
ties.  Otherwise,  constant  values  are  set  up  within  SINIT. 

(b)  Selection  of  constants  for  TV  belt  in  use 

NTVB  no  default 

1  700: 1  model 

2  2000: 1  model 

3  5000:1  model 

(c)  ILS  on/off 

LSILS  no  default 

0  ILS  off 

I  ILS  on 

(d)  Type  of  ILS  beam 

ILSFLG  no  default,  value  obtained  from  XILSFLG  (see  section  7.5.2) 

1  3°  straight  beam 

*  o  , 

2  6  straight  beam 

3  6°  changing  to  3°  at  height  HKINK  (qv) 

(e)  Discretes  in  and  out  (see  also  section  6.7) 

IAD4SL  AD4  sense  lines  (default  1  set  in  DSCRT1) 

IAD4CL  AD4  control  lines  (default  1  set  in  DSCRT1) 

IDSCFL  discretes  in  to  Sigma  (default  1  set  in  DSCRT1) 

IPCOFL  discretes  out  from  Sigma  (default  0  set  in  DSCRT1) 

0  do  not  execute  at  all 

1  execute  in  DSCRT1,  ie  in  fast  loop 

2  execute  in  DSCRT2,  ie  in  slow  loop  (if  any) 

(f)  Wind  shear 

ISHR  default  1,  set  in  SINIT  ' 

1  no  shear 

2  logarithmic  profile 

3  linear  profile 

The  value  of  ISHR  is  obtained  from  XISHR  set  up  in  the  Initial  Conditions 

file, 

(g)  Control  of  random  number  generation  for  turbulence 
LSEED  default  0,  set  in  SWIND 

0  seeds  for  random  number  generation  derived  from  time  of  day 

)  constant  seeds  employed,  enabling  a  repeatable  turbulence 

sequence  to  be  obtained 
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External  switches 

Desk 

switch 

Name 

Purpose  when  set  on 

1 

Calls  PV100  during  INITIAL 

5 

LCYCLE 

TV  belt  cycles  continually 

6 

LXTVIC 

Resets  TV  x  position 

7 

LYTVIC 

Resets  TV  y  position 

8 

LTURB 

Turbulence 

TV 

slew 

LFWD 

TV  belt  (and  x  position)  moves 
at  maximum  rate 

forward 

TV 

slew 

LBACK 

TV  belt  (and  x  position)  moves 

backward 

at  maximum  rate 

7.5  Initial  conditions  and  initial  values 

7.5.1  Initial  conditions 

The  standard  SL1  model  embodies  12  integrations,  each  of  which  requires  an 
initial  condition.  These  are  given  in  the  table  below  as  * SL 1  names'.  Many  of 
these  variables  do  not  form  a  natural  or  convenient  set  for  the  user,  who  is  more 
interested,  for  example,  in  specifying  an  initial  airspeed  rather  than  components 
of  speed  relative  to  the  ground,  or  an  initial  angle  of  attack  rather  than  pitch 
attitude.  Hence  SINIT  performs  various  calculations  to  transform  a  user's  initial 
values  into  those  that  the  system  requires. 

One  particular  feature  is  that,  regardless  of  the  wind,  the  ground  speed 
components  are  derived  implicitly  to  satisfy  the  specified  initial  heading  PSIDIC 
so  that  any  datum  wind  or  wind  shear  conditions  are  not  revealed  inadvertently 
to  the  pilot.  More  detail  is  given  in  the  description  of  the  SINIT  routine  in 


Appendix  A. 

SL1  names 

Source 

Definition 

WKNIC 

Calculation 

Velocity  comp  rel  to  ground. 

north 

WKEIC 

Calculation 

Velocity  comp  rel  to  ground. 

east 

WKDIC 

Calculation 

Velocity  comp  rel  to  ground. 

down 

SXIC 

XIC 

X  position 

SYIC 

YIC 

Y  position 

SHIC 

HIC 

Height 

PPIC 

PDIC, 

default 

zero 

Rate  of  roll,  body  axes 

QQIC 

QDIC, 

default 

zero 

Rate  of  pitch,  body  axes 

RRIC 

ROIC, 

default 

zero 

Rate  of  yaw,  body  axes 

PPHIC 

PHIDIC,  default  zero 

Roll  attitude  angle 

TTHIC 

Calculation 

Pitch  attitude  angle 

PPSIC 

PSIDIC 

Heading  angle 
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7.5.2  Initial  values 

Initial  values  are  parameters  relevant  to  the  current  'run1  and  are  not  in 
general  true  initial  conditions  for  integration.  However,  those  initial  condi¬ 
tions  listed  above  that  are  not  calculated  are  obtained  from  a  set  of  initial 
values  created  by  the  user  in  a  file.  This  file  is  read,  during  the  INITIAL 


phase,  by 

the  routine  RDICFILE 

and  its  contents  stored 

in  a  20-element  array 

(A( 1 91 )  - 

A ( 2 1 0) )  which  is  part 

of  system  COMMON. 

At 

present  only  14  elements  are  actually  used  (as 

defined  by  the  value  of 

NVALS) . 

These  are 

Name 

Element 

Definition 

Units 

BETADIC 

A(  1 9 1 ) 

Sideslip  angle 

degrees 

ALFADIC 

A( 1 92) 

Angle  of  attack 

degrees 

CAMDIC 

A( 1 93) 

Climb  angle 

degrees 

VKTIC 

A( l 94) 

Airspeed 

knots 

XIC 

A( 1 95) 

X  position 

feet 

YIC 

A( 1 96) 

Y  position 

feet 

HIC 

A( 1 97) 

Height 

feet 

W 

A( 1 98) 

Aircraft  weight 

pound 

XCG 

A( 1 99) 

Aircraft  centre  of 
gravity 

fraction  of  ref  chord 

XILSFLG 

A(200) 

ILS  flag  (actually 
used  in  integer  form, 
ILSFLG) 

XISHR 

A (20 1 ) 

Wind  shear  flag 
(actually  used  in 
integer  form,  ISHR) 

VWKT0 

A(202) 

Datum  wind  speed 

knots 

PSIWD 

A(203) 

Wind  direction 

degrees 

PSIDIC 

A(204) 

Aircraft  heading 

degrees 

Spare 

A(205)  -  A(2 10) 

Additional  items  could  be  read  by  increasing  NVALS,  e.7  by  inputting  a  new 
value  from  the  semi-permanent  changes  file  (see  section  7.7)  read  in  SYSCOM. 

7.6  Modification  status 


To  keep  track  of  the  modification  status  of  each  routine,  it  is  recommended 
that  the  user  allocate  an  identifier  per  routine  and  assign  a  value  to  identify 
the  version.  For  example,  with  the  routine  TOTF 
DATA  MTOTF  /  1  / 
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to  indicate  version  1,  or 

DATA  MTOTF  /  171176/ 

to  indicate  the  date  of  the  current  version.  These  values  should  be  changed 
whenever  the  routine  is  modified.  By  placing  these  status  identifiers  in  user 
COMMON,  the  current  versions  of  each  routine  can  be  readily  checked  on-line. 

This  technique  has  been  adopted  for  all  the  SESAME  routines.  The  name  of 
the  version  identifier  is  the  routine  name  preceded  by  K  ,  eg  for  routine  SALFBET, 
the  identifier  is  KSALFBET  and  contains  the  date  in  integer  form.  These  identi¬ 
fiers  are  held  in  system  COMMON  (see  Appendix  E). 

Identification  of  a  new  version  of  the  executable  program,  or  load  module, 
is  desirable.  The  SL1  model  program  has  a  version  number  (IVERSION)  for  its 
source  code,  but  this  is  not  changed  if  a  new  load  module  is  built  without  chang¬ 
ing  the  SL1  source.  Identifying  individual  routines  as  above  will  help  but  a 
global  load  module  identifier  could  also  be  created  in  a  'history'  routine,  the 
only  purpose  of  which  is  to  set  a  version  number  for  the  load  module,  eg 
DATA  LHARRIER  /  240776  / 

and  perhaps  also  to  include  a  modification  history  in  'comment'  form. 

A  recent  modification  to  the  loader  utility  in  the  RBM  operating  system 
now  places  the  date  and  time  of  load  module  creation  into  the  program  header . 

This  information  can  be  inspected  from  the  keyboard. 

7.7  Changes  to  parameters 

It  is  possible  to  introduce  lasting  parameter  changes  into  a  completed  and 
working  aircraft  model  program  without  recoding. 

On-line  changes  can  be  achieved  from  a  keyboard,  during  any  mode  of  a  simu¬ 
lation,  through  the  program  PVIOO*^.  The  user  has  the  option,  at  the  time  of 
making  the  changes,  of  declaring  that  the  changes  be  'retained',  in  which  case 
they  are  automatically  copied  into  a  file  and  read  back  again  (by  PV300)  at  the 
start  of  each  subsequent  run.  These  changes  are  initially  name-orientated,  but 
PV300  relies  on  absolute  address  so  that  if  alternative  aircraft  model  programs 
are  used,  the  retained  changes  file  needs  to  be  reset. 

Semi-permanent  changes  to  a  user's  parameters  held  in  USER  COMMON  may  be 
achieved  by  creating  (from  cards)  a  file  called  USERCHNG  which  is  read  each  time 
the  routine  USERCOM  is  executed.  Changes  to  parameters  in  SYSTEM  COMMON  may  be 
achieved  by  a  similar  process,  via  the  file  SYSCHNG.  All  these  changes  work  only 
by  parameter  name,  as  explained  in  the  description  of  SYSCOM  in  Appendix  A. 
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7.8  Operational  considerations 

The  binary  versions  of  all  SESAME  routines,  in  relocatable  form  (ROM),  are 
held  in  files  on  disc  and  may  be  incorporated  when  building,  via  OLOAD,  the  final 
executable  program.  To  ensure  that  the  size  of  user  defined  COMMON  is  correctly 
set  up,  it  is  essential  that  a  user  ROM  containing  the  definition  of  USER  COMMON 
be  loaded  before  SDAC. 

If  data  logging  is  in  use,  the  correct  run  number  may  be  obtained  by 
inspecting  the  variable  NRUN. 

8  CONCLUSIONS 

A  system  of  equations  and  associated  computer  subroutines  has  been  developed 
to  make  it  easy  to  create  an  aircraft  model  program  for  real-time  flight  simula¬ 
tion.  Several  simulations  have  been  completed  and  it  has  been  found  to  be  very 
useful.  The  objectives  of  the  scheme  have  been  achieved.  In  particular,  the 
resulting  model  program  is  more  efficient  during  execution  in  terms  of  time 
required  and  the  modular  structure  has  also  given  the  desired  benefits  at  the 
program  development  stage. 

The  facilities  of  SESAME  are  described  as  they  exist  at  the  time  of  publi¬ 
cation  of  this  Report.  Extensions  and  enhancements  may  easily  be  added  as  new 
needs  arise. 

Furthermore,  the  routines  as  described  in  this  Report  are  not  specific 
either  to  real-time  simulation  or  to  the  computer  on  which  they  have  been  devel¬ 
oped.  They  are  written  in  Fortran  and,  supplemented  by  appropriate  integration 
techniques,  could 
aircraft  dynamics 


be  used  on  other  computers  for  general  purpose  calculations  of 
in  six  degrees  of  freedom. 
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DESCRIPTION  OF 

INDIVIDUAL  SESAME  ROUTINES 

Section 

Routine 

Source  file 

Binary  file  (ROM) 

A.  I 

SACCBOD 

SACBS 

SACBR 

A. 2 

SACCLIN 

SACLS 

SACLR 

A.  3 

SACCROT 

SACRS 

SACRR 

A. 4 

SALFBET 

SALFS 

SALFR 

A. 5 

SCOUNT 

SCNTS 

SCNTR 

A. 6 

SDCOS 

SDCSS 

SDCSR 

A. 7 

SEULER 

SELLS 

SEULR 

A  •  8 

SILS 

SILSS 

SILSR 

A. 9 

SINIT 

SINTS 

SINTR 

A.  10 

SPATH 

SPTHS 

SPTHR 

A.  1  1 

STV 

STWS 

STWR 

A.  12 

SVELOC 1 

SVLIS 

SVLIR 

A.  13 

SVELOC2 

SVL2S 

SVL2R 

A.  14 

SWIND 

SWNDS 

SWNDR 

A. 15 

SYSCOM 

SYSCS 

SYSCR 

A. 16 

DSCRT1 

SCR  IS 

SCR  1  R 

A. 16 

DSCRT2 

SCR2S 

SCR2R 

A.  17 

ISDSCR 

ISDCS 

ISDCR 

AO 
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A,  1.2  Call 
A.  I  .  3  Inputs 


A.  !  Subrout ine  SACCBOD 

A.  1 . 1  Purpose  Calculates  accelerations,  in  body  axes,  of  aircraft  centre 

of  gravity  (eg),  pilot  and  various  accelerometer  locations. 

A,  1.2  Call  CALL  SACCBOD 

A.  1 .3  Inputs  P,Q,R 

l’DOT ,  QDOT,  ROOT 
RG,  W 

FTX,  FTY,  FTZ 
S 1  3 ,  S23,  S33 
XP,  ZP 

XI,  Yl,  Zl;  X2 ,  Y2 ,  Z2;  X3,  Y3,  Z3;  XA,  YA ,  ZA ;  X5,  Y5,  Z5 

A.l.A  Outputs  AXCG,  AYCG,  AZCG 

AXACG,  AYACG ,  AZACC. 

AXP,  AYP,  AZP 
AZAP,  AZCG I 

AY  I,  AZ2,  AX3 ,  AYA ,  AZ5 
A.  1.5  Data  required  XP,  ZP 

XI,  Yl,  Zl;  X2 ,  Y2 ,  Z2;  X3,  Y3,  Z3;  XA ,  Y4 ,  ZA ;  X5,  Y5,  Z5 
A.  1.6  Description  (NB  all  accelerations  in  units  of  g  )  -  see  also  section  3.9 


A.l.A  Outputs 


A.  1.6. I  Specific  accelerations  at  aircraft  centre  of  gravity  (as  measured 
by  an  accelerometer). 

a  ■  F_v/W 

xcg  TX 

a  *  F  /W 

yeg  TY 

a  -  1T7/W 

zcg  TZ 

Note  that  in  steady  level  flight,  with  0  =  0  ,  a  'normal'  accelerometer 

at  the  centre  of  gravity  registers  -I  g  and  in  a  pull-up,  for  example,  registers 
-1.5  g. 

A.  1.6. 2  Absolute  accelerations  at  aircraft  centre  of  gravity 


a 

ss 

a 

■f 

s  . 

xacg 

xcg 

1  3 

a 

m 

a 

•f 

yacg 

yeg 

23 

a 

m 

a 

4- 

zacg 

zcg 

33 

are  direction  cosines. 

*9 

Sl  3  ' 

-  sin  6 
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A.  1.6.3  Accelerations  at  pilot 


xp 


yp 


r.p 


axcg  -  [yq2  +  r2)  ■  ypr  +  q)]A 

aycg  *  [ypq  +  r*}  +  Zp(qr  ■  ^]/8 
a2cg  *  [XP(pr  “  q)  ~  VP‘  *  q2)]/li 


A. 1.6.4  Absolute  normal  acceleration  at  pilot  (or  'manoeuvre  | 
motion  cues) 


a  -a  ♦  S 
zap  zp  33 

A. 1.6. 5  Acceleration  relative  to  1  g  for  recording 


zcgl 


a  ♦1.0 
zcg 


A. 1.6. 6  Acceleration  for  slip  ball  (located  at  Xj,V|,Zj) 

ay,  °  aycg  *  [xl(pq  "  "  yl(p"  *  *  Vqr  "  p) 

A.  1.6.7  Accelerations  for  g  meter  (located  at  x,,y,,z,) 

az  ”  azcg  *  [x: (pr  “  q)  +  y2(qr  ”  p)  ”  7-2(p"  +  q  ) 
A. 1.6.8  Accelerations  for  three  independent  accelerometers 


A 

A 


xcg 


ycg 


-  j^x.3(p'  ♦  r")  -  y3(pq  ~  t)  ~  z3(pr  ♦ 

♦  |"x^(pq  ♦  f)  -  v^(p"  ♦  r“)  ♦  z,  (qr  -  p)j^g 
f  ♦  j^x5(pr  -  i'J)  +  (qr  ♦  \\)  -  z^p-  ♦  q")^g 


A.  1 .  7  Init ial isation 


Nil 


A.  1.8  Subroutines  used 


Nil 

A . I . 9  Remarks 

9 

A.  1.9.1  See  F.tkin  (Dynamics  of  Atmospheric  Flight)  pp  122-12 
tion  of  acceleration  at  an  arbitrar>  point. 


9  I 


g'  used  for 


for  deriva- 
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A.  1.9. 2  'Flat  earth'  is  assumed. 


A.  2  Subroutine 

A. 2. 1  Purpose 
A. 2. 2  Call 
A .  2 . 3  Inputs 

A. 2. 4  Outputs 


SACCLIN 

To  compute  linear  accelerations  in  earth  axes 

CALL  SACCLIN 

FTX,  FTY,  FTZ 
RXMASS,  G 

SU  -  S33 

FTN,  FTE,  FTD 
VKNDOT,  VKEDOT,  VKDDOT 


A. 2. 5  Description 

A. 2. 5.1  Transforms  total  forces  in  body  axes  (FTX,  FTY,  FTZ)  to  earth 
axes  vFTN,  FTE,  FTD),  using  direction  cosines  S^  etc. 


TN 


TE 


TD 


S1 1FTX  +  S21FTY  +  S3IFTZ 


S12FTX  +  S22FTY  *  S32FTZ 


S13FTX  +  S23FTY  +  S33FTZ 


F  /m 

KN 

TN 

m 

F__/m 

KE 

TE 

KD 

FTD/m  +  8 

A. 2. 5. 2  Calculates  accelerations  in  earth  axes 

V. 

A. 2. 6  Initialisation 
Nil 

A. 2. 7  Subroutines  used 
Nil 

A. 2. 8  Remarks 

User  must  supply  total  forces  in  body  axes  (FTX,  FTY,  FTZ). 

A . 3  Subroutine  SACCROT 

A. 3. I  Purpose  Calculates  angular  accelerations  in  body  axes 

A. 3. 2  Call  CALL  SACCROT 
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A. 3. 3  Inputs  XLLTOT,  XMMTOT,  XNNTOT 

P.  Q,  R 

CII  -  CI10 

A. 3.4  Outputs  PDOT,  QDOT,  ROOT 

A. 3. 5  Description 

Given  total  moments,  angular  velocity  components  and  inertia  coefficients, 
calculates  angular  acceleration  components  in  body  axes. 


• 

P 

-  C^L  ♦  CI2N  ♦  CI3pq  ♦  Cl4qr 

q 

-  CI5M  ♦  CI^(r2  -  p2)  +  CI?rp 

f 

-  CIgN  +  CI2L  ♦  CIgPq  ♦  CIj0qr 

A. 3.6 

Initialisation 

Nil 

A. 3. 7 

Subroutines  used 

Nil 

A. 3. 8 

Remarks 

A.3.8.1  CII  etc 

are  calculated  in  SINIT  from  inertia  information  supplied 

by  the 

user. 

A. 3.8.2  The  user  must  provide  the  total  moments  in  body  axes  (XLLTOT, 

XMMTOT 

,  XNNTOT). 

A. 3. 8. 3  An  xz 

plane  of  symmetry  is  assumed. 

A. 4 

Subroutine 

SALFBET 

A. 4. 1 

Purpose 

Calculates  angles  of  attack  and  sideslip 

A. 4. 2 

Call 

CALL  SALFBET (DT) 

A. 4. 3 

Inputs 

UB,  VB,  WB,  RADTOD 

A. 4. 4 

Outputs 

ALFAR,  ALFAD,  SALFA,  CALFA,  ALFADOT, 

BETAR,  BETAD,  SBETA,  CBETA,  BETADOT 

A. 4. 5 

Argument 

DT  is  frame  time  of  loop  in  which  SALFBET  is  used. 

A. 4. 6 

Description 

A. 4. 6.1  Given  the  body  axes  components  of  airspeed,  calculates  a  and  8 
4 

from  (see  Hopkin  section  6.2) 
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a  -  tan  1  WB/UB,  with  a  taking  the  sign  of  WB  and  being  in  the  range 
-it  to  +n  . 

0  -  tan"'  VB/(UB~  ♦  WB~V 

being  in  the  range  -tt/2  to  +n/2  ,  and  taking  the  sign  of  VB. 

A. 4. 6. 2  Calculates  the  sine  and  cosine  of  a  and  0  . 

A. 4. 6. 3  Converts  a  and  0  to  degrees  and  stores. 

A. 4. 6. 4  Estimates  the  rate  of  change  of  the  incidence  angles  (a  and  6) 
from  the  present  and  two  preceding  values,  eg 

a  -  (a.,  -  4aj  +  3a)/(2DT) 


where  DT  is  the  time  step. 


A. 4. 7  Initialisation 


During  the  first  entry  to  the  routine  the  two  preceding  values  of  a  and 

•  •  •  • 

r*  used  xn  estimating  a,  0  are  set  to  the  current  values  of  a  and  0  ,  thus 

•  • 

making  a  =  0  =  0.0  for  two  time  steps  after  'compute'. 

A. 4. 8  Subroutines  used 

Library  mathematical  functions. 

A. 4. 9  Remarks 


For  UB  -  WB  -  0.0 

,  a  »  0.0. 

If  VB  is  zero  as 

well,  then  8  «  0.0. 

A. 5 

Subroutine 

SCOUNT 

A.  5. ! 

Purpose 

Resets  pass  count  for  derivative 

A. 5. 2 

Call 

CALL  SCOUNT 

\.  ?  .3 

Inputs 

JCOMP  (from  COMMON  'CONTROL') 

A. 5. 4 

Outputs 

NCPASS,  JJCOMP,  NRUN 

A. 5 .5 

Description 

section 


A. 5. 5. 1  Resets  pass  count  NCPASS  to  zero.  NCPASS  is  used  in  the  derivative 
section  to  keep  track  of  sub-steps  of  integration  routine. 


A. 5. 5. 2  Picks  up  JCOMT  from  COMMON  'CONTROL'  and  reassigns  to  JJCOMP  in 
i.OMMON  SYSTEM'  to  enable  JJCOMP  to  be  used  as  an  ' init ial isat ion  complete'  flag. 
JCOMP  (and  hence  JJCOMP)  gets  set  to  1  once  'compute'  is  press,  d. 
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A. 5. 5. 3  Calls  DLRUN  to  pick  up  current  run  number  held  in  mailbox. 

A . 5 . 6  Initial isation 
Nil 

A. 5. 7  Subroutines  used 
DLRUN 

A. 5.8  Remarks 

A. 5. 8.1  Since  DLRUN  accesses  the  foreground  mailbox,  any  program  using 
SCOUNT  or  DLRUN  must  be  run  in  the  foreground. 

A. 5. 8. 2  Run  number  is  only  obtained  if  data  logging  is  not  inhibited. 

A. 5. 8. 3  Because  SCOUNT  is  not  called  during  initialisation,  the  run 
number  is  not  obtained  until  after  ’compute’  is  pressed. 

A. 6  Subroutine  SDCOS 

A. 6. 1  Purpose  Calculates  direction  cosines  from  Euler  angles 

A. 6. 2  Call  CALL  SDCOS 


A. 6. 3  Inputs  PSIR,  THETAR,  PHIR,  RADTOD 

A. 6. 4  Outputs  PSID,  THETAD,  PHID, 

SPSI,  CPSI,  STHETA,  CTHETA,  SPHI,  CPHI 

A.b.5  Description 


A. 6. 5. 1  Calculates  sines  and  cosines  of  the  three  Euler  angles 


4)  -  heading  0-360°  (or  ±180) 
0  -  pitch  ±90° 

4>  -  bank  ±180° 


A. 6. 5. 2  Calculates  nine  direction  cosines  used  in  transformation  from 

.  4 

earth  to  body  axes  and  vice  versa  (see  Hopkin  section  5.6). 


SI  1 

ss 

COS 

0 

cos 

m 

*1 

<2 

S12 

39 

cos 

Q 

sin 

s 

S13 

m 

-  sin 

0 

m 

43 

S21 

m 

sin 

* 

sin 

0 

cos 

-  cos 

♦ 

sin 

0 

39 

S22 

• 

sin 

sin 

0 

sin 

♦ 

+  cos 

♦ 

cos 

* 

m 

S23 

• 

sin 

♦ 

cos 

0 

m 

m3 

S3 1 

a 

cos 

* 

sin 

0 

cos 

sin 

♦ 

sin 

39 

ni 

S32 

- 

cos 

* 

sin 

0 

sin 

* 

-  sin 

♦ 

cos 

m 

n2 

S33 

■ 

cos 

♦ 

cos 

6 

■ 

n3 

)08 
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stores 


A.  6 .6 


A. 6. 7 


A. 6. 8 


turns. 


A. 7 


A. 7. 1 
A. 7. 2 


A. 7. 3 


A. 7. 4 


A. 7.5 


A. 7. 6 


A. 7. 7 


A. 8 


A. 8.  I 


A. 6. 5. 3  Calculates  sec  6  and  tan  6  . 

A. 6. 5. 4  Converts  input  angles  (ip,  0,  <f>)  from  radians  to  degrees  and 


Initialisation 

Nil 

Subroutines  used 

Library  mathematical  functions,  eg  sine ,  cosine 
Remarks 

Makes  no  provision  yet  for  0  =  90°  ,  nor  for  large  angles,  ie  multiple 


Subroutine  SEULER 


Purpose 

Call 

Inputs 


Outputs 


Calculates  Euler  attitude  rates. 
CALL  SEULER 
P,  Q,  R 

SPHI,  CPHI,  TANTHT ,  SECTHT 
RADTOD 

PHIDT,  THETDT,  PSIDT 
PD,  QD,  RD 


Description 

A. 7. 5.1  Calculates  Euler  attitude  rates  from  body  rates  and  Euler  angles 


4>  =  p  +  (q  sin  <j>  +  r  cos  $)  tan  0 

6  =  q  cos  <t>  -  r  sin  $ 

\p  -  (q  sin  <f>  ♦  r  cos  4>)  sec  0 

A. 7. 5. 2  Converts  body  rates  to  degrees/s  and  stores. 
Initialisation 
Nil 


Subroutines  used 

Library  mathematical  functions. 

Subroutine  SILS 

Purpose  To  calculate  glide  slope  and  localiser  indications 

for  conventional  instrument  landing  system  (ILS). 


008 


Appendix  A 


47 


A. 8. 2  Call 


CALL  SILS 


A. 8. 3  Inputs 


JJCOMP,  LSILS,  ILSFLG 
X,  Y,  H,  RADTOD,  DECTOR 


A. 8. 4  Outputs 

A. 8. 3  User  options 

LSILS 

0 

1 

ILSFLG 

1 

2 

3 

A. 8.6  Descript  ion 


EGS,  ELOC,  HSLOPE,  RGS,  RLOC,  RKINK 

XGS,  XLOC,  YLOC,  SGS,  SLOC,  U SLOPE,  BSLOPE,  HKINK, 

CILSI,  CILS2,  CILS3,  CILS4 ,  XKINK 

ILS  on/off  flag 

off,  glide  slope  and  localiser  indications  set  to  zero 
on 

select  beam  type 
3°  straight  beam 
6°  straight  beam 

6°  changing  to  3°  at  height  defined  by  HKINK 


A. 8. 6.1  If  ILS  is  required  (set  by  LSILS),  calculates  conventional  angular 
ILS  beam,  with  no  beam  noise,  assuming  receiver  in  aircraft  is  at  centre  of 
gravity. 


A. 8. 6. 2  If  ILS  is  not  required,  glide  slope  and  localiser  indications  are 
set  to  zero  and  no  other  calculations  are  performed. 

A. 8. 6. 3  Localiser  model 

Beam  origin  is  at  XLOC,  YLOC.  (X  »  0  is  at  runway  threshold,  Y  =  0  on 
runway  centre  line). 

RLOC  is  range  from  transmitter. 

Angular  error  ELOC  is  returned  scaled  to  ±1.0,  ie  as  a  fraction  of  the 
defined  beam  width. 


Beam  width  is  ±SL0C  degrees. 

A. 8. 6. 4  Glide  slope  model 
Glide  path  origin  is  at  XGS. 

RGS  is  range  from  transmitter. 

HSLOPE  is  height  of  beam  centre  line  at  range  RGS. 

Angular  error  EGS  is  returned  scaled  to  ±1.0,  ie  as  a  fraction  of  the 
defined  beam  width. 
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Beam  width  is  +  SGS  degrees. 

USLOPE  is  upper  segment  slope  (if  applicable)  in  degrees. 

BSLOPE  is  lower  segment  slope  in  degrees. 

For  a  single  segment  beam,  USLOPE  -  BSLOPE. 

If  a  2-segment  slope  is  used  (selected  by  ILSFLG),  the  kink  occurs  at  the 
defined  height  HKINK. 

Range  and  position  from  kink  (RKINK,  XKINK)  are  derived, 

A. 8. 7  Initialisation 

A. 8. 7. I  During  initialisation  (when  JJCOMP  ■  0)  beam  constants  are  set  up 
according  to  the  choice  of  beam  defined  by  ILSFLG,  Initialisation  occurs  regard¬ 
less  of  the  state  of  the  ILS  on/off  flag. 

A . 8 . 8  Subroutines  used 

BOUND 

Library  mathematical  routines. 

A. 8. 9  Remarks 

A. 8. 9.1  SILS  is  used  by  SINIT,  if  requested,  to  place  the  aircraft  exactly 
on  the  beam  as  an  initial  condition. 


A. 8. 9. 2  If  the  user  wishes  to  create  his  own  ILS  routine  (or  similar  guid¬ 
ance)  he  could  use  the  same  name  and  ensure  that  it  will  interface  correctly  to 


SINIT, 

ie  receive  a 

value  X  and  return  a  value  HSLOPE. 

A. 9 

Subroutine 

SINIT 

A. 9.  1 

Purpose 

Various  initialisation  calculations 

A. 9. 2 

Call 

CALL  SINIT 

A. 9. 3 

Inputs 

XIX,  XIY,  XIZ,  XIZX 

HIC,  XIC,  HSLOPE 
DENR,  SPSNDR,  TEMPR,  PRESSR 
SPSL,  RHOSL,  TEMPSL,  PRESSL 
VWKT0,  PSIWD,  SHRFAC 

PSIDIC,  GAMDIC,  VKTIC,  BETADIC,  ALFADIC 
SWREF,  SPAN,  W,  XCG,  ZCG,  XCGREF,  ZCGREF 
KISA,  ISHR 
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A. 9.4  Outputs  CI1  -  Cl  1 0 

x 

ROOTSIG,  SPSOND,  RHO,  TALPHA,  PALPHA 

CPSIW,  SPSIW,  VWNL0,  VWEL0,  VWN,  VWE,  VWD 

VKNIC,  VKEIC,  VKDIC,  VN,  VE,  VD,  VK,  VT 

THETADIC,  PHIDIC 

SB2,  XMASS,  RXMASS,  DXCG,  DZCG 

PDIC,  QDIC,  RDIC 

RADTOD,  DEGTOR,  FPSTKT,  XK2FPS,  G,  RG 


A. 9. 5  Data  to  be  supplied  by  user 


SWREF,  SPAN,  XCGREF,  ZCGREF,  ZCG,  XIX,  XIY,  XIZ,  XIZX 


A. 9.6  User  options 


KISA 


flag  to  control  variation  of  atmospheric  properties 
with  height 


0 


sea  level  conditions  at  all  times 


1  standard  atmosphere 

-1  fixed  values  appropriate  to  initial  height  (HIC) 

ISHR  wind  shear  selection  flag  (see  SWIND  for  full 

description) 


A. 9. 7  Description 


A. 9. 7.1  Calculates  inertia  coefficients,  according  to  equations  (20) 


CI0 

1.0/(XIZ*XIX  -  XIZX*XIZX) 

CI1 

= 

XIZ*CI0 

CI2 

= 

XIZX*CI0 

CI3 

= 

XIZX* (XIX  -  XIY  +  XIZ)*CI0 

CI4 

= 

(XIZ* (XIY  -  XIZ)  -  XIZX*XIZX)*CI0 

CI5 

= 

1 .0/XIY 

CI6 

= 

XIZX*CI5 

CI7 

= 

(XIZ  -  XIX)*CI5 

CI8 

= 

XIX* CI0 

CI9 

s 

(XIX* (XIX  -  XIY)  +  XIZX*XIZX)*CI0 

Cl  1 0 

= 

-XIZX(XIX  -  XIY  +  XIZ)*CI0 

A. 9. 7. 2  If  HIC  is  set  negative,  calculates  initial  height  to  put  aircraft 
exactly  on  ILS  beam. 

A. 9. 7. 3  Sets  up  initial  values  of  atmospheric  properties.  If  KISA  =  0, 
sets  up  sea  level  conditions.  Otherwise  uses  SATMOS  to  derive  properties  appro¬ 
priate  to  HIC. 

A. 9. 7.4  Calculates  local  true  wind  components  from  datum  wind  speed 
(VWKT0J  and  direction  Olfy) 
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VWNL0  “  ”  VWKT(J  COS 
VWEL0  “  ”  VWKT0  Sln 

and  allows  for  wind  shear  to  give  initial  values  for  V._,,  V,_. 

WN  Wh 

A. 9. 7. 5  If  initial  condition  is  not  at  rest  derives  velocity  relative 
to  ground,  and  then  components  for  IC  values,  given  wind  speed  and  direction, 
desired  initial  airspeed  and  heading.  (Note  that  VK  in  this  routine  is  resultant 
velocity  relative  to  ground,  whereas  VK  obtained  in  SPATH  is  truly  the  ground 
speed,  ie  in  the  horizontal  plane.) 

An  iteration  technique  is  used  to  ensure  that  the  initial  track  and  head¬ 
ing  are  compatible  with  the  specified  wind.  The  initial  heading  will  thus  not 
reveal  to  the  pilot  anything  about  the  applied  wind. 

A. 9. 7. 6  If  initial  condition  is  at  rest  (defined  by  VKTIC  <  1.0)  then 
ground  speed  is  set  to  zero. 

A. 9. 7. 7  An  initial  value  for  9  is  calculated,  assuming  0  =  0,  from 

0  *  Y  +  a  . 

A. 9. 7. 8  Aircraft  related  constants  are  calculated  from  data  supplied  by 
the  user 


SB2  =  SWREF*SPAN*SPAN 
XMASS  =  W/G 
RXMASS  =  1.0/XMASS 
DXCG  =  XCG  -  XCGREF 
DZCG  =  ZCG  -  ZCGREF 


A. 9. 7. 9  Various  system  constants  are  given  default  values  via  DATA  state¬ 
ments,  as  detailed  below. 

Default  values 

RADTOD  57.295780 

DEGTOR  0.0174533 

FPSTKT  0.5921053 

XK2FPS  1.688889 

G  32.174 

RH0SL  0.0023769 

SPSL  1116.45 

RG  0.03108 

TEMPSL  288.15 

PRESSL  14.69597 

DENR  1.0 

SPSNDR  1.0 

TEMPR  1 . 0 

PRESSR  1.0 

PDIC  0.0 
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QDIC 

0.0 

RDIC 

0.0 

PHIDIC 

0.0 

KISA 

0 

ISHR 

1 

A. 9. 7. 10  A  software  switch  is  also  set  on  (IHSH10  =  1 )  to  enable  the  high 
speed  hybrid  input/output .  It  is  communicated  via  labelled  COMMON  HSHIO. 

A. 9. 8  Initialisation 


Whole  routine  is  for  initialisation  only. 

A. 9. 9  Subroutines  used 

SATMOS,  S1LS,  WSREAR  and  library  functions. 
A. 9. 10  Remarks 


Will  not  initialise  correctly  if  initial  *  is  non-zero. 
A. 10  Subroutine  SPATH 


A.  10.1  Purpose  Calculates  velocity  over  ground  and  flight  path  angles 

A.  10.2  Call  CALL  SPATH 


A, 10.3  Inputs 


A. 10.4  Outputs 


A. 10.5  Descript  ion 


VKN,  VKE,  VKD 
RADTOD,  FPSTKT 

VK,  VKKT, 
PSIKR,  PSIKD, 
GAMMAR,  GAMMAD 


A. 10.5.1  Calculates  velocity  over  ground  (ground  speed) 


VK  -  (VKN2  +  VKE2)* 


A. 10.5.2  Converts  VK  to  knots  and  stores. 
A. 10,5.3  Calculates  climb  angle 


Y  =  tan~ 1 (-VKD/VK) 

in  the  range  -n/2  to  +tt/2  . 

A, 10.5.4  Calculates  track 

X(*  V  "  ta,r'(VKE/VKN) 


in  the  range  0  to  2n 
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A. 10.5.5  Converts  y  and  x  to  degrees  and  stores. 
A.  10.6  Initialisation 
None 

A .  1 0 . 7  Subroutines  used 

Library  mathematical  functions. 

A. 10.8  Remarks 


A.  10.8. 1 

If 

VKD  = 

VK  =  0.0  ,  y  =  0.0  . 

A.  10.8.2 

If 

VKE  = 

VKN  =  0.0  ,  x  =  0.0  . 

A. 10.8.3 

For 

VKE 

=  0.0  ,  x  =  0  for  V^j  positive  but  x  =  180°  for 

Vjqj  negative. 

A.  1  1 

Subroutine 

STV 

A.  1  1 .  1 

Purpose 

Computes  TV  positions  and  rates,  and 

controls  belt 

A.  1  1 .2 

Call 

CALL  STV 

A.  1  1 .3 

Inputs 

VKN,  VKE,  VKD,  VSHIPKT,  XK2FPS 

X,  Y,  H,  XIC 

XP,  ZP,  CTHETA,  STHETA,  CPSI,  SPSI, 
JJCOMP ,  NTVB 

LFWD,  LBACK,  LYTVIC,  LXTVIC,  LCYCLE 

PSIDT, 

THETDT 

A.  1  1 .4 

Outputs 

XTV,  YTV,  HTV,  XIC,  XIHX,  XIHY 

XDOTM,  XDTV ,  YDTV ,  HDTV 

SXPL,  SXMIN,  SXPLUS ,  SXMINUS,  SYPL, 
XDTMX 

IDSYNC 

SYMIN, 

HTVLIM, 

A.  1  1  .5 

Data  to  be 

supplied  by  user 

NTVB  TV  belt  flag 

1  700: 1  model 

2  2000:1  model 

3  5000:1  model 
VSHIPKT  ship  velocity,  knots. 

A. 11.6  Description 

A. II. 6.1  Initialisation. 


L_ 


A.  11.6.2  Belt  slew  control. 

If  forward  slew  is  selected,  X  rate  is  forced  to  maximum  positive  value 
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If  backward  slew  is  selected,  X  rate  is  forced  to  maximum  negative  value, 


In  both  cases,  the  aircraft's  position  keeps  in  step.  Positional  calcula¬ 
tions  can  also  be  referred  to  a  ship  moving  at  constant  velocity  VSHIPKT. 


A. I  1.6. 3  Integration  control  for  X  and  Y  (belt  logic). 


If  Y  reset  is  selected,  Y  position  resets  to  YIC. 


If  X  reset  is  selected,  X  position  resets  to  XIC  as  originally 


defined  in  ICEILE. 


If  X  reset  is  not  selected,  additional  logic  is  invoked  to 


control  the  belt  (x) : 


if  X  is  not  near  the  belt  join  (defined  by  SXMIN,  SXPL  -  see 
belt  diagram),  the  program  continues  to  the  position  module. 
If  X  is  near  the  join  but  ' LCYCLE'  is  not  set,  the  program 
continues  to  the  position  module. 

If  'LCYCLE'  is  not  set,  the  program  continues  to  the  ■ osition 


module . 


If  'LCYCLE'  is  set,  X  is  forced  to  cycle  between  SXMINUS  and 
SXPLUS.  IDSYNC  is  set  to  1  to  uncouple  the  belt  position 


feedback . 


If  'LCYCLE'  is  set,  but  the  demanded  X  position  is  less  than 
the  negative  cycle  limit  SXMINUS  (as  can  happen  if  X  had  been 
previously  reset),  the  X  integrator  is  reset  and  XIC  is  set  to 


SXMINUS.  This  will  force  the  belt  to  start  to  move. 


If  X  is  between  SXMINUS  and  SXPLUS,  the  program  continues  to 
the  position  module. 

When  X  reaches  SXPLUS,  it  resets  to  SXMINUS,  before  going  on 
to  the  position  module. 


A.  I  1.6. A  Position  module  -  TV  positions  and  rates  in  three  axes  (equa¬ 
tions  (22)  to  (25)). 


Calculates  pilot  position,  and  rate  of  change  of  his  position, 
allowing  for  his  location  away  from  the  aircraft's  centre  of 


gravity . 


If  TV  belt  is  trying  to  position  before  minimum  X  possible  on 


the  belt  (SXMIN),  or  if  X  reset  is  demanded,  then  X.^,  is  set  to 


zero  to  stop  hunting.  This  is  not  done  if  LCYCLE  is  set. 


If  H  is  above  the  ceiling  for  the  model,  h,^,  is  set  to  zero 


to  stop  bouncing. 
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(d)  Limits  are  finally  applied  to  x.^,,  y^v  and  h.^  to  keep  within 

model  confines.  These  limits  do  not  affect  computed  aircraft 
position  and  navigation. 

A .  I  1 . 7  Ini t ial i sat  ion 

A.  I  1.7.  I  During  ini t ial isat ion,  the  set  of  TV  constants  appropriate  to 
the  particular  belt  are  selected  from  the  array  TVCON(8,3).  These  constants 
define  the  limit  of  movement  as  follows: 


SXPL 

TVCON 

(1. 

NTVB) 

Maximum  value  of  x^.^,  allowed,  also  determines 
when  'near  belt  join'. 

SXMIN 

TVCON 

(2, 

NTVB) 

Minimum  value  of  Xjy  allowed,  also  determines 
when  'near  belt  join'. 

SXPLUS 

TVCON 

(3, 

NTVB)  ] 

x  cycles  between  these  two  values  if 

SXMINUS 

TVCON 

(4. 

NTVB) J 

'cycle'  mode  is  selected 

SYPL 

TVCON 

(5, 

NTVB) 

Maximum  value  of  y^,  allowed 

SYMIN 

TVCON 

(b. 

NTVB) 

Minimum  value  of  y_,^,  allowed 

HTVLIM 

TVCON 

(7, 

NTVB) 

Ceiling  for  TV  height 

XDTMX 

TVCON 

(8, 

NTVB) 

Belt  velocity  when  slewing. 

The  x-related  quantities  are  illustrated  in  Kig  A4. 

STV  holds  three  values  of  each  parameter  in  TVCON,  the  appropriate  one 


being  selected  according  to  the  TV  belt  in  use.  Because  of  this,  on-line  perma¬ 
nent  changes  to  SXPL  etc  can  only  bo  achieved  by  changing  the  associated  element 
of  TVCON. 

A. II. 7. 2  The  initial  value  of  X  obtained  from  the  1C  file  is  stored  in 
a  temporary  location  (XICF)  so  that  a  demand  to  reset  X  will  reset  X  to  this 
value  regardless  of  whether  cycle  has  previously  been  selected. 

A. I  1.8  Subroutines  used 

BOUND 

A . I  1 . 9  Remarks 

A.  I  1.9.1  The  array  of  TV  constants  is  accessible  to  1’VIOO  and  so  can  be 
overwritten  if  new  values  are  found  to  be  more  appropriate  to  the  hardware. 
Alterations  can  also  be  made  via  the  semi -permanent  changes  file  (see  SYSCOM) . 
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A. 11.9.2 
control  desk  and 


The  logical  controls  are  set  via  special  buttons  etc  on  the 
so  are  hardware  dependent . 


LCYCLE 

desk 

switch  5 

LXTVIC 

desk 

switch  6 

LYTVIC 

desk 

switch  7 

LEWD 

slew 

toggle  patched 

LIVACK 

slew 

toggle  patched 

1DSYNC 

software  flag  contro 

o  ADA  sense  line  DCS  03  (line  A) 
o  ADA  sense  line  DGS04  (line  5) 
s  ADA  control  line  DOC  I  7  (line  16) 


A. 11.9.3  To  start  the  aircraft  at  some  position  outside  the  limits  of 
the  TV  model,  set  up  X1C,  Y1C  to  suit.  When  the  aircraft  position  comes  within 
the  confines  of  the  belt,  it  will  begin  to  move  with  the  aircraft.  X  and  Y  posi¬ 
tions  can  be  reset  using  their  appropriate  buttons  without  having  to  reset  the 
whole  computation. 


A. I  1.9. A  If,  after  a  reset  to  some  distant  location,  immediate  movement 
of  the  belt  is  required,  put  'X  reset'  off  and  'cycle'  on. 

A. I  1.9.5  No  provision  has  been  made  for  large  bank  angles  (>l  rev),  nor 
for  cloud  or  visibility  control. 

A. 11.9.6  Successful  use  of  this  routine  to  drive  the  TV  model  is  depend¬ 
ent  on  analogue  computer  circuitry  incorporating  position  feedbacks. 

A. I  1.9.7  Do  not  initialise  with  'CYCLE'  on  if,  subsequently  in  the  run, 
you  may  want  to  reset  X  outside  the  confines  of  the  belt.  This  is  because 
initialisation  occurs  twice  and  so  XIC  is  overwritten. 


A.  12 

Subroutine 

SVELOC 1 

A. 12.  1 

Purpose 

Calculates  body  axes  velocity  component 

A.  12.2 

Call 

CALL  SVELOC 1 

A.  12.3 

Inputs 

VKN,  VKE ,  VKD 

VWN ,  VWE,  VWD 

SM  “  S33 

A.  12. A 

Outputs 

UB,  VB,  WB 

VN,  VE,  VD 

A.  12.5 

Descript  ion 

A. 12.5.)  Given  the  components  of  velocity  relative  to  ground  (VKN,  VKK, 
VKD)  and  the  components  of  total  wind  speed,  including  turbulence  (V'WN,  VWK,  VWD) , 
calculates  components  of  velocity  relative  to  the  air  in  earth  axes  (VN,  VK,  VD) , 
then  resolves  to  body  axes. 


5b 


Append i x  A 


A . 1 2 . b  Init ial isat i on 
None 

A .  12.7  Subroutines  used 
None 

A.  12.8  Remarks 

A. 12.8.  I  Assumes  wind  velocity  component  is  positive  in  the  name  direc¬ 
tion  as  the  component  relative  to  the  air.  Thus,  for  a  head  wind  and  for  the 
aircraft  on  a  northerly  heading  (t|i  =*  0)  ,  VWN  is  negative. 


A. 12.8.2  VN, 

VF,  VD  are  initially  calculated  in  S1NIT. 

A.  13 

Subrout ine 

SVEL0C2 

A. 1 3. 1 

Purpose 

Calculates  total  airspeed,  equivalent  airspeed.  Mao' 
number  etc  using  appropriate  atmospheric  properties 

A .  1  3 . 2 

Call 

CALL  SVEL0C2 

A. 13.3 

Inputs 

UB,  VB,  WB 

11 

SWRKK ,  CREK,  SPAN,  SB2 

KISA 

KPSTKT ,  RHOSL,  SPSL,  TEMPSL ,  PRESSL 

A. 13. A 

Output  s 

VT,  VTKT ,  VEAS,  VK.ASKT ,  XMACH 

QllYN,  HLFROV ,  QDYNS,  QSCREF,  QSSPAN,  QSB2IY 

ROOTS IC,  RHO ,  SPSOND,  DENR,  SPSNDR,  TEMPR ,  PRFSSR 
TRATIO,  A  RAT  10,  DRATIO,  PRATTO,  TALPHA,  PALl’HA 

A . 13.5  Data  to  be  supplied  by  user 

Note  that  certain  geometric  parameters  for  the  aircraft  must  be  supplied 
bv  the  user. 


SWRKK 

CREF 

SPAN 


reference  wing  area 
reference  chord 
reference  wing  span 


SB2  is  a  combination  (SWREF*SPAN*SPAN)  calculated  in  S1N1T. 

A .  13.6  User  options 

The  flag  KISA  selects  the  way  in  which  atmospheric  properties  varv.  It 
must  only  be  changed  in  INITIAL,  before  SIN1T  is  executed. 

KISA  -  1  standard  atmosphere 

*  0  fixed  values  appropriate  to  sea  level 

■  -I  fixed  values  appropriate  to  initial  height  MIC 


Otis 
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A. 13.7  De. cription 

A.  13.7. I  Calculates  total  airspeed,  converts  to  knots  and  stores. 


/  2  2  ,J\- 

^*B  +  V13  WbJ 


A.  13.7.2  If  the  flag  KISA  is  set  to  1,  obtains,  for  the  current  height, 
new  properties  (density  ratio,  speed  of  sound  ratio,  temperature  ratio  and  pres¬ 
sure  ratio)  of  a  standard  atmosphere,  using  the  routine  SATMOS. 

If  KISA  is  not  set  to  1,  atmospheric  properties  are  kept  constant,  as  set 
up  by  SINIT. 

A. 13.7.3  Calculates  properties  of  the  atmosphere 


air  density 

P  =  pslo 

speed  of  sound 

a  =  aSLar 

ambient  temperature 

Ta  "  TSLTr 

ambient  pressure 

Pc  *  PSLPr 

Calculates  equivalent 

airspeed,  converts  to  knots  and 

V 

EAS 

=  VTa*  . 

Calculates  dynamic  pressure  and  related  terms 

QDYN 

Jpv^ 

QDYNS 

>»VTSW 

QSCREF  » 

KVref 

QSSPAN  = 

Kv 

QSB21V  = 

!»VV2) 

Calculates  Mach  number 

M  = 

Vm/a 

T 

Calculates  compressible  adiabetic  flow  relationships 

initialisation  purposes,  these  ratios  are  given  datum 

T 

ratio 

(1  +  0.2M2)'1 
<Tr..io>! 

a 

ratio 

0  .  “ 
ratio 

(a  .  )5 

ratio 

P  .  • 

ratio 

a  .  T  . 

ratio  ratio 

58 


Appendix  A 


A. 13.8  Initialisation 


In  addition  to  those  items  mentioned  above,  datum  values  of  atmospheric 
properties  are  provided  in  SINIT  (qv) . 

A. 13.9  Subroutines  used 


SATMOS 

A. 13. 10  Remarks 


Subroutine  ATMOS,  an  alternative  to  SATMOS,  is  available  to  calculate 
atmospheric  properties  for  hot  days.  It  is  not  included  as  an  option,  however, 
as  it  is  very  much  larger  than  SATMOS,  the  subroutine  for  'standard'  days. 


A.  14 


Subroutine 


SWIND 


A.  14. 1  Purpose 
A. 14.2  Call 
A.  14.3  Inputs 


A.  14.4  Outputs 


A.  14.5 


User  options 

LTURB 

0 

I 

ISHR 

1  (default) 

2 
3 

LSEED 

0  (default) 
I 


Controls  turbulence,  wind  and  wind  shear 
CALL  SWIND 
H,  HIC,  VKTIC, 

XK2FPS,  FRAMET1 ,  SFRACG,  SRDECAY 
VWNL0,  VWEL0,  VWDL0,  CPSIW,  SPSIW 
USIG,  VSIG,  WSIG 
ISHR,  JJCOMP ,  NG,  LTURB,  LSEED 

UTURB,  VTURB ,  WTURB 
VWN ,  VWE,  VWD 
VWNL,  VWEL ,  VWDL 
SHRFAC 


Turbulence  switch 

no  turbulence  calculated  or  output 
turbulence  calculated  and  output 
wind  shear  selection 

no  shear  applied.  Returns  immediately  with 
SHRFAC  =  1.0.  Time  overhead  is  minimal. 

ARB  mean  shear  profile  (logarithmic) 

linear  shear  up  to  333  ft,  thereafter  constant 

random  number  seed  selector 

seeds  in  GUSTS  routine  obtained  from  time  of  day 

constant  seeds  enable  a  repeatable  turbulence  sequence 
to  be  obtained 


USIG,  VSIG,  WSIG  rms  settings  for  each  component  (never  exactly  achieved 
it  is  a  random  process!).  Can  be  changed  in  flight. 

VWKT0  datum  wind  speed,  true  (knots) 
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PSIWD  datum  wind  direction  (degrees).  PSIWD  =  0.0  is  a  wind 

from  the  north.  VWKT0  and  PSIWD  are  set  as  initial 
conditions  and  are  used  by  SINIT  (qv)  to  calculate 
VWNL0,  CPSIW  etc. 


For  the  remaining  parameters  (except  for  VWDL)  changes  are  only  effective 
if  applied  during  initialisation. 


VWDL 

SRDECAY 

SFRACG 


NG 


wind  component  in  'down'  sense  (default  =  0.0) 

decay  parameter,  should  not  be  changed  (default  =  0.7) 

intermittency  parameter,  in  range  0.0  to  0.99 
(default  =  0.0).  Should  never  be  set  equal  to  1.0. 

For  minimum  intermittency,  set  SFRACG  =  0.0.  A  value 
of  0.7  is  probably  as  extreme  as  will  ever  be  needed. 
Refer  to  Ref  12  for  discussion. 

Controls  the  three  gust  ramp  lengths  used  in  the 
turbulence  generation  routine  GUSTS1^,  via 
NGUSTS  =  NG*  TIL^i.aj_sp,eed  _(f  t/s2  . 

NGUSTS,  an  important  parameter,  required  by  the  GUSTS 
routine,  defines  the  number  of  gusts  per  second.  Thus, 
if  the  initialisation  speed  is  200  ft/s  (118  kn) , 

NGUSTS  =  NG  =  A  (the  default  value  for  NG)  and  so  the 
shortest  gust  gradient  distance  is  200/A  =  50  ft  and 
the  longest,  determined  in  the  GUSTS  subroutine,  is 
16  times  the  shortest,  or  800  ft.  Initialisation  at 
speeds  much  below  100  kn  may  require  NG  to  be  increased. 
If  the  initial  speed  is  zero,  NGUSTS  takes  a  minimum 
value  of  1 .  A  further  limitation  is  that  NGUSTS  must 
be  less  than  the  repetition  rate  of  the  simulation, 
eg  for  a  frame  time  of  50  ms  (20  solutions/s)  and  NG  =  A, 
the  maximum  initial  speed  must  not  exceed  1000  ft/s 
(592  kn) .  This  constraint  is  checked  and  imposed  in 
SWIND. 


A.1A.6  Description 

(a)  Derives  shear  factor  (SHRFAC)  according  to  height,  using  WSHEAR 

subroutine . 


008 


(b)  Calculates  local  wind  velocity  components  VWNL,  VWEL  at  height, 

including  shear  effects. 

VWNL  =  VWNL0  *  SHRFAC 
VWEL  =  VWEL0  *  SHRFAC 
VWNL0,  VWEL0  are  calculated  in  SINIT. 

(c)  If  turbulence  is  required  (set  by  LTURB) ,  calculates  three  raw 
components  of  turbulence  (using  GUSTS  subroutine).  The  gust 
velocities  obtained  are  scaled  by  nominal  rms  values  for  each 
component  (USIG,  VSIG,  WSIG) ,  set  by  the  user,  to  give  turbul¬ 
ence  components  UTURB,  VTURB,  WTURB . 
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(d)  If  turbulence  is  not  required,  turbulence  components  are  set 
to  zero. 

(e)  The  calculated  turbulence  components  are  added  to  the  mean 

wind  components  to  provide  three  components  of  total  wind 

(VWN,  VWE,  VWD)  in  earth  axes. 

VWN  =  VWNL  -  (UTURB  *  CPS1W  -  VTURB  *  SPS1W) 

VWE  »  VWEL  -  (VTURB  *  CPSIW  +  UTURB  *  SPSIW) 

VWD  =  VWDL  “  WTURB 

(f)  A  flow  chart  is  shown  in  Fig  A9. 

A. 14.7  Initialisation 

During  initialisation,  the  gust  generation  routine  GUSTS  is  called 
regardless  of  whether  the  turbulence  switch  LTURB  is  set  on  or  off.  This  is 
necessary  to  perform  the  required  preliminary  calculations.  However,  the  turbu¬ 
lence  components  returned  are  set  to  zero.  GUSTS  uses  its  own  flag  (IFLAG)  to 
control  initialisation,  but  this  is  tied  to  the  system  flag  JJCOMP  in  SWIND. 

A. 14.8  Subroutines  used 

WSHEAR,  GUSTS,  FDC ,  RANDU,  VSTEP,  SEEDVAL 
Library  mathematical  functions. 

A . 14.9  Remarks 

A.  14.9.1  The  calculated  turbulence  components  are  added  to  the  mean  wind 
components,  which  only  later  become  resolved  into  body  axes.  UTURB  is  along  the 
datum  wind,  a  positive  gust  increasing  the  wind  (crudely,  increasing  the  airspeed); 
VTURB  is  across  the  wind  (a  positive  gust,  crudely,  increasing  the  sideslip); 

WTURB  is  vertical,  positive  up  (to  increase  incidence). 

A. 14.9.2  Wind  shear  and  attenuation  of  vertical  turbulence  with  height 
work  from  altitude  rather  than  height  above  local  terrain. 

A. 14.9.3  Details  of  the  gust  generation  process  are  contained  in  Ref  12. 

A. 14.9.4  Wind  shear  is  only  in  magnitude,  not  in  direction. 

A. 14.9.5  The  turbulence  control,  LTURB, is  set  by  a  desk  switc!'  (currently 
number  8). 

A. 15  Subroutine  SYSCOM 

A.  15.1  Purpose  Sets  up  communication  with  system  variables. 


A. 15.2  Call 


CALL  SYSCOM 
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A.  15.3  Inputs  As  set  up  in  'changes'  file. 

A. 15.4  Outputs  Nil 

A. 15.5  Description 

A.  15.5.1  Includes  in  the  labelled  common  area  SYSTEM  all  the  names  of 
system  variables  so  that  they  are  set  up  in  the  Namelist  table. 

A. 15.5.2  Calls  PV200  in  the  defined  way  to  set  up  the  Namelist  table 
address.  The  defined  way  is 

DIMENSION  NAME  (2) 

DATA  NAME/ ' SYSCOM' / 

NAMELIST 

CALL  PV200(NAME) 

GO  TO  3 

2  INPUT  101 

3  CONTINUE 

A. 15.5.3  Provides  for  'semi-permanent'  changes  to  be  read  from  a  file, 
using  the  self-identified  input  facility 

DATA  ISYS/50/ 

REWIND  I SYS 
INPUT  (ISYS) 

REWIND  ISYS 

A. 15.6  Initialisation 

This  routine  is  only  called  in  the  INITIAL  region,  so  it  is  all 
initialisation. 

A. 15.7  Subroutines  used 

PV200  plus  library  routines.  Note  that  to  test  SYSCOM  in  the  background, 
a  dummy  PV200  should  be  supplied. 

A.  15.8  Remarks 

A. 15.8.1  Unit  50  used  for  inputting  semi-permanent  changes  should  be 
assigned  at  model  load  time  to  a  file  in  the  user's  area,  eg 

: ASSIGN  (F:50,  D3,  SYSCHNG) 

A. 15.8.2  Prior  to  running  the  aircraft  model  program,  the  'changes'  file 
should  be  filled  with  the  desired  changes  in  the  form 

XI  =  23.0 

SWREF  «=  560.0 
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The  concluding  asterisk  (*)  is  essential.  If  no  changes  are  required  the  file 
nuat  contain  an  asterisk  otherwise  a  run-time  error  occurs. 

A. 15.8.3  The  positioning  of  SYSCOM  means  that  IC  values  cannot  be  over¬ 
written  since  they  are  read  after  SYSCOM  is  executed.  However,  changes  introduced 
through  SYSCOM  can  be  altered  on-line  since  retained  changes  are  applied  (by 
PV300)  after  SYSCOM. 


A.  16 

Subroutine 

DSCRTI /DSCRT2 

A.  16.  1 

Purpose 

To  control  execution  of  routines  for  discrete  input/ 
output . 

A. 16.2 

Call 

CALL  DSCRTI  (CALL  DSCRT2) 

A. 16. 3 

Inputs 

IAD4SL,  IAD4CL,  IDSCFL,  IPCOFL 

A.  16.4 

Outputs 

By  calling  other  routines,  the  arrays  ISLR(16), 

ICLR( 16) ,  IP (32 ) ,  IDS (32 )  ,  ICO(32)  are  filled. 

A.  16.5  Data  to  be  supplied  by  user 


Four  flags  must  be  set  to  control  which  of  the  basic  routines,  if  any, 
is  executed  in  which  real  time  loop. 


1AD4SL  Execution  control  flag  for  ADA  sense  lines  (READSLR) 

1AD4CL  Execution  control  flag  for  AD4  control  lines  (SETCLR) 

IDSCFL  Execution  control  flag  for  discretes  in  to  the  Sigma 

(READSCR) 

IPCOFL  Execution  control  flag  for  discretes  out  from  Sigma 

(SETDSCR) 

Each  flag  can  take  the  values: 


0  Do  not  execute 

1  Execute  in  DSCRTI,  ic  in  the  fastest  loop 

2  Execute  in  DSCRT2,  ie  in  the  second  loop. 

\ 

Default  values  (set  in  DSCRTI)  and  typical  values  are 


Flag 

Routine 

Defaul t 

Typical  (possible) 

LAD4SL 

READSLR 

1 

1  (2) 

IAD4CL 

SETCLR 

1 

1  (2) 

IDSCFL 

READSCR 

1 

2 

IPCOFL 

SETDSCR 

0 

2  (1) 

A. 16.6  Description 


Calls  appropriate  routines,  according  to  flag  settings.  Routines  which 
are  called  by  DSCRTI  are  not  then  called  by  DSCRT2,  and  vice  versa. 


008 


Appendix  A 


63 


A. 16.7  Initialisation 
Nil 

A. 16.8  Subroutines  used 

READSLR,  SETCLR,  READSCR,  SETDSCR 

A.  16.9  Remarks 

Execution  timings  are 

SETCLR  0.16ms 

READSLR  0.22  ms 

SETDSCR  0.34  ms 

READSCR  0.87  ms 

Total  1.59  ms 

Hence  a  millisecond  or  more  could  be  saved  in  the  main  loop  by  executing 
READSCR  and  SETDSCR  in  the  second  loop. 

A.  17  Function  subroutine  ISDSCR 

A.  17.1  Purpose  To  read  a  single  specified  Sigma  8  sense  line 

A.  17.2  Call  I  =  ISDSCR  (ISET,  LINE) 

ISDSCR  is  returned  with  a  value  0  or  1 ,  depending  on 
whether  the  specified  line  LINE  is  off  or  on. 

A. 17.3  Description 

This  function  tests  the  status  of  an  individual  sense  line,  either  from 
the  logic  patch  panel  in  rack  N(ISET  *  0)  or  from  the  control  desk  switches 
(ISET  =1). 

A. 17.4  Initialisation 
Nil 

A. 17.5  Subroutines  used 
READSCR 
A.  1 7 .6  Remarks 

This  function  is  intended  for  use  during  the  INITIAL  region  only,  when 
time  is  not  critical.  It  is  not  a  fast  routine. 
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i 
j 
• 
ft 
ft 
7 
0 
9 
10 
1  1 
12 
i  3 

1  % 
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1  6 
i ; 
1  8 

19 

20 

2  1 
22 
23 

2A 

2ft 

2* 

27 

28 

29 

j0 

31 

32 

33 
3 ♦ 
3ft 
3* 
37 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 
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c 

c 
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SUBROUTINE  SACCLlN 


•  «»ITT£N  *t:  D.eLOMCtO 

•  d*u:i?.§-?a 

•  OlSC»I*TlON:  T*  I S  Su**®0T1N£  CalCjLATeS  T*«|  TOTAL 

•  FORCE  COMPONENTS  TP  THt  £a«Tm  Axis  GIVEN  Th£ 

•  Comr0NEnT6  TO  TH|  boot  axes  and  The  DIRECTION 

•  COSINES.  TNt  ACCELERATION  IS  also  CALCULATED 

•  FROM  TN|S  INFORMATION  ANO  THE  A  1 H  a?  T  MA$S 


12-8-7*  ST  ST  fO  Common  £Nla«QEO 

CONMON/srtTlM/A| 1000  »,L I AOO  » 

VERSION  identifier 

OATA  XSACCLlN/l?0*78/ 

E0uIvAlFnCEixSaCClIN,l(?5AII 

•••  VARIABLE  INPUTS  ••• 

#  A  )  19)1# 

# A (  22)1# 

#  A  1  25))# 

#  A  (  1  38  I  >  # 


#81139)) # 
# A ( 2% 0  I  ) 

•••  constants  ••• 


EC  U 1 VAlfNCC) Si  1 

1  I  S2  i 

2  (  S3 1 

3  )  E  T  x 


#  A  I  1  7  )  I  #  I  S 1  2 

#  A l  2C)I#iS22 
# A  t  23 ) ) # ( S3 2 

#  A  I  139 ) I # t ETV 
iRxmaSS  , A  I  1)1 


# A  I  18  I  ) # ) Si  3 
# A  (  21 ) I#  |  S23 

#  a (  2*  I  )  t i S33 

#  A  )  13b  I  )  #  )  ETJ 


•••  VARJABLF  OUTPUTS  ••• 


COjl VALENCE l ETn 
1  IVXNOOT 


♦  4L1 J7 1 1  *  lETE 
.  A  1238  I l # i VKCOOT 


|AU3«M»)TT0 
#  A  1239  I  )  » I vkOOOT 


lAfilu 


EQUIVALENCES  #  A l 1 %5 I  I 

c 

C  CALCULATE  TOTAL  FORCES  Tfl  EARTH  AXES 

ETN  •  si 1»ETX«S21aETv«S31aET2 

ETE  •  S12aETxaS22aETv«S32aET2 

ETO  •  S13aETX*S?3aETv«S33aET2 

C  CALCULATE  ACCELERATION 

vknOOTm  ETNaRXMASS 
VKEOOTm  ETCaRXMaSS 
VXOOOT.  ETOaRXMaiSaU 
RETURN 
ENO- 

)M 161#  33«$ACR6  I 


I  EON 
EXT. 


m Tran  l v#  version  coo 
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1  • 

subroutine  saccrot 
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C 
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• 
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• 
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»• 

c 

• 
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c 
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c 
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c 

• 
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12* 

c 

• 

13« 

c 
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1*. 

c 

lb. 

c 

12-1-71  SVSTC*  CKMnBN  thiA"OCD 

1*. 
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1/. 

c 

1*. 

c 

VERSION 
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c 
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c 
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c 

•••  VARIABLE  INPUTS  ••• 

2*. 

c 
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c 

•••  variable  outputs  ••• 
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c 
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c 
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c 

calculate  anoular  accelerations 

37. 
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3f.  ROOT  • 
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A|.  RETURN 

A*.  (NO 
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C 
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c 
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I* IM *4|*0*l« ANv»vN*| t*l»ll  IV  IV  I 
•4  144  •  AI4NVI *N,*  I 
•  AVt  V  A  i>  •  4I|  I  AN»MA|H*U 
•*4?  A  •  It'*,  H|  f  All 
4*414  •  4** I *1 1  AN  I 
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1  •  I 
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C  C4|4UiAl|  |.|H|4Af|»l  *V  At  V  A  AM*  It  |  I  4,  I  *«|  N  ||  t  |.f  **|  «  ||i(i*  *41 

C  VIM  N|4|  ININ1* 

c 
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C  I  *N4>  PIN  |  M|  |N|  »  |  Ai  |  *A|  |  INI 

c 
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4 
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VIII/* 
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VIPI/* 
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1  • 

SUBROUT  I  Nt  SCOUNT 

2* 

C 

3  • 

c 

•  ••  • 

•  • 

c 

* 

• 

9* 

c 

•  0E6CR IPT I  ON  :  IMIS  KOUTJMt  SETS  FASS  COUNT  Fp« 

• 

ft* 

c 

•  DERIVATIVE  SECTION. 

• 

/• 

c 

•  WRITTEN  bY  ;  b«N* toil inson 

• 

ft* 

c 

•  DATE  :i?*ft-76 

• 

5#* 

c 

•  MODIFIED:  31-5-77  TO  InCLUOE  Ol«UN 

• 

iu* 

c 

• 

• 

1 1  * 

c 

•  •♦• 

w* 

c 

13* 

c 

1ft* 

c 

VERSION  IDENTIFIER 

lb* 

c 

1ft* 

0»T»  KSCOU NT/31057// 

3 1  Ob  7  7 

W* 

tUUl  VALENCE  <  XSCOUM  *H2b9U 

1«* 

c 

l2-ft-7ft  SYSTEM  COMMON  fcNLARQED 

1¥* 

C0NP0N/SYSTEN/4C 1000 l,L( *001 

?o* 

c 

PIC*  OP  JCOMP  (SET  IN  SwA  1  T  » 

21* 

COMMON /CONTROL /JCOMP 

22* 

c 

JCOMP  -  i  AFTER  COMPUTE  PRESSED 

23* 

c 

2ft* 

c 

variable  inputs  ••• 

2b* 

c 

2ft* 

c 

2/* 

c 

•  ••  VAR  I AOLf  OUTPUT*  ••• 

2ft  • 

c 

2$  • 

tool  VALENCE  (NRUN  /Clift** I 

310577 

30* 

l001*»L£NCtlNCP»66  <U  1IIH  JJCONP  ,  L  1  B  1  1 

31  * 

c 

NC  PA  SS  IS  COUNT  OF  NuMbt  H  OF  PASSES  THnCuum  OtRIVATlVE 

SECTI ON 

32* 

c 

txcotro  80  F»H  |N  ONF  INTtOHBTION  SUP.  StT  to  ZF«0  In 

Tm  IS 

33* 

c 

NOOTINt  prion  TO  001N0  lNTt6x»TI0N. 

3ft* 

NC  PA  SS  •  o 

3b* 

c 

SET  J  JC  On  P  #  I N  SYSTEM  CONHON/TO  ENABLE  STATUS  TO  bE 

3ft* 

c 

USE  0/ £ • G*  IN  6 1 LS  «  S  T V 

37* 

JJCOMP  •  JCOMP 

3«* 

C  PICK  OP  RUN  NUMBER  FROM  MAILBOX 

3** 

CALL  DlRUN  (NRUN) 

310577 

•  0* 

RE  TURN 

ftl* 

END 

IAS616N  |h|Sl#03#8l>CSS  I 


fFORTRAN  LB*N§ 

EX  T »  fORTRAN  l  V,  VERSION  £00 
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10. 
11* 
15. 
13. 
1  ft. 

ia. 

1  6. 
17. 

1  ft. 

la. 

id. 

*!• 

ii* 


c  J  • 

?ft. 

fa. 
f  ft. 
it* 
in* 
fa. 

30. 
31  . 
3f  . 
33. 
i«i 
jb. 
3ft. 

3  7. 


C 

c 

c 

c 

c 

c 

c 

t 

c 

c 

c 

c 

c 

c 

c 

c 


Sue  ns 


c 

c 

c 

c 

c 

c 


*"I  tTih  hr  .  - 

l'*H  :  l;.>  •  u,0l-“>  ilt-C 
Dtscnjurjos’.* 

T"!  0J  ncriM,  r  SU!,‘'0‘>nNt 

"Ull  *  .  1  Ct>l«lS,i|J „  “‘•Ljt*’*.S 
SfcC  TutTi  * '■PUT  **Gl£s  Ik.^  t 

c.lcUi.,(.p;k°  '-et.  ,„t 


Le  rs  ,i  i  cut 


'"■‘'“Cl  IDtM,rUl> 


04  7 ft  K  si  c ac 

C“«*0CBS 

—  >«IM(  lK,urs ... 

twul*.UNCl,„s 

•••  V»u1,BLf  „UTPUI#  a< 

,  “''‘‘■‘"Ctiosjo 

l  'SPSI 

3  *  CT*«£Tft 

4  ,s*CTmT 

b  ‘  *if 

6  1  Sf? 

•S3f 


.  ft  I  b  / 


I  J* 


•••  CPNST.S,,  ... 


»'I.(C»SI 

*  1  *  1 » I 

#4<  *  *  1  *  # I S I  3 

#4‘ 

** 1  **  *  >#  I S3  J 


#A| 

;«".,sr„( 

**  com 

'*•  *»".<SU 

*V  «M.,S3, 

#  4  •  f  a  1 1 


.ft  < 

.ft  i 


*> » ». 
1 1  ii# 
!♦»»# 
17||  # 
50  n, 
53  m. 


3ft  < 

3a. 

♦  o. 
M  • 
ftf. 
*J. 

♦  ft. 

♦  a. 

*ft. 

**. 

♦  a. 
bo. 
bi. 
bf . 
b3. 
bft. 
ba. 

bft. 

b/. 

bft. 

ba. 

ftO. 
fti  • 
ftf. 

63. 
ftft. 
ftd » 


s::  -'""'tir1 

PI  •  COSl UH ] u I 

HI*  -  CB“n* 

:  :  ss;;.-« 

'l  •  suHi.sr,,tt 

■  S.nMCThtr.  SI»C»Ni.t(,s, 

sssTBSs? . - 

w^safigsi.. . 

PHJO  •  PnL  ’"‘“’OO 
•<tru*N  p"1n*««otou 

(hn 


ftft 

1‘kSlGN 


•<tru*N 

'"'Ol.UJ.S^utS 


LS.N. 

**Tb  *u*rmN  . 

V‘PSl5h  I0u 
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|  •  »UftRJUT|N|  tlUUN 

rf*  c 

J*  C  ••••••••••••••••• . . 

«•  c  •  • 

ft*  C  •  HN|?1|K  M*  •  ptpiu»  ULO  • 

ft*  C  •  OATl  I  lfA»7ft  • 

/.  c  •  OtfcCRIPTlON  :  THlft  guftROUflNl  C  Al COLA  Tl S  • 

•  *  c  •  Tut  *lMTuOt  Rftlftft  0  I V 1 N  TM|  IOMMON|NH  • 

9*  c  •  0»  T»*t  AsaulAN  VtiOCITV  IN  IH1  MOO*  AMft  • 

0*  C  •  • 

1 1  •  t  . . . . . . . 

i  •  C 

J*  t  |  ff  •••7ft  StSTl"  C  MMWPN  iNLA**atO 

ft*  CfftM0N/SV8T| M/ At  1000 » #L • ftOO » 

ft*  C 

ft*  C  VlNtie*  IP(NT1F|(M 

/*  C 

•  *  PATA  HgluLlR  / 1 f 0 • 7 ft / 

9*  tVJlVAlINCl'KSkUllH  ,117110 

0*  C 

1*  C  •••  VANlAftLl  JNMuTi  »M 

i •  C 

J*  tUOl  VAtlNCit  S7’Ml  *  A  I  llllnCfNl  #  A I  |ftll#(9CClHT  *A<  1MI 

ft*  1  iTANTMT  (  A  l  1  ft  H  «  I  R  #  A  I  4>ft>>,<0  •  A  l  ?7»» 

ft*  f  IN  » A  I  i§  •  I 

ft*  C 

7*  C  •••  VAftIAt*l(  Of  TP*  J  T  ft  ftftft 

ft*  C 

9*  tuul  valincci  Rmiot  #a<  tmc  tot  *ai  jo  i  » » ( *>s  n*  ?  »ai  jim 

0*  |  (RO  #  A I  JfH»tUO  $  A  (  J3M»(ND  »A!  Jft  I  I 

1*  c 

J •  C  ftftft  CONSTANT!!  ••• 

J*  C 

ft*  tUPl  VALINCUHAOTOP  fAllftlll 

ft*  c 

ft*  C  CALCUIATC  ATTlttofft 

7*  RHJOT  •  RftTANTHTftl Oft^RMl ♦NftCKHl I 


3ft  • 

TMCTPT  •  U *C RN 1 *R« MRh I 

3»* 

Rtior  •  ttCTHTftl QvSRmI ftMftC*Hi » 

AP* 

c 

CALCUlATf  ANOULAft  VELOCITY 

ftl* 

RD  •  RftRAPTOP 

ft«* 

UP  •  OftRAOTOO 

ft  J* 

HO  •  MftRAOTOO 

ftft* 

c 

ftft* 

Ml  TURN 

ftft* 

•  NO 

I  ASS  ION  IMI|,0),tlUI  > 

1 7  OR T MAN  L  6#  Nft 

(AT*  7  ON  T  N  AN  1 V  i  VfcNSlON  ( 00 


«  »  ,  00  t  mr  *  #1  ;  %  , 


«  #  ,  t  mf  *  #1  ;  t  ,  t%4 

•  •4  .  *»,'»  ••• 

'  #*  »•♦  <  •  #•  ;».  .  -  •*  ;5; 


*♦//»»  »"'»  > .  v>  ,  •  ,  „:?■»  *.  •  •  Ifc  ;  •  • 

'  •••  «»*?*♦.'  »,•»,**  ••• 

fvi  .«  ifi  #  c.y.  •»  tf<  .  .4  ;m  . 

t  »•'!/•  <#**l9C  #*  •  1 7*  «  • »  »«  ;•*  «i  }7« 

l  •««  N*!***  i  r  v<  *C  ^  ••• 

•  VW  <»  idkK  i«IWI<  imV^K  *•  }79  lt 

I  '  V#%  #«  I  ltd  »  >  •  '  K9C  #4«l«»#»#l-4i4*4  »  I  1  %l  I  I  # 

»  Jl?*!  ft  II*  mill*'  »*«;**•>» 

*  #41  !•*>  UlCIlW  #•  I  1*7 »  » «  I  Clw*i  .AilMi  »» 

•  -  CUM  «i<  |||M 

ta 

%  «•«  Mt 

c 

%-4  »•  .1  I  »  1  |  «  J  H  »  '  !/f  1  T  9*  »  *  I  |  *2  I  I 

C 

i*  •  >  X4*»<io*o  iq*  »« 

!•  >1 *ll««l  v.  I  *0  ?0 
t  M  III  *1  OwfMO 

*  **%•*  VOf  »0  »0 
«#0  » *  90 

Of  *  0»i?  !*»•/# 

C  Ml  '  »U  *»**i  t  /0  I*  f  f  («i  I  M  T  t  •)*  *<M  C0fM*lt1k 
'**  M 

C  UIK4 

C  I  **|H  »UNH«Vf1  010  *1.0*1 

t  »  H4|N  Hu»,t ,4»,*  MU  bu«*t 

t  I  W4|N  *iQMl  ■*! 
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1  LONTlNut 

C  It  V  CONSTANTS  f  OW  MAIN  NU*«AV#J  Otd  ULIOl  *tO*>t 
■US  •  1000*0 

■toe  *11000*0 

TtOC  •  0*0 

SUS  •  0.7b 

*t  oc  •  ? .  o 

0*LP**t  •  J.o 

ttStP'C  ■  j.o 

HR  {NR  ■  0,0 

UO  TO  10 
t  CONTJNut 

C  11*  CONSTANTS  70»i  WAIN  RuNmAT,*  010  UCIOk  StPKt 
■US  ■  1000.0 

■LOC  •  11000.0 

TtOC  •  0.0 

SUS  •  0.70 

*LOC  •  2.0 

UVtfl't  •  4*0 

»StO*k  ■  4.0 

MR  IN*  •  0.0 

UO  TO  10 
3  CONUNut 

C  Its  CONSTANTS  70 Hf  StUMfM  |*/j|  TO  MAIN  N^.Ai 


100* 

■  US 

•  1000.0 

101* 

■tot 

•  uooo.o 

102* 

TtOC 

•  O.o 

10J. 

SUS 

•  0.75 

to*. 

stoc 

•  8*0 

10*. 

UStOHC 

•  4.0 

104* 

SStOHi 

•  3.0 

10/* 

HR  IN* 

•  500*0 

104* 

10  COnTINui 

109* 

cotttcrivt  constants 

110. 

Cltsi 

•  RaDTOD/SOS 

Ill* 

Cltst 

•  NAOTOD/SL0C 

112. 

CItSJ 

•  Tam  WSLflFl*0F0T04» 

Hi* 

Cits* 

•  TAM  UStOFI  •OfGTOKl 

IASS 


1 4 •  **IN*  •  RUS  •  MR!*K/C!t*) 

ft.  C  HA  IN  HOUTINI 
4.  «0  COHTIHUt 

/.  C  HkRt  iHHtOlATttV  IF  |  H  'CO^OTf 
*•  HR  I HR  •  •( R  •  kR | HR  I 

»•  Httt  •  *|R  •  kOS) 

FO*  *»tOC  •  •  •  Rt©C> 

FI*  HStOn*«at*C  IttJ  •  O.S*«RRlN«  ♦  a0S<RR|Nr|  >«cils« 

'*•  tos  •  *0o*0<M.0#l*0»C|ttl*|M  •  HStOFCl/RUSt 

'»•  kto*.  •00UU4H*l,0#»*0#Clt»i*IT  •  V10CI/RL0CI 

F4*  90  COhTIHUC 

FS*  Hi  TURN 

*•  t*»0 

)H  IH |SI#  OJ* 1 1  NTS  I 


|  F  OK  T  N  AN  IS,  NS 

CRT.  F  OH  TH  AH  |V|  Vi»Sl<JN  ( 00 


08 
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I  •  SUftROUTUt  8  I N  I  ? 

t*  c 

3*  C 

4*  C 

8*  C 

*•  c 

>•  c 

».  c 

9.  c 

|0«  c 

II*  c 

1  2.  c 

14*  C 

|4.  C 

lb*  C 

lb*  C 

!/•  C 

18'  C 

18*  t 

20*  C 

21*  C 

C 

24*  C 

24.  C 

2b*  C 

28.  C 

2/*  C 

28*  C 

28*  C 

30.  C 

31*  C 

32*  C 

33*  C  1 1*8*74  8V8TC*  C0HH0N  (NiaROCO 

3%*  C  17*10*77  Cl 3  COttRCCHo 

38*  C  8*  2*78  FLAG  f  OH  MjUH  SPEED  UAC  C0NTNOL  AOOiU 

38.  C  8*8*78  CQHH0N  HSMI0  ENtAHuED  To  2  *0*08 

37*  C0nH0N/8V8TiN/A( 1000 |#L( AOOI 


taMHU*  M  J  O*0Lf>F|tlD 
DATE  17*10*77 

DISCRETION  :  Thu  *U8R0UTlNt  CALCULATE  S 

THt  IMTUl  vALOtb  OF  THk  I  N£  AT  I A  COEFFICIENTS 

GIVEN  II«IV|W«Wl 

THt  INITIAL  VALWl*  OF  VAMIA8LL8  LUTEO  81L0* 

FR0H  THt  ST  AT  £0  INPUT  VA«|A8Lfc8  F  HOH  Tml  uS(R 

081 H  PNOVtOtO  INF OHNATIONJ. 

8CTA01C#AL^AulC«0An0IC.VKT  1C#*1C#V1C# 
MIC#1SMA#V*A!O#PSIa0#PS1r01C#8aBEF#SPAN#A 

INITIAL  VALUL8  ootput:- 

M8L0PC#C^8lH.«P8lM«VMNL0#VNtL0»6HPFAC# 

VPN. VAk#  VAO#  Ot NR#8P8N0R# T4.HPP# PRES8M $ 
R0OT8IG# VT#VA# VANIC# VAklC# VA0IC# VN#Vk# 
VO.PSIOtC. THtT aoic.haotoo.oloton, 
FP8Ti(T#aT0FP||#Q#  RmOSL#  SP8L  #PUIC#QDIC# 
RDIC#Pm1DIC#xHAGS#NXHaSS#SR2 

N.B 

HIC  NtOATlVC  PLALk*  AIRCRAFT  ON  US  Ok  AH 
AT  SPtClFltP  BANQt  *IC 

A  ISA  •  0  (OkfAOLTi  SELECTS  SkA  LEVEL  aTHOSPmENE 

VATIC*  0*0  St  T8  8T8TCH  UP  FON  ZkNd  QWOuNO  8PCE0 
OT  Hf  Wt»  I6L  fet  TS  UP  FOP  Al  ASPE CO  1 1 A8  I •  VA  T I C 


18.  C0iH0N/HfHl0/lH8H|0« tOUP 

10*  C 

fO*  C  VkPSlON  I OC NT  IF  It A 

1*  c 

•2*  OATA  A81N1T  7050878/ 

•3*  kWUl  VALENCE  I  ASIMT  UI283M 

>4.  C 

5*  C  •••  INPUT  VANU8LC8  ••• 

•  8.  C 


/• 

EOUI<<ICNCE< *lX 

,*<  »»  1  l.l  KIT 

#At  49  I  1 # | Xl 2 

# A  1  boil* 

8. 

i « m 

,*  i  oi  ii, lie  a 

0  A  1  98II#(XCUNEF 

#  A  1  99  II# 

9. 

\  ZCQNEF 

, *110011. |SP*N 

#  A  1 103  1  l#(S*HkF 

#  A ( 1 04  1  1  0 

0* 

1  mSlopc 

,*I1T3II,I«CT*0IC 

»*USI  1  1  ,  1  *L  FAOI  C 

• » u  *i  1 1 , 

!• 

IGA HO  1C 

,* <|»JI l,< VKTIC 

#  A  1  1  94  1  1  #  (  X I  C 

, *  <1 *6 1 1 . 

2* 

1  HlC 

,*itt;n,iKco 

#A ( 199 1 1 • ( VARTQ 

>*<202  1  1 » 

4. 

IPS1KD 

#  A  !  203  H  #  t  PS  ID IC 

#41204 I  1  •  < ShhF AC 

#  A ( 228  1  1# 

IN 

lAUMH 

5* 

lUOl VALENCE 

1  18  HR 

#LI  ft  1  I # 1  A  I 8A 

#L  (  111 

4. 

c 

/• 

c 

OUT  NUT  *»*l»l»LC»i 

8. 

c 

9. 

CUUi <*EtNCEI «XN»St 

,A|  |II,|KN*«« 

«•<  IH.lClt 

#  A  1  Jftll# 

o. 

1 

ICI2 

•  A  1  34)I#<CI3 

0  A (  37 1 1 # ( C I 4 

#  A  (  3811# 

I* 

2 

(CIS 

,«i  jan.ict* 

,«<  *UI 1, 1  Cl  7 

•Al  4111# 

2* 

3 

1  Cl  8 

0  A  (  4fll#|CI9 

#AI  43  1  1 # ( Cl  1 0 

#  A (  4411# 

4* 

4 

1  VAN 

0  A  1  ftft  M  # | VAC 

#  A  t  ft4  1  1 • 1 VAO 

>*<  57  l  1 , 

4. 

ft 

1 PH06L 

0  A  I  79H#|8PSL 

#  A  1  80II#(DXCG 

•  AII01  II# 

8* 

4 

I02CG 

. A ( 102  t 1 • i 802 

•  AUOI  1  >*l  X 

#41149)1# 

4. 

7 

1 VAN1C 

*A 1211 II  0 | VAEIC 

#A  (212I|#|  VAO!  C 

#4(213)1# 

7. 

8 

TPHIOIC 

»AI21«H»{  THCTAU1C 

,*  <2|b  i  I , 

8. 

9 

\  VN 

#41217  »»•  |  VI 

#  A  !  2 1 8  I  I  #  1  VD 

#  A  1  2 1  9  I  | 

9. 

c 

o. 

tuul  <*ttNCCl  PC  ic 

»  4 ( 220  M  # | 00  1C 

#4(221  1  \ 0  (ROIC 

•  A  1 222  1  I# 

!• 

1 

1  CHS J  w 

, *122*11. <SP81» 

#  4 ( 22ft  1  >•( VANLO 

#  A (224  I  | # 

»• 

t 

T VAtLO 

#A  1227  II 0 | RHO 

# 4 (  72  M  #( SP80NO 

,*<  TJ 1 1, 

3* 

3 

IP00T6U 

0  A  1  71  1  1#  1  VT 

,* I  *711 » l <« 

#  A  1  Aft  1  |  # 

A. 

4 

1 T»lpn* 

.*<  «»  I  1  •  <  P«.PH» 

#  A  <  84  1) 

»•  C 


I 
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•  ••  CONSTANTS  Pth  l  v  cl)  •  •• 

AWUt  VALINCCt  RaDTOC  »A  l|«l  I  I  •  |OCOTO»  # A < | •* H » < Fpgttf T  •«dOll« 
|  #AI  |%b  )l  •  IDfNN  ill  Mill 

i  tSPSNUW  «AI  .tl  «l>l«|MMiSSN  #  A I  B? I  I  » 

9  tTtHPSL  » A I  /AMilWlSSk  »AI  77ll#  |RU  #AI|A4ll 

INlTULlSt  CONSTANTS 
OAT A  HaOTOO#DAUT0H/9?.29S7«U#U*O17*493/ 

DATA  f PSTKT# MS/0.V»?103J*1 .A  BA SB 9/ 

1)4  T A  G#RmOSl»SPSL/3?*1 7*#0*00tf J749«l||A*«b/ 

DATA  Ho/ 0 • 0 J i OB/ 

DATA  TinPSL»PN(.SSL/?M*|4#  1**49497/ 

DA  TA  DA  NR# SP SNOB « TtNPR# PRC SSh/ |*0# 1*0# 1*0#  l  •  J/ 

DATA  P01C«Q0IC#H0IC#»hI01C/u*u#0*0«U*0#0*U/ 

DATA  MSA/O/ 

DATA  lSHB/i/ 

SAT  FUAD  FOR  Hi GH  9P|tD  HYBRID  OUTPUT 
DATA  1  MS  H 1 0/ 1 / 

CIO  •  I  *0/ < * IZ*» I IZ**x I Z* » 

CalCulata  Initial  valuaS  of  Inartia  coafficiants 
Ull  •  C 1 0* A 1  2 

ci?  •  cio**iz* 

CIJ  •  CIO*  »xlx«xlv*AiZ  IMU> 

Cl*  •  C10*IAIZ*I  XI  Y»vl?|«Kl2lAXUl  > 

CIS  •  1 • 0 / X 1 T 
C 1 4  •  C 1 4*  x 1 2x 
Cl  7  •  C15*  <XI2*X|X  I 
C1B  •  CIO*Xlx 

C 1  9  •  CIO*  I  »  *1  x*xi  Y  I  •  *  lk»xUl»l|  ZX  » 

CI10  •  C I0*XIZ *• (.*1 «*»| r.»l<! 

ADJUST  h  FOR  A/C  ON  04  A R  |N<0» 

if  iric*gt*o*o>  go  to  io 
x  •  Ale 
CAUL  BUS 


MIC  •  mSLOPA 
10  CONT iNuf 

IF  «<1SA*C0*0)  GO  TO  ? A 

SAT  uP  INITIAL  AThoSPhIRIC  VALUtS 
call  S AT  HOS I H| C  «  Oi NM  t SP&nUR# ! A NPH # PH A 8SR I 
{«  ROOTSIG  •  SuRTlOCNR) 

SPSOyO  •  6PSL*SPSNDH 
HHO  •  RHOSLaDANN 
TALPHA  •  TtPPSL*Tt*F'R 
PALPHA  •  PHI  SSL* PR AS 9R 


CALCULATA  LOCAL  iRUt  *|n U  LOPPONANTS  ttiVAN  OATU*  "  1  NO 
SPAIO  ANO  0IRACTI9N 

CPSlR  •  COS(PS|BO*DAQTOR ) 

BPS  I «  •  SlN|PS|NO*DAaTOR I 
YHNLO  •  •VMRT0*X*?FPS*CP*1« 

VBALO  •  •VhnT0*XX2FP9«8PS1« 

allobanca  for  b i n r>  shcah 
call  bShCARi 16hN#hIL#ShNFaL) 

VBN  •  VbNL0*SHRF4C 
Yb  A  •  VBtLD*SHRFAC 
YBO  •  0*0 

•AT  STARTING  value  OF  Track  tHUAL  To  HIAOINv* 

PSIKJIC  •  PS  ID |C 

ANO  8T0RI  OCS  IN AD  haaDInO  FOR  NafAMANCI 
PS101C0  •  PS  ID IC 


CALCULATA  SIN  AND  COS  OF  UA"HA 
CUV  •  COSI GAHDlC*OiUTOR) 

CUM  •  COST psikdic*daoton I 

SUV  •  S!N< UARO|C#riuTOR» 

SUM  •  S1N|PS1N01C*DAGT0P I 
IF  I VKTIC*Ll •! *0  I  GO  TO  <0 
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ib*?. 

IbJ. 

lb*. 

Ibo* 

lb6. 

lb/« 

Jb«. 
!**• 
IbU. 
161  • 
162. 
16  J« 
1  t*» 
16b* 
166* 
16/* 
160* 
16** 
i  /o. 

in* 

!**• 

1/3. 

!/•• 

17b. 

176* 

!>/• 

1?B* 

17*. 

lbo* 
161  • 
lb^* 
lbj* 
lb*. 
Ibb. 
1  8  6  • 
lb/* 
lba* 
lb!#* 


C  it  T  UK  VELOCITIES  FfR  Jt6l«§0  AIRSHktL 
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lib* 

8  llTjKB  #IOSl 

lib. 

C 

11b. 

c 

*0M  CALL  Tp  byiOO 

117* 

DIMENSION  NAME i i | 

1 1  *  • 

oaTa  name/ »svsce«' / 

119. 

nanelIST 

1*0* 

c 

Ok >■  INC  INPUT  ONI?  FOR  StfTEH 

UI* 

OATA  ISvS/bO/ 

Wi» 

c 

Ui* 

c 

Ub* 

CALL  *»V200INAH£  1 

Ub* 

00  TO  3 

Ub. 

?  IM*U  T  I  101  1 

u/* 

3  CONTINUE 

us* 

c 

INPUT  SCNUPfNNANENT  CNANOfS 

U9. 

HlNlNO  ISVS 

130* 

INPUT  ilSVSl 

131* 

*1*1*0  1SVS 

IU« 

NitUMN 

Ui* 

bNO 

I S  TS 


ussjbN  ir  IS1«03#SCN1S  I 


ifoktman  Lb* ns 

EAT*  MINIMA*  lv*  VlMSlON  £  00 
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>0 

id 

i  0 
*1 
it 

sj 


dUdHOUTlSt  CSCMTt 

c 

c  . . . . . . . ••♦••••• 

C  • 

C  •  AJTHPR  :  ».N.TO«t ISkON 

C  •  D AT  t  :  ?!'«•'» 

C  •  FjHSOSl  :  T  H| $  MQutiki  «AnOLCS  UlSCNtlik  IN  FHAHl  |. 


C  •  Th(  FLAGS  UD«»l*  UOaCL#  IOSLFL*  IHCfFl 

C  •  A  Mi  SiT  (aTCNMAcLV  AM)  Oi  Tl  MN  |Nt  MMlCM 

C  •  NOoTlNtt  AH  i  10  |*I  I  ALLOT  ID 

C  • 

c  . . . . 

c 

CONNON/SVSUH/AI 1000 ),L I AOOl 

c 


C  VIHS10N  I  01  N  T  1  f  1  t  H 

c 

OATA  H0SCNT1  /ilOd'A/ 

iUOl  FALINCH  hOSCHTJ  ,U  ?7H 

c 

iUtfl  VALINCll  UOaSL  #L  I1S0>  »#(  IaOaCL  iUWOIIi 
I  (IDSCFL  #1  I  I90H#  I  l*»COFL  »Ui)OH 

OATA  UO«SL«  UOaCLaIDSCFl#  l^LOFL/1#  l»  l#0/ 

C 


IF  llAD«SL*fcL«!  (CALL  SfALSL* 
iF I  I AO«CLtib*l  ICALL  SITClH 
1M  1  JSCT  L*LO*  IlCALL  HCADSC* 
IF  I  I**C0FL*(L«1  ICALL  SFTDSC* 
Hi TUHN 
LNO 


IASSIGN  (NISI# 03#  SCM£S  I 


IFOHTMAN  LS#NS 

i  k  T  •  FOMTMAK  |v#  VtRSlON  (00 


vuoh out  t  nc  dscht? 

c 

c  . . . . . . . . 

AJTHPH  i  B.S.TOHL iNStfh 

o ati  :  ?«•••/* 

FUHFP81  :  This  HOUT|N£  HANOLrS  OlSCNiTiS  IN  FHAHi  {, 
TNI  FLAGS  UOaSl#  UOaCL#  lOSLF  l*  IHCPH 
A  Hi  SiT  fkTlR'tAl.LV  AM)  Oi  T|  NMNt  HN|CM 
H  0U  T | Nt  §  AM  £  10  Mi  (MiCOTiO 


li*  c 

I  A •  COHHON/SVSTCN/AI  10OO )#C) AOU) 

Id*  C 

Id*  C  ViMSION  1 01 NT  IF  1(1 

l/*  c 

IS*  OATA  KOSCR  T  2  /l*Q*?*/ 

IS*  iUJlwALlNCK  aDSCHT*  «u  7S  )  ) 

?o*  C 

21*  iuul  VALiNCU  IADaSl  •  LII*OM«t  UOaCL  iLHHHi 

H»  I  IIOSCFL  »L  I  ISO)  )»  |  1SC0FL  #LIF/L)I 

si*  c 


I  ASS  l1 


IF  I  I A0«6L*iU*i ICALL 
IF  l|A0«CL*L6*i)CALL 
IF ( |0SCFL*IS*2 )CALL 

if \ |scofl*Cw*2 (Call 

Hi  turn 

LNO 

«*ISl#03# 1S0CS > 


If  OH 
Ik  T  * 


ONTMAN  Ivi  ViSSlON  too 


Hr ADSLH 
S(  Klh 
SCAOSCH 
SITDSCH 
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1  • 
2* 


b* 


ft* 

!»• 

1U* 

II* 

w* 

1 J* 

lb* 
1ft* 
1  /• 
1ft* 
111* 
«o* 
21* 
22* 
2  J » 
2ft* 
2b* 
2ft* 
27* 
2«* 
21. 
10* 
21* 
32* 

IF  In 


FUNCTION  iftUSL«i  ISiT.t INI > 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

c 

C  VtNftlOK  I  01  NT  IF  lift 

c 

OftTft  ft  ISOftC*  /|90»7ft/ 
tWUlyftV.CNCl«NlSL'SCN  #L(  ;yn 
Ol*lN»jON  1M32»#  IUb(  J2» 

tQUl VAlINC€(  t&L* It  »  »L I Ibl » I « ,  ICL* < l  I  »L <  1  M  n * 

1  *L  I  1*1  »  »  .  I  I0S(  |  I  *u»22>  >  »  .  »  ICPl  |  I  ,L  ic/1  || 

CftwL  «C*OSC« 

C 

ISOSCR*  1 P  C  L  IN£  l 

if i i sc r *iu* i i i sotcN«  losiuiNti 

NCTUNN 

4.H0 


AJTNOft  !  0  •  N  •  T  (>*l  t  Nft  ON 
OftTt  l*-ft->* 

FjwPObl  •  Th|B  FUNCTION  NfcftOS  S10«A  DlftLNtUS  IN  ].£. 

P  ft  T  C**  Ab  l  f  LtUlCftTkO  Dis>»  b»UC**tb. 

1*  lb  IsTiNJi.^  *t*N  Cbk.CUTlC'*  Wc*l»«G  '  IM)  Itl1  • 
ft  NOUF*(  N  T  s  : 

1SCT  SftlCTS  PftlCnftHLi  0  UCftt  T  tb  <  0  I  On 

UtU ICftTfcO  D*j»*  bft i T  C*k  S»  I  |  >  • 

LIM  U£f  INtS  twi^CNT  IN  fthMAt 
ISUbCN  IS  NfcTv»N\tO  0  ON  l 


CONHON/StSTlN/fti  10UU  i  , l t  ftco l 
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7 

• 

9 

10 

n 

a 

u 

is 

13 

15 

w 

It 

\9 
20 
2 1 
22 
23 


32 

3% 

33 
35 
3  7 

35 


PROGRAM  GENERAL  HQOCL 

C  SIT  VERSION  NUMBER 

OAT  A  I V EH  SI  ON/ 050£ 7 1/ 

12 • t • 7 A  SYSTEM  COMMON  ENLARGED 

23 1 R • 7 a  01  SC RE  TC  6  INCLUDED 

3*10*74  FRAME  TIME  VARIABLES  TlOlEO  UP 

21*12*74  TOT N  ANO  TOTF  SEPARATED  AND  SOURCE  CLEANED  UP 

5*  2. 7b  short  name  list  specif ieo  to  cut  storage 


control 
mode  ■ 

ENDC 

initial 


RE  alt 


01  HENS  JON, COMMON, OAT  A  OE CL AN at  IONS 

COMMON/STSTEM/A i loooi #l ( aoo  » 

DIMENSION  RjCVALSteoi 
DATA  Nl PASS / 1 / 

equivalence  i «icvals 1 1 ), a <  un  > 

COMMUNICATION  between  • INITIAL*  and  SYSTEM  COMMON 


24 

EQUIVALENCE (VKNIC 

#  AI2U  1  it  1  VKtIC 

•  A  I  2 1 2 

it 

23 

IX  1C 

•  Al E9o 1  l»  IYIC 

• Al 194) 

)• 

24 

IPDIC 

, A 1 22u ) ) a I0U1C 

• Al 221 

it 

27 

(PMlUlC 

, A ( 21 4 1  1# ITMETAOIC 

• Al 219 

1  • 

25 

IFRAMEtj 

• Al  9gi 1  •  •  • 

29 

ioeotor 

, A 1 142) it I0ELT1 

•  Al  91 

) 

EQUIVALENCE (NCPAS9 

•  LI  21  I# IN1PASS 

•  LI  3 

1  • 

U 

II6MH 

•LI  9I> 

#A(2i3  n*  •  • 

•  A  I  1 97  >  > •  •  • 

•  A 1 222  )  I  •  •  • 


REAL  Time  PARAMETERS 
DATA  FRAMETl/60*0/ 
data  TIMLEF Tl/0*0/ 
data  DELT 1/0*050/ 
NUMBER  OF  1C  VALUES 
DATA  NVALS/1%/ 


39 

40 

41 

42 

43 

44 

43 

44 
47 

45 
49 

30 

31 

32 

33 

34 

33 

34 
37 

35 

39 

40 

41 
A2 
53 
44 
A3 

44 
#7 

45 
49 

70 

71 

72 

73 

74 

73 

74 


C 

C  SET  UP  COMMUNICATION,^  ENABLE  VARIABLES  IN  SLI  LABELLEO 
C  COMMON  TO  Bt  ACCESSED 
C 

DIMENSION  Name l 2 1 
DATA  NAME/'MODEL*/ 

NAMELIST  C I  NT  1 ,  UlOI,  JALGl*  INK,  IHY,  HP,  1 V  tN  S I  ON  *  M00£  ,  .  . 
NSl«NVALS#S*#SY#HM,SkiL*SYIC,SHlC«T#TlMLCFTl 
LABL1CALL  PV  200  I  NAME ) 

GO  TO  LABO 

LABBJINPUTU01  > 

LABO : continue 

Call  syscom 

CALL  USERCOM 

c 

C  CHEC*  ACCESS  'ONI  AND  INITIALISE  A 04 
CALL  INTLAD4 

C  READ  IC  CONDITIONS  from  MASTER  FILElf U IaPLACE  IN  SYSTEM  COMMON 
CALL  RDICFILC INVALSfRICVALS I 
USFIG  •  RICVALSI  10) 

1SHR  •  RICVALSI  11  ) 

c  •  input  retained  Changes  • 

CALL  PV 300 
C  INTERPRETER 
c  MMMMM* 

1  TP  •  1 

C  CALL  PviOO  IP  OE  SR  Switch  1  ON 

IFt  ISDSCRIIM  )*EU.1ICALL  wvlOO 

c 

C  STEP  SIZE  If  OR  REAL  T]Ml,  •  FN AM£T 1/ lOOO *0 ) 

C  SLI  CALCULATES  INTEGRATION  STEP  FROM  ClNT/NS* 

C  FOR  USER  ACCESS, SET  C|NT  FROM  DCLt 
CINT1  •  DEL  T 1 

c  initialise  system  functions. 

CALL  SINIT 
C  INTERPRETER 


05027B 


050275 

050275 


38 
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7/ 

C 

79 

l  tp  •  a 

79 

C 

CAbL  TwlOO  IT  DCIA  Switch  1  ON 

•  0 

IT! IS OS C"ll. | t.CQ.l ICALL  wvlOO 

•  1 

C 

•  2 

c 

SET  UP  INITIAL  VALUES  TOR  INTEGRATIONS 

•  3 

C 

B% 

V  VR  NIC  •  VKNIC 

09 

VVREIC  •  VWtIC 

09 

VVWD1C  •  VADIC 

9/ 

STIC  •  «lc 

09 

SYIC  •  Vic 

99 

•Hie  •  Hie 

90 

pwic  •  PDIC*DEUT<JR 

91 

ootc  •  aoic%otaTfR 

92 

RR I C  •  RD1C*OLOTOR 

93 

PPHIC  •  PMlDlC%0f OTOR 

9% 

TTMIC  •  ThctAOIC*P£(jTOR 

99 

PPSIC  •  PS  10 IC  %0E0  TOR 

99 

c 

9/ 

C 

for  control  of  *post-inteuRat ion*  CalCULATI' 

90 

c 

suit  integration  algorithm  iaiqor. default  i 

99 

IF(  UlQ1*EU%%  IHIPiSS  •  2 

loo 

1 T I  IALOI .tu.5 IN1RASS  •  « 

101 

c 

TO  ENSURE  That  SUC"  • R03T.I wICORA Tl ON '  ROUT 

102 

c 

INITIALISATION  pass  through  « 20001' 

103 

NCPASM9 

10% 

c 

109 

c 

INTERPRETER 

109 

c 

107 

ITP  •  Z 

ioa 

c 

call  Rvioo  it  pis*  switch  i  bn 

109 

ITl ISDSCR 1 1. 1 1 .LW. 1 ICALL  WV100 

110 

c 

111 

1VIT  •  0 

112 

in i t  ? : cont i nue 

113 

INIT  w  IN  I T  •  t 

ID 

11* 

11/ 

ll« 

119 

120 
121 
122 
12  J 
12% 
129 
129 
12/ 
12a 

129 

liu 

131 

132 

133 
13% 
139 
139 
13/ 
139 
139 
1«U 
1%1 
1  %2 
1  ♦  3 
1%% 
1%9 
1  %9 
1%/ 
l%a 
i%9 

190 

191 

192 


C  INITIALISE  TWICE  ao  THAT  'UDvn'  etc  SET  up  CORRECTLY 
IF*  INIT.EUM  »  00  TO  I  Nl  T2 

c 

C  INTERPRETER 

C  MMMMM* 

I  TP  w  9 

c  call  pvjoo  if  desk  switch  i  on 

IFllSDSCftl 111  l.LU.l  ICALL  PVlOO 

c 

C  START  CLOCK  for  REAL  TINE 

ir<«ooc»  call  startc 

OYNAMIC 

C 

c  OYNAMIC  »M 

c 

1FI MODE  I  CALL  E*lT8 
DERIVATIVE  L00P1 

c 

c  DERIVATIVE 

c 

c 

c 

C  •  • 

C  •  SYSTEM  FUNCTIONS  • 

C  •  • 

c  COMMON  STATEMENT  TO  allow  CONTROL  PARAMETER  'MOPE ' 

C  TO  0E  ACCC8S60  BY  TMt  DERIVATIVE  SECTION 
Co««ON/*9»iio/fei>t 
PECLARC  LOGICAL  "001 
C  INTEORATION  CONTROL  8T»T£"l"tB 
*»«I»BLl  T  •  0.0 

CINTERVAL  CINTI  •  0.00 
NSTCRS  N$|  •  1 

A  LOOM  I T  Mf.  IALOI  •  S  <  jALUt  •  8 
HIlBUifKAnill.  I>LAUl/rl"ctFT1.1.0#l  I 

rrpcepural 

C  1?-*-?*  STS  It  H  CO  ""ON  (NLAROEO 
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153 

COMMON/ ST ST E"/X 1 10001  it l*UU  1 

15* 

C 

155 

c 

••  INPUT  VXRlXBLES  »• 

15* 

c 

15? 

CBUIvXlENCEiPNIDT  ,X|  gtt  1#  IT  Hfc  TOT 

•  A  |  30)  I,  IPS10T#AI 

3111m* 

15« 

iPOOT  <xi  «»>i»iouot 

•  A  l  «*))#IR00T 

•  Al 

*7  m  •  •  • 

159 

| VKNOOT  • A( 23») !• I V«EDOT 

• A  I  239  1  )$  1 VKUOOT« Ai 

2*Q 1  l«  • 

1*0 

ixlc  >*1 1»9I  i» ixuarn 

•  Al 159)  !•  1  DEL  T 1  •  A  1 

91  1  )••• 

1*1 

IXlHX  f R|IS»I l« IXlMV 

•  At  1*01  1 

1*2 

c 

1*3 

CQUIVAlENCEiNCPASS  «LI  *>>/tNlPASS 

•LI  311 

1** 

c 

1*5 

c 

••  OUTPUT  VXM (xBLES  •• 

1** 

c 

1*? 

EQUIVALENCE |PMIN  •  A  |  / 1  )  t  1  TMC  TAR 

•  A 1  51  1# IPSIR 

•  A  1 

3  I  1  •  •  • 

!*• 

1 P  • Al  2*1)/ (0 

•  Al  271  1 •  IR 

•  A  1 

illlXI 

1*9 

IVXN  /X(  521  1# 1 VKE 

•At  5311* IVKO 

•  At 

»*))».• 

170 

IV  • Al 1*91 1/ 1 t 

<XI  |t>0  1  It  IN 

•  Al 151  )  »••• 

171 

(TlNi  t X)  »•>> 

172 

£N0 

173 

c 

1  7* 

c 

•  • 

179 

c 

•PHt.lNTtQMAL  CALL 9  • 

17* 

c 

»  • 

177 

c 

179 

c 

179 

c 

*IA0  AnO  pRITE  OltCHcTC  LINE* 

190 

CALL  08CRT1 

191 

c 

192 

CALL  SCOUNT 

193 

c 

19* 

PROCEDURAL 

199 

OATa  IF  IRST/ j / 

19* 

IF  1  IF IRST  *N£  « 1 )  00  TO  NOT? 

197 

c 

INSERT  CALLS  Hf  R£  F0R  ROuUneS  TO  B£  EXECUTED  aFTIR  AOC 

BUT 

BEFORE 

199 

c 

•MX  T I  ON  »XND  »MICH  XRE  Nor  PXRT  OF  EVXLUXTION  LOOP  FOR 

DERIVATIVES 

199 

c 

only  execute  these  Routines  once  per  frame 

190 

c 

1)  pre  derivative  s/r  call* 

191 

C 

192 

CELL  CONTROLS 

193 

c 

CALCULATE  MINO 

19* 

CALL  StolND 

190 

c 

19* 

c 

RESET  FLAG 

197 

1 F 1 RST  •  0 

»»• 

NO  TF ICONTINUC 

»»» 

E  NO 

200 

C 

201 

PROCEDURAL!  •  AOCUll 

202 

END 

203 

c 

204 

C 

209 

c 

A 

20* 

c 

•  ahoular  «ot: 

207 

c 

• 

209 

C 

209 

c 

210 

c 

RTMTUDE  ENGLE S 

211 

PPSI  •  1NTEUIPPSIDT.PP81C 1 

212 

TTHET  •  INTEOlTTNOT.TTNlC l 

213 

PPM!  .  INTEOlPPHlOTxPPNIC) 

21* 

C 

219 

PROCEDURES  •PPSI/TTmCT/PPMI  i 

21* 

c 

EQUATE  VARIABLES 

217 

PSIR  •  PPSI 

219 

TNETXR  •  ITMEt 

*l» 

PNU  •  PPHl 

220 

c 

CXLCULXTE  DIRECTION  COSINES 

221 

t 

222 

CALL  BDCOS 

in 

C 

22* 

ENO 

229 

c 

22* 

PP  •  INTCOl PPDOT/PPIC 1 

)8 
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r 


227 

02  •  IfcUGIOCiOflT'OQlC  » 

22* 

*4  •  mTt<i(«KOOT*«R  iC  t 

22* 

C 

230 

PROCEDURAL*  PPSiDT#  TthOT^P^hiuT-PP/WU#*.*  » 

231 

C 

19U»U  v»*l»BLe«  • 

232 

p  •  pp 

233 

Q  •  OQ 

23* 

R  •  RR 

23* 

c 

CALCULATE  ATTJTuOE  «AT£6 

23* 

CALL  SEULEK 

237 

PPSIDT  •  PSIDT 

23# 

TTHOT  •  T  H£  TO  T 

23* 

pphiot  ■  phiot 

2*0 

END 

2*1 

c 

2*2 

c 

2*  J 

c 

•  • 

2** 

c 

*  L I Ni AN  nol ION  • 

20 

c 

•  • 

2*# 

c 

2*7 

c 

2** 

c 

VIL9CJT1ES 

2** 

W*N  •  IKtEQt  VVHH07  #VV*NlC  » 

2*0 

W*E  •  INTtGl  W*tOT#V  VAElC) 

2*1 

vvao  •  JNTtai vv«oot»vvaoici 

2*2 

c 

2*3 

PMOCEOUPAL  «  *v VAN;  V V K£  $  V  VK 0 1 

2** 

c 

laUATC  VARlABLCi 

2** 

VAN  ■  VV#N 

2*6 

VAi  •  V  V *  fc 

2*  7 

V  AO  •  VVAO 

2*« 

c 

CALCULATE  8  00  V  A  ALE  ViL*ClTlt* 

2** 

CALL  SVEL0C1 

2*0 

c 

CALCULATE  ALPHA  A  HI)  BETA  AN  0  ASS0C1ATL0  *  INIS /C  OS  INfcS 

2*1 

CALL  SALEBET(OELTi > 

2*2 

LNO 

2*3 

PROCEDURAL*  HHQOTi SVQOT#  » 

2** 

c 

2** 

C 

EQUATE  VARIABLE* 

2** 

HHOOT  •  • VA  Q 

2*7 

8 TOOT  •  VKE 

2*8 

6 *00 T H  •  #00 T N 

2*9 

END 

27U 

C 

271 

C 

CALCULATE  POSITIONS 

272 

8#  •  HO0lNT(SK0OTH,*AlC#i#v*lR**T » 

273 

8  T  •  H001NT (SVOOT*SVIC» l»u* l«t# T ) 

27* 

hh  •  1nT£G<mhO0T#ShICI 

27* 

c 

27* 

PROCEDURAL 

277 

c 

USER!  ROUTINES 

27* 

c 

CALCULATE  total  FORCES  in  BOW  AXES 

279 

c 

2*0 

CALL  TOTf 

2*1 

c 

2*2 

End 

2*3 

c 

2** 

procedural <PPDOT/Uwoer«rtRD0r  •  pp»gq»rri 

28* 

c 

uoer'S  routine  re  calculate  wents 

28* 

call  TOTN 

08027S 

2*/ 

c 

288 

c 

CALCULATE  ANGULAR  ACCELERATIONS 

289 

call  saccrot 

08027S 

29V 

c 

EQUATE  VARIABLES 

291 

PPOOT  •  PDOT 

0802  7 S 

292 

QQDOT  •  UDOT 

080278 

293 

RROOT  •  ROOT 

080278 

29* 

End 

29* 

c 

29* 

c 

297 

PROCEDURAL*  VVANOT*VvkCOT»vvkuOT») 

298 

c 

CALCULATE  LINEAR  ACCELERATIONS 

299 

CALL  SACCLlN 

30V 

c 

equate  VARIABLES 

301 

VVRNDT  •  VKND0V 

302 

VVREDT  •  VREDOT 

008 


Appendix  C 


OJ  VtfftOOT  •  VKD09T 

0*  1*0 
0*  POOCCOUOftLl  •  SX#SY#HM| 

o*  c 

0/  c  •  • 

0«  C  •PO»T.|*Tia*UL  CALLS* 

o»  c  •  • 

10  c 

11  *NCPASS«NCPASS«| 

it  c  ncpass  •  count  or  nuhoc*  or  passes  executed  »o  fan  in  one  inteooation  step 

13  C  NIPAOO  ■  CONSTANT#BeT  UP  IN  |NI Tl AL#FOH  Ny«»i*  OF  PASSES  UtCO  OV 

I A  C  INTCOOATION  OOUTJNL. 

10  IF|NCPaSS.LT.NJPASS  100  TO  SXPASS 

i*  c  inocot  N£Ot  calls  to  »post»intcooation»  ooutines#to  oc  executeo  only 

1/  C  ONCC  P£0  STEP# 

10  c 

19  C  Floor  PICK  00  TIMC 

to  TIN!  •  T 

tl  C  CQUATc  vaouolco 

tt  R  •  IK 

t 3  V  ■  OY 

t4  M  •  HH 

to  C  CALCULATC  BOGY  AXES  ACCCLCO at  IONO  F  00  NOTION  I 'OUTPUT*  ONLY  I 

tA  CALL  OACCOOO 

if  c 

t«  C  CALCULATE  INOTOUNLNt  0CA01NU0  I  'OUT out  *  onlyi 

to  C  CALCULATE  1LS  I  'OUTPUT*  ONLY) 

30  CALL  OILS 

31  C 

3t  C  LOGIC  ANO  POSITION  FOG  VFA  CONTOOL  I ‘OUTPUT'  ONLY  I 

33  CALL  STY 

34  C  EQUATE  VAGIABLCB 

30  GXIC  •  MIC 

30  I  NX  •  XIHX 

3/  INY  •  XlHY 

30  OKOOTN  •  XPOTN 

39  C  CALCULATE  CAS#NACH#0YNANIC  PNCOOUOE  ETC* 

40  CALL  OVCLOCt 


1  C  CALCULATE  FLIGHT  OATH  A HOLES 

t  CALL  SOATH 

3  C 

4  C  31  OOST  0E0IVAT1VI  0/0  CALLS 

0  C 

A  C  0/0  'OUTSO*  HANDLE!  OUTPUT  FUNCTIONS 

/  C 

5  CALL  OUTSO 

9  C 

0  CALL  SDAC 

1  C 
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3  IFIOST  •  1 
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0  CNO 

A  C 
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•  OACl  A 1 • OACl  7)#DACt  St#OACt  9l#DACilO># 

9  OACl || I#OaC| it  >#UACl 13  I# OACl 14 i# DACl 10  I# 

0  OACl  I A)* OACl  l7)#UACl  IS)#  OACl  19  I#  OACl  tO  I# 

1  OACl tl I# OACl tt )#0AClt3)#0AClt4 )#OACl tB I# 

t  0ACltA>#0AClf7 J#OACltSI#0AClt9l#UAC|30l# 

3  OACl 31  I # OACl 3t ) f  UAC l 33 )# OAC I  34 ) # OAC l 30  )# 

•  OACl 3A  *# OACl 37 l#UACl 3S l#OAC| 39  l#UACl 40  I# 

A  OACl A| )«0AC|4tl#UACl43l#0ACl44l#0ACl4»># 

A  DACl 4A )*  DACl 47 )#  UACl 4t )• ) 

7  ENO 

•  C 

9  C  OOOCCOS  AUTONAT 1 C  HOLD 

U  CALL  OV 700 

I  IFINOOE)  call  TENnnaTEI 'MNCOAFT*  I 

t  END 

3  END 

4  TCONINAL 

b  ENO 

A  END 
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Pt/lIC 

imn» 

•  INIT 

ivueci 

II 

ItM 

Total  ain  AOkkO 

knot* 

ivuecr 

Hi 

vtuHt 

connoncnt  or  TUNOULtNLk 

P  T/HC 

•  ■INI) 

•  Hi  NO 

OTOTtX  CSNNCM  .  W|AL  V  AN  1  A  At  t  A  |A>J>Tt 


ALHNHUTIC  OH  01  1 

lbi«»  hah 

tlf'79 

HAOl 

ILL***? 

PONTOAN 

QUANT  Jt t  OaOCNIHTION 

UNIT! 

CALCuL ATlO 

uuu 

NO* 

NANf 

IN 

IN 

II 

VHD 

VfLOCITT  Of  H|NO*06hN 

P  T/*|C 

•  ■INO 

B vCLOC | 

#41 

VHOL 

H(AN  hINO  COHO  at  Ml  |<»*t  t  oo*n> 

PT/ltC 

•hind 

SH|NO 

ffJ 

VhOLO 

•|N0  VJLOCI TV lOOHNI AT  HlldHt 

M/ltC 

•  HIND 

•  ■l  NO 

1* 

V  •! 

VtLOCltT  Of  W |NO»i Alt 

PT/ltC 

•  HIND 

•VlLOCl 

DO 

VH(L 

Hf  A  N  NINO  LOHH  AT  M(|OhT((A«TI 

PT/»|C 

•  NINO 

•  h|NO 

ID 

VHtLO 

H 1 NO  VC  LOCI  TV  <t*iT|AT  HilQNT 

FT/«|C 

•  |N|T 

•  ■INO 

•  INIT 

tot 

VHNTO 

UaTUN  NINO  •Otto 

KNOTS 

ICP  ILC 

•  INIT 

i» 

VHN 

VfLiC IT  V  Of  «|NO*NOHTH 

PT/ltC 

•  HIND 

•vtLiCl 

fff 

VHNL 

H(AN  Hi  NO  COHO  at  NtlNHTtNTNI 

P  T/UC 

•hind 

•  hind 

Hi 

VHNLO 

PT/ltC 

•  INIT 

•hIND 
•  INIT 

in 

N 

LO 

ICP  ILC 

•  INIT 
SACCBOO 

•0 

Hi 

Vfi.  CONN  h%W  Ti  A|«»10QT  AK|« 

PT/«tC 

•VtwOCl 

•VCLOCt 

•alpoit 
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us 


I 


Ik'NlMTIC  9IIHN 

ili*tNT  PONfkAN  OUANTJTV 

NO «  K*Nt 

i*f  h»)U 

!J«  NTuM 

l«t  * 

#NP  111 

#N©  It 

#•4  kj 

t»l  M 

<•1  k» 

l*t  kcvl 

•»  kci»*t|k 

!»•  *UOT* 

!••  kpt  m 

Ibft  kofv 


»>»U"  iOt'Mi  •  kt  a 
lui*»  NAk  I 

Pt»c«|PTie* 

NMf  pi  TUNSOLtNCt 
LflHkCNINT  Of  KitfftOttfcCt 

k  peime>i  of  c.u. 

tppm  IKAT  |  %  Of  SLIP  PALL  |AV|) 
COONOlNkUb  OF  a  "tUN  <A<#| 
COOkO*  Of  ACCCltMOHt  UH  4*01 
koewo*  Q*  Accciikonu*  iav«i 
COONO*  Of  ACCIC1M0H1  UN  iklOl 
k  Ci6«  l©laT10n  In  aialmaM 
k  LOCATION  OF  Nlf  P  T 
N  AT  1  Of  k*i|NU  Of  A/C'b  A  fOt,N 

®LI  *  k«u 
k  Vltdtlt?  fOM  Tv 


UfUk  lM*tt 

\% 

usto 

IN 

•  NINO 

•  NJNO 

•  U 

lilt 

•Acctoo 

•ACCOOQ 

•accooo 

•  AC  CO 00 

•ACCOOO 

NlN|f 
NlMt 

Oik  Mill 

•  TV 

•  Tv 


uMTv  CAtkutAUO 

IN 


M  /  V(  C 

UNA  k 

M/UC 

•  ■JNO 

P  T 

01 k • T lVi 

P  T 

until 

f  1 

U»4.k 

f  T 

USlk 

P  T 

UNI  k 

P  ? 

UNI  k 

ICP  Hi 

UNik 

M/»|C 

•  TV 

f  1 '  Si  c 

•  TV 

f  T  /  si  C 

•  TV 

•  VtH*  CONNOS  •  k| At  V  AN  1  A  Hi ! •  JJ.J.If 


okotk 

lows  nan 

1  J*  *  f 

PAll| 

NO  # 

NANI 

QUANT  ITT  l>l»C«l»IIPN 

UNITS 

CALCUL AT(0 
IN 

U«(0 

IN 

kOt 

IlklOl  »k9*l  0A1UIN 

*  T 

•  UP 

•  ILS 

Iff 

kIC 

ini t i ik  leuiTieN 

P  T 

ICP  IL1 

•  INIT 

isinii 

•  Tv 

IM 

*IC» 

P  T 

•  TV 

•  TV 

»*» 

•  |MR 

CANTNOU  VAN  POk  k  tMtOMftTlAN 

•  TV 

OlkMIVi 

wo 

klHT 

k  Is  TN  9k  V  IM  1 9*  T  INTltflltlOK 

•  TV 

0|*>Tl«l 

too 

■  KWUI 

HITkia  In  >»*»k'  »»»» 

icriki 

INITIAL 

toi 

MtNk 

iimn  in  ‘nil*  ’■*" 

|CP ILC 

INITIAL 

•  Ik 

NfWfNT  91  |N€»T|l.l9kk 

•  CUO  P T* 

ult  k 

•  IN  1 1 

At 

«IT 

ItMM  0*  islNT  |IiI|TCn 

•cut  PT* 

u*v« 

»1N>T 

to 

111 

k On l N T  OP  |NlkT  |Af  VAN 

•CUO  FT# 

utlk 

•  INIT 

tl 

■  111 

PIONINT  OP  |N|kl  |A#PkOOUCT 

•LUO  PT# 

UM 

•  INIT 

III 

kkfPPt 

CONSTANT 

•  INIT 

•  nINO 

•  TV 

•  INIT 
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cbxxsx  -  mi  «uimn 


Ai*»A04T|C  0*04* 

|«l«l  MA 

IS# 'TO 

*A0& 

Ikinkf 

NO* 

rMfiK 

•MNC 

OUANT l V T  0tOC*t*t|O* 

UNIT* 

CALCULATfP 

IM 

uuo 

i* 

III 

uiw 

*OOITION  0»  OtA*  »1N* 

M 

•  Its 

oite 

»4| 

•kkl*» 

TffAL  NlK|MrfMU,»IOT  AfeLO 

lo»kt 

uot* 

•accoot 

!>• 

MC 

C9CUIIM  0«I9!N  III 

9  1 

oite 

tit* 

|M 

unu 

talk  »»»  lain  x.M.t. 

9  1 

uot* 

To 

MACK 

x«Ch 

ovitecf 

•vtLOC* 

i 

MAN 

Otuo 

unit 

UNIT 

III 

IWTfV 

TOTAL  NONtNT,M!C*«*oeO>  lid 

tO»fT 

uot* 

SACCMOT 

IK 

•NNTOT 

t*Tak  iKiini  ml 

lo»kt 

uot* 

SaCCOOT 

iso 

M 

»  IK  0»  »lk»t  «kk  TB  RlfCii 

9  T 

uot* 

•tv 

no 

MCI 

•tec  a#  Near  *kt  to  cu 

•  TV 

•  TV 

SACCOOO 

III 

•  TV 

9  T 

•  TV 

•  TV 

too 

T 

t  position  OK  c*o* 

9  T 

OCN'TIVC 

•  TV 

•  us 

ill 

u 

coeooitiAf m  ok  otie  OAtt  iavii 

M 

use* 

SACCOOD 

OVOTIN  CONNOV  •  N(AL  VANlAOtCS 

1 9  •  3  •  Tl 

AtPHAOtMC  ONOt* 

losoa  nan 

It# ' TO 

*AQt 

NO* 

K0*T*A* 

NAN| 

OUANTITT 

Ot»C*|OT|ON 

UNIT* 

CAtCut ATIO 
IN 

usto 

IN 

*» » 

vt 

CBBIDIXaTaa  or  li  xk  TIN  lAlfl 

Kf 

UtCN 

SACCOOO 

IA« 

Tl 

c»e*o. 

OA  ACCCLLNONt It*  IAMJI 

kt 

U*l* 

SACCOOO 

to/ 

TO 

Cxeao. 

B*  aCCikiMBntTia  iat.1 

kt 

UStN 

SACCOOD 

•  TO 

TO 

Cbbbo. 

0>  »CClkk*6xtUR  i»2Si 

T  T 

ust* 

•accooo 

♦oi 

TAOClII 

AMOAT  0 

>  uticakio  inpuis 

•  AOC 

use* 

AAA 

taoCiaa  i 

•  Mat  a'  mtcakto  iNBui* 

SAOC 

ust* 

100 

TOT  V 

t  vctociTv  roe  tv 

KT/ttC 

•  TV 

too 

TIC 

initial 

*00  1  T  ION 

9  T 

ICKItf 

INITIAL 

III 

TtOC 

toe  alio 

I*  0*1 01 N  (VI 

KT 

•  Its 

tlk* 

III 

TTV 

T  fOOIT 

ft 

•  TV 

•  TV 

too 

ii 

CBBBOINara*  er  kip  aatk  iat|i 

KT 

ust* 

•accooo 

tot 

ii 

coeooiNATto  or  •  nc tin  iuii 

KT 

Util 

•ACCOOO 

too 

21 

C 80*04 

OK  ACCtttNONt  UN  1  AM)  1 

KT 

ust* 

SACCOOO 

too 

20 

ceeao. 

OK  ACCf  Lt NONC TIN  IAVAI 

9  T 

use* 

SACCOOD 
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•ton*  common  •  Hi  At 

VAHUOtlS 

IB*)* 

AlPMAtfiriC 

ONOiN 

tout  MAH 

19* 

«  ?f 

AiHl 

|U"|NT 

MO* 

reNtoi* 

nan* 

OUANt  |tV 

OliC"|*tlO* 

UNIT! 

CAtCot At|0 
IN 

WSIO 

IN 

*M 

i» 

c«e»o> 

0»  ACCtt|«0*lU«  (Alftl 

f  t 

UtlH 

oiccooe 

90 

ica 

1 

CtOa 

location 

UOIH 

0  IN  1 1 

IOO 

1CIM’ 

1 

t«CA 

U«4H 

•  INI? 

no 

IP 

1 

IOC 

9*  H|  tOt  Hit  to  Mifta 

r  t 

UOIH 

otv 

tot 

ip ta 

1 

tec 

or  *»i  tot  nit  to  co 

•  tv 

•  Tv 

ktCCISV 
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Appendix  E 

INDEX  TO  SYSTEM  COMMON,  INTEGER  VARIABLES 
E.l  Numeric  order 

'  <.•**©*  •  iHuuk©  vaniaolki 


hU*t*U 

0*0  IN 

lotto  nan 

lie • 79 

OAQt  l 

0U4NTITT  U*»C«l»TlON 

UN|T» 

CALCULATED 

U9I0 

NO* 

NAN| 

IN 

IN 

1 

MOA 

bf  liC  T  T4N14U3N 

u*c* 

•viloco 

SIN1T 

• 

• 

NCOAOO 

Count  04  niiu  tn*b 

INITIAL 

OCO'TlYC 

•COUNT 

1 

NIOAOO 

count  or  Him  to  sun  intco 

INITIAL 

OtNMlTS 

• 

A 

UIM 

•CLlCTt  U»  QCANittC  MLt»L«» 

INITIAL 

til* 

9 

l*M« 

»|N0  *T*4  4  4  |T*4*i  I  Oi  i  IISKHI 

INITIAL 

•  «1nO 

SINIT 

0 

IOYO 

UNIT  NO  T  0«  STS  UN  CN4NUCS 

uoio 

•YOCON 

7 

NTYO 

Ty  OlcT  OalCCT 

0910 

•  TV 

• 

JX  ONN 

In|  Tt  41 1*41  ION  CONOitU  »L*0 

•COUNT 

•  U9 

•  TV 

0*1  NO 

9 

NUNO 

U9I0 

SOAC 

10 

NUNV 

NO. or  V  4* IU4CS  Oi  SiT  TNON  14* 

utio 

•  OAC 

tl 

NOACt 1 1 

reiNTC*  Tg  V4NI40H  TON  D4C  | 

utt* 

SOAC 

••••••••  N0H  NIUIM  (II^CNT  NUNKM  •«••• 


*u"ln  ic  oouco 

ILiNiNT  FONTOAN  QUANTITY 

NO  •  NANf 


•  V  ST IH  CONNOS  •  ISHOtO  VANlAOLCS  |4.0J.?« 
lotto  NAN  I 1* 1 79 

OaOCNINMON  unit*  CACCUCATiO 

IN 


OAQi  2 

utco 

IN 


••••••«  NOT!  NIOOINQ  (CfNCNT  NUNtfOO  ••••• 


70 

N04CI44 1 

nointc*  To  vi 

IOIAOLC 

ON  Oal  *0 

UU« 

II 

MAOC 

VCOtlON  III  - 

OOUTINl 

•AOL 

•  AOC 

70 

•  00  AC 

«(NS|ON  III  - 

OOUTINC 

•OAC 

•  OAC 

77 

■  MCNTI 

VIMION  III  • 

ooutinc 

OICNTI 

0«C0T| 

70 

ooocort 

v(N«|  On  111  • 

OOUTINl 

OtCNT  i 

OOCOTt 

79 

Kiooocn 

VCOtlON  la  • 

OOUTINC 

ItUtCN 

| to ICO 

•0 

•  1 

NAOCl!) 

4004T  or  4DC 

CNANNCL 

N0N0LN9 

uoco 

NO  NO  (00  MM* 

too 

NAOCl 001 

4004T  or  4IIC 

CNANNCL 

NUNOCNt 

utco 

tot 

100 

107 

•  OAC 


SAOC 


SAOC 
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V 


1  19 


Stmt"  10HH8N  -  IMtlit*  vtHUBLIS  l**O3>/0 

NUHLHlL  OKDIH 

10130  nan  13# • /f 

*AG 

lUi»»lNT 

FO*T*A 

N  QUANTITY  INSCRIPTION  UN  1 T 6 

ClLtUL  »UP 

US£0 

NO. 

NAH£ 

IN 

IN 

!•* 

LS1L6 

C  On  T  h  0c  *  IKCUTICI  Of  tlL. 

use* 

OILS 

1*9 

100 

l AD*§L 

C  «N  TBBL  i  tltCUTIBK  0*  HC*PkL* 

USt  * 

0  SC  *  T  l 
09C*T  3 

131 

1  8L  N  |  I  1 

AO*  6CN6L  LlNC  1 

REAOSLR 

U  S£  * 

HI  86 lNQ 

CLCHCNT 

NUH»CR6  ••••• 

1 A  A 

ISLHl 1.) 

*0*  6CNSI  L 1 N (  18 

*£ A06L* 

US£  * 

18/ 

NO 

8CALC  FACTO*  I**  NQU6T6 

use* 

S-INO 

188 

L  St  ID 

LOnTNOLS  N*NDOH  NUH6|N  SitO 

use* 

SrtINO 

189 

NNUN 

CURRENT  *UN  NO.  »«QH  HAlLOOK 

SCOUNT 

1/0 

1 A08CL 

CONTROL*  UCCUTION  Of  bETCLR 

use* 

DSC  *T 1 
DSC  *  T  3 

1/1 

1  Cl  r i l > 

* D*  CONSUL  LINC  1 

use* 

6CTCL* 

H1681N0 

CLCHCNT 

188 

ICLNI  18  1 

»P«  CBNtHoL  LINt  1* 

use* 

8CTCLN 

SVST£H  COHHON  •  INTEOCR  VAN  I  ABLE S  16*03. 7» 


NUH£NIC 

ONOCR 

10:30  has  1 3# ' 79 

CLlHkNT 

FORTRAN  QUANTITY 

OesCHlPTlON  UNITS 

calculated 

NO# 

NAH£ 

in 

16/ 

186 

189 

190 

106CFL 

CONTROLS  execution  Of  NLAOSCR 

use* 

191 

IP(  1  1 

PATChAOLC  UISCRLTE  1  TO  61QHA 

REAOSCR 

>•••••  NOT  C  Nl.SINO  CLCHCNT  NUnBCR 

333 

1*133) 

H*TCH»»LC  ulSCRtIC  32  TB  klQHt 

R£*06CR 

hj 

1 08 1 1  I 

OfS«  SnITcn  1  TB  Slknt 

REAOSCR 

•  ••••#••  NOTE  HI  86 INQ  CLEHENT  HOI 


388 

106133) 

DESK  8b  ITCH  33  TO  SIGH* 

333 

KSACC800 

VERSION 

ID  •  ROUTINE 

SACC800 

338 

R8ACCL1N 

VERSION 

ID  -  ROUTINE 

SACCLlN 

it? 

R8ACCR0T 

VERSION 

Id  -  routine 

6ACC  NOT 

REAOSCR 

saccbod 
s  ac  Cl  l  n 
9ACCR0T 


use  0 

IN 


DSC«Ti 

0SC*T3 

I  SO SCR 


ISOSCN 

initial 

ISDSCH 


ISOSCH 


•8 
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•  >m»  ta»nc*  •  iKTkum  utmiim  i«.or. ;> 


nunc*  1C 

0«OC" 

!0;*o  **•  ii» 

•  AOC 

U(Kk^ 

•0* 

QUANTITY 

Oi®C* !• T l ON 

UNITS 

CALCULAT'D 

IN 

useo 

IN 

*•• 

RIU7NT 

▼  (••1  ON 

lu 

•OUT |N( 

•■(.fact 

SAtFOCT 

»» 

■  SCOuNT 

▼(••ION 

10 

•OUT |N( 

•C Oust 

•COUNT 

*•  0 

■  10  CO  9 

▼  (Oil  ON 

10 

•OUT 1N( 

tocos 

cocos 

»*4 

Rl(UUN 

▼  COO I  ON 

10 

•our  Inc 

•cut  IN 

scute • 

fftf 

MILO 

▼(••ION 

10 

•OUT  INC 

out 

sits 

Iftl 

•  UNIT 

vtutlSx 

10 

•OUT IN( 

•  l N|  T 

•  IN  I T 

2l« 

■  ••AT* 

▼cooiOn 

10 

•OUT  INC 

••■T« 

snath 

Iftf 

■  •TV 

▼(••JON 

lo 

•OUTINi 

•TV 

STV 

MlUKl 

▼(••ION 

10 

•OUT  INC 

•▼CCOCl 

svcteci 

*•? 

MVllKI 

▼(••ION 

10 

•OUT 1N( 

BVCLOCr 

tvCLOC? 

!»• 

■  ••  INO 

▼(••ION 

lo 

•OUT  INC 

•  NINO 

•  *  I  NO 

!•» 

■  •▼•COT' 

▼  (••JON 

10 

•OUT  INC 

ST  SC ON 

STSCON 

*70 

IOC  Oft 

controls  ■■(Cut ion  of 

SCTOSCR 

use  • 

OSCNTl 

DSCNT2 

*71 

icon  1 

CHtNQC 

0«tN 

I  OUT  OF 

•  IONA 

uses 

6CT0SCN 

•TOT (N  CON NON 

•  INTCOCN  VARlAOtCS  1 A  *03* 70 

NUNCNIC 

ONOCN 

io:oo 

NAN  13* • 70 

•  AOC 

CtCNCNT 

FORTRAN  QUANT  ITT  OiSCRJPTION 

UNITS  CAtCUtATCP 

use  d 

NO* 

NANC 

IN 

IN 

*0* 

icon*  1 

ChaNOC  OVCN 

)*OUT  OF  S1QNA 

use  • 

SCTOtCN 

JO* 

AO  A 

*00 

tlOACll » 

OaC  OOuNCc 

useo 

•  DACN 

TC  N1S0INQ  (LCNCNT 

NUN0(00  ••••♦ 

*5* 

CIDACIAO) 

DaC  tOUNCe 

•OINTCNOaL  00 Pj 

U8CN 

SOaCN 

11) 

LfOACll » 

OAC  SOuNCc 

USCR 

SDaCN 

▼00 

LfOAClAOl 

Otc  cauxct 

•0 I N  TC •• *  t  OOP 2 

USC  N 

SOaCN 
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E .2  Alphabetic  order 


•  ▼fit"  COHHC*  •  IhttUfN  VAMlABLtS 


AlRMAttt TIC 

:  9B01* 

10  :bJ  MAM  |)| • 79 

**AOC 

ILAHaNT 

QUANTITY  OikCBIBTlON  UNITS 

CALCULATED 

usto 

no. 

N«H 

IN 

IN 

WU 

l aoaCl 

controls  baCCution  o#  Bit  cl  ii 

USER 

OSCRTI 

OSCRTff 

1*0 

l ADaSL 

CONTROLS  tiCCUTION  Of  NEAOBlR 

UBtR 

08CRT) 

DSC** 

in 

ICLRl  1 1 

A0«  CONTROL  LIH  1 

UBIR 

bctclr 

IBB 

ICLRl  1*  1 

AOa  CONTROL  LINE  Ik 

UBtR 

sctclr 

tn 

ICO 1 1  I 

Cnanoc  o%an  i  Out  Of  sigma 

USER 

BtTOBCR 

JQi 

I  CO  Off  1 

C manoi  oy i m  iff o ut  of  siGMt 

U8t*» 

setoscr 

ill 

i Pi t |  l 

UcSR  BbITcm  1  TO  S 1  OH A 

RtAOSCR 

INITIAL 

1S0SCR 

**• 

1 OS  Off » 

DESK  BN  IT cm  Iff  TO  • Iona 

MtAOBCR 

I  SO BCR 

isa 

l  JSC*  L 

COnTRX*  laCCuTION  Of  heaoscr 

UBtR 

OSCRTI 

OSCRTff 

• 

lCSflU 

•ELECTS  ill  BCAM<B(i  AlLBfLO) 

INITIAL 

BUS 

1*1 

1R|  |  1 

PaTChAOLI  UlSCRLTE  1  TO  BiuMA 

RtAOSCR 

I  SO  SCR 

Hi 

I  *  Off  1 

PATCMAQlI  ultCNiTC  Jff  TO  blQHA 

RtAOSCR 

ISO SCR 

BTBTtM  COMMON  •  iNTtUl 

R  VARIABLES  Ib.OJ./S 

ALPHABETIC  OR 01 R 

10  :so  MAR  |9#  ’ 

79 

RAGE 

tLt*tNT 

NO* 

fORTRAN 

NAME 

QUANTITY  DtSCHlRTlON 

UNITS 

CalCulateo 

IN 

U  St  D 

IN 

170 

IRC  Of  L 

CONTROLS  tlECUTlON  Of  BETUBCR 

USER 

OSCRTI 

OSCRTff 

B 

IBhR 

NINO  SMEAR  SMARttSCE  AlbHMl 

INITIAL 

S  M  1  NO 
SINIT 

l»! 

lit HI t 1 

B0%  SENSE  LINE  1 

REAOSLR 

USER 

IBB 

1SLRI IB  1 

*0*  SENSE  l!*E  1* 

REAOSLR 

USER 

B 

lots 

unit  no  Ron  system  changes 

USER 

SySCOm 

• 

J  jC  OMR 

Initialisation  conrlete  fLAo 

SCOuNT 

SlLS 

STY 

S  A  1  NO 

7} 

ROBCRTI 
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DESCRIPTION  OF  PROGRAM  COMMLIST 

F . I  Introduction 

The  program  COMMLIST  produces  an  index  to  the  variable  names  used  in  a  large 
labelled  COMMON  area.  A  data  record  is  prepared  on  cards  for  each  variable  name, 
giving  a  brief  definition  of  the  variable,  its  units  and  where  it  is  calculated 
and  used.  A  collection  of  such  cards  (which  may  be  held  on  a  disc  file)  may  be 
sorted  according  to  a  number  of  criteria,  of  which  the  most  important  are  alpha¬ 
betic  order  of  name  and  numeric  order  of  element  number. 

For  simulator  purposes,  two  large  COMMON  areas  are  in  current  use,  one 
defining  system  variables  which  are  the  same  from  simulation  to  simulation,  eg 
height,  speed  and  the  other  defining  user  variables  which  are  entirely  specific 


to  the  particular  aircraft  being  simulated. 

F.2  Input  data 

F.2. 1  Main  data 

record 

Column 

Purpose 

Length 

1  -  4 

Element  number 

4  digits 

5-16 

Variable  name 

12  characters 

18  -  47 

Description 

30  characters 

49  -  58 

Units 

10  characters 

60  -  67 

Name  of  routine  in  which  variable  is  calculated 

8  characters 

69  -  76 

Name  of  routine  in  which  variable  is  used 

8  characters 

78 

Continuation  indicator 

1  character 

blank  =  no  continuation 

C  =  continuation  . 

In  general 

,  no  field  delimiters  are  required.  However, 

it 

may  help  human 

reading  of  data 

cards  to  insert  a  comma  between  fields.  This 

is 

allowed,  but  is 

optional.  The  element  number  should  be  right  justified  in  its  field.  A  minimum 
of  2  data  cards  must  be  input. 

F.2.2  Continuation  cards 

If  the  continuation  indicator  is  set  (C  in  column  78),  the  next  card  is  read 
for  additional  ’used  in'  subroutine  names.  The  format  of  this  continuation  card 
is  up  to  8  fields  of  8  characters,  separated  by  commas,  a  blank  field  signifying 
the  end  of  the  list  of  subroutine  names  and  a  field  containing  Cbbbbbbb  signi¬ 
fying  continuation  on  the  next  card. 
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F. 3  Complete  job  deck 

F.3.1  Original  data  on  cards 

! JOB  SIM, BARRY  LIST  COMMON  VARIABLES 

! ATTEND 

.'COMMLIST 

Title  card 
Selection  card 

.  Data  cards 

!EOD 
!FIN 

F.3.2  Original  data  on  disc  file 

It  is  assumed  data  cards  have  been  copied  to  a  file  by  :COPY 
! JOB  SIM, BARRY  LIST  COMMON  VARIABLES 
! ATTEND 

! ASSIGN  (F: 100,DC,RCOMMON) 

! COMMLIST 

Title  card 
Selection  card 

Additional  data  cards  (optional) 

!EOD 

!FIN 

F.3.3  Remarks 

Note,  in  example  section  F.3.2  above,  that  although  the  main  input  file  is 
taken  to  be  a  disc  file,  data  can  still  be  read  from  the  card  reader  in  the  same 
run.  This  enables  additional  data  cards  to  be  added  to  and  merged  with  the  main 

file,  thus  avoiding  the  need  to  recreate  the  main  file  for  every  small  change  in 

* 

,  data. 

Note  also  that  if  original  data  is  on  cards  only,  as  in  example  section  F.3.1, 
.  device  100  does  not  need  to  be  assigned  to  the  card  reader. 

!E0D  is  essential. 

F.3.4  Title  card 

The  80  characters  on  this  card  are  printed  out  as  a  heading  on  each  page 
of  the  output  listing  exactly  as  they  appear  on  the  card.  This  enables  the  user 
18  to  identify  the  nature  of  the  information  being  listed.  If  the  contents  of  the 

title  card  are  centred  on  the  card,  they- will  be  centred  on  the  listing  page. 
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F .  3 . 5  Selection  of  listing  options 

Options  are  specified  on  a  selection  card,  containing  up  to  20  words  of 
4  characters.  The  options  are 


Abbb 

Dbbb 

Nbbb 

CbbbSSSSSSSS 

UbbbSSSSSSSS 

Xbbb 


Alphabetic  order  of  variable  name 
Data  listed  as  stored 

Numeric  order  according  to  element  number 
Lists  all  variables  calculated  in  subroutine  SSSSSSSS 
Lists  all  variables  used  in  subroutine  SSSSSSSS 
Read  next  card  for  additional  instructions. 


The  blanks  (b  =  blank)  are  essential.  There  are  no  field  separators. 

As  many  options  as  desired  may  be  specified.  If  they  will  not  fit  on  one 
selection  card,  continuation  is  indicated  by  the  X  option. 

F.4  Output  listings 

The  output  listing  reproduces  the  input  information  spaced  out  for  better 
legibility.  An  extra  column  labelled  'quantity'  is  left  blank  for  the  insertion 
of  mathematical  symbols  associated  with  the  variable. 

For  the  alphabetical  option,  all  blank  records  ( ie  those  with  only  an  ele¬ 
ment  number)  are  ignored  and  not  printed. 

For  the  numeric  option,  a  warning  message  is  output  when  an  unused  element 
number  is  encountered. 

For  options  C  and  U  the  'calculated  in'  or  'used  in'  columns  are  omitted, 
as  appropriate. 

In  the  early  stages  of  building  a  COMMON  list,  it  is  worth  including  data 
cards  containing  only  an  element  number.  These  will  be  listed  by  the  numeric 
option  and  enable  details  of  new  variables  to  be  written  in  on  the  listing. 

Appendices  D,  E,  G  and  H  were  produced  by  the  COMMLIST  program. 

F.5  Arrays 

It  is  recommended  that  the  first  and  last  elements  of  an  array  are  included 
as  data  items,  in  order  to  make  it  clear  in  the  listings  which  items  in  the  com¬ 
plete  list  are  genuinely  unused.  COMMLIST  will  handle  arrays  of  more  than  one 
dimension,  provided  that  an  element,  eg  X(12,10)  can  be  fitted  into  the  allowable 
12  characters.  If  this  recommendation  is  adopted,  the  user  will  find  the  warning 
message  'NOTE  MISSING  ELEMENT  NUMBERS'  interspersed  between  the  first  and  last 
elements  in  the  numeric  order  listing. 
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F.b  Operational 

considerations 

F.b.l  Devices  used 

DCB 

Device 

Purpose 

J00 

defaults  to  card  reader 

data  file 

101 

console  TTY 

opening  message  (HCOMMLIST) 

105 

card  reader 

title  card 

105 

card  reader 

selection  (control)  cards 

105 

card  reader 

extra  data 

108 

line  printer 

all  listings 

F.b. 2  Messages 

F.b. 2.1  H COMMLIST 

This  message  on  the  console  teletype  announces  the  start  of  execution  of 
COMMLIST. 

F.6.2.2  ***  NOTH  MISSING  ELEMENT  NUMBERS  *** 

This  is  a  warning  message  to  draw  attention  to  gaps  in  the  sequence  of 
element  numbers. 

F.b.  2. 3  COMMLIST  ABORTED  -  TOO  MANY  INPUT  RECORDS 

This  message  is  output  on  the  teletype  when  too  many  input  records  have 
been  read.  It  is  repeated  on  device  108,  supplemented  by  the  last  record  read. 

F. 7  Limits  on  data 

This  program  uses  many  (23)  arrays  to  hold  the  input  information.  These 
arrays  are  dimensioned  at  present  to  cater  for  500  input  data  cards.  However, 
since  the  DO  loops  are  in  the  range  1  to  1000,  the  number  of  input  data  cards 
could  be  increased  to  1000  by  altering  the  DIMENSION  statements  in  the  program 
and  by  increasing  JMAX,  the  limit  on  the  number  of  records  read.  Any  increase 
will  greatly  enlarge  the  core  store  needed  to  run  the  program,  so  should  not  be 
made  unless  really  necessary.  This  present  limit  of  500  applies,  of  course,  to 
individual  variable  names,  so  the  program  will  cope  with  a  user  array  of  2000  or 
more  locations  if  many  constituents  are  themselves  large  arrays,  for  which  it  is 
adequate  to  include  only  the  first  and  last  elements,  as  already  discussed. 

F .8  Location  of  system  COMMON 

The  raw  data  for  system  COMMON  variables  are  held  in  files  on  the  DC  area 
of  the  disc.  Real  variables  are  in  file  DC,  RCOMMON  and  integer  variables  in  fil 
DC,  ICOMMON.  An  index  listing  may  therefore  be  made  at  any  time  using  the  techni 
que  of  section  F.3.2,  with  the  job  card 
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IASSIGN  (F : 1 00 ,DC , RCOMMON ) 
or 

IASSIGN  (F: 1 00, DC, I COMMON) 

included  as  required. 
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INDEX  TO  SYSTEM  DATA  TO  BE  PROVIDED  BY  THE  USER,  REAL  VARIABLES 
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SEE 

LI0ACI4BI 

o*c  sounCe  nointessyloopi 

soacn 

HI 

LfOAClII 

Sac  sounce  sointessyloone 

SDACN 

too 

•BTOM  0 

LtDAC (4Q| 

OAC  sounce  NOINTESSyLOONE 

SOACN 

I  MB 


•8 


LIST  OF  SYMBOLS 


a 

a  f  a  f  a 
x  y  z 

b 

c 

Cref 
Cl |  etc 


F  F  F 
TN’  TE’  TD 

ftx,fty,ftz 

g 

h 

I  ,1  ,1  ,1 
x  y  z  zx 

L,M,N 

M 

^1  *  ^2*  ^3,ml ,ni •  etc 
m 

P 

P.q.r 

R 

xp 

S.  . 

S 

SW 

T 

u,v,w 

V 

V., 


VVVD 

V  V  V 
KN’KE*  KD 

V  V  V 
WN’  WE*  WD 

VWNL0,VWEL0,VWDL0 

VWKT0 

V  V  V 
WNL,  WEL’  WDL 

W 

X 

x.y.z 


speed  of  sound 

accelerations  along  x,  y,  z  body  axes 

reference  wing  span 

mean  chord 

reference  chord 

inertia  constants 

shear  factor 

applied  forces  in  earth  axes 
applied  forces  in  body  axes 
acceleration  due  to  gravity 
height  of  centre  of  gravity 
moments  of  inertia 
total  moments 
Mach  number 
direction  cosines 
aircraft  mass 
atmospheric  pressure 
angular  velocity  components 

vertical  distance  of  pilot's  eye  above  centre  of  gravity 

horizontal  distance  of  pilot's  eye  from  centre  of  gravity 

direction  cosines  (S^  =  etc) 

transformation  matrix 

reference  wing  area 

atmospheric  temperature 

velocity  components 

equivalent  airspeed 

ground  speed 

true  airspeed 

components  of  airspeed  in  earth  axes  (north,  east,  down) 
components  of  velocity  relative  to  earth 
components  of  total  wind  velocity  relative  to  earth 
components  of  datum  mean  wind  speed  in  earth  axes 
datum  mean  wind  speed 

components  of  mean  wind  speed  at  height 
aircraft  weight 
force  along  x-axis 
positional  coordinates 


LIST  OF  SYMBOLS  (concluded) 


X  ,2 

P  P 


X  ,2 
peg  peg 


Ax  ,Az 
eg  eg 


a 

0 

Y 

e 

p 

a 


X 

* 


Suffices 

a 

B 

eg 

G 

L 

P 

r 

SIG 

SL 

T 

TV 

a 


coordinates  of  pilot's  eye  point,  relative  to  reference 
centre  of  gravity 

coordinates  of  pilot's  eye  point,  relative  to  actual 
centre  of  gravity 

centre  of  gravity  displacements  from  reference  position 

angle  of  attack 
angle  of  sideslip 
climb  angle 
pitch  attitude 
air  density 

atmospheric  density  ratio 
bank  angle 
track  angle 
heading  angle 
wind  direction 


absolute 
body  axes 

at,  or  of,  the  centre  of  gravity 

gU8tS 

at  arbitrary  location 

pilot 

ratio 

rms 

sea  level 
turbulence 

television  (for  outside  world  display) 
ambient 


A  dot  over  a  variable  denotes  differentiation  with  respect  to  time. 
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Fig  1  Block  diagram  of  translational  equations 
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Block  diagram  of  rotational  equations 
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Fig  3a&b  Alternative  structures  for  initialisation 
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Fig  6  Information  flow  -  utilities 
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Fig  7  SL1  modal  program  structure 
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* 
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Fig  A3  Flow  diagram  for  ‘initial  flight  path'  phaaa  of  SINIT 
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Fig  A7  Flow  diagram  of  TV  bait  logic  in  STV  subroutine 
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